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Abstract 

It is suggested a metric ansatz describing matter in an expanding universe, hence interpolating 
between the Schwarzschild metric near central bodies of mass M and the Friedman-Lemaftre- 
Robcrtson- Walker metric for large radial coordinate, given by 

ds 2 = Zc 2 dt 2 - ~ (dn - ^-p- + 5 cdtj -rfdCl , 

where Z = 1 — 2GM/ (c 2 ri), G is the Newton constant, c is the speed of light, H = H(t) 
is the time-dependent Hubble rate, c?fi = d6 2 + sin 2 6 dip 2 is the solid angle element and we 
are employing Schwarzschild expanding coordinates n (also known as physical coordinates for 
expanding space-time). For constant exponent a — it is retrieved the isotropic McVittic 
metric and for a = 1 it is retrieved the locally anisotropic Cosmological-Schwarzschild met- 
ric, both already discussed in the literature. It is shown that, for constant exponent a, the 
event horizon at the Schwarzschild radius r\ = 2GM/c 2 is only singularity free for a > 3 
and space-time is asymptotically flat for a > 5 which excludes these known cases. Also it 
is shown that, to strictly maintain the Schwarzschild mass pole at the origin r\ = without 
the presence of more severe singularities, hence describing a complete space-time with finite 
total mass-energy within a shell of finite radius, it is required a radial coordinate dependent 
exponent a.{r\) = «o + «i 2GM / '(c 2 r\) with a negative coefficient a.\ < such that at the 
event horizon, for o>q— \ai\ > 3 space-time is singularity free and for ao— |ai| > 5 space time is 
asymptotically flat. This metric may solve the long standing puzzle of describing local matter 
distributions in an expanding universe firstly addressed by McVittic. 

The curvature, curvature invariant and stress-energy tensor are analyzed in detail being derived 
the allowed bounds for the parameters Qo an d cxi that allow simultaneously space-time to be 
singularity free (except for the Schwarzschild mass pole at the origin) and the mass-energy 
density to be positive definite outside the event horizon. It is shown that, although space- 
time is locally anisotropic near the mass M, isotropy at spatial infinity is maintained. This 
characteristic is qualitatively consistent both with the experimental evidence of local anisotropy 
due to matter structures and global spatial isotropy. The modified Newton law for this metric 
is derived being shown that for planetary scales the usual General Relativity Newton law is 
approximately maintained for the full allowed range of the parameter ao while for galaxy 
scales and large values of the parameter a > c 2 / \f<lo{GM H ) 2 there is a significant deviation 
from this law which may contribute, for instance, for the flattening of galaxy velocity curves, 
hence allowing, at least partially, to describe dark matter effects, here interpreted as due to 
the universe expansion. Are derived solutions for planetary orbits both in the circular orbit 
and perturbative static elliptic orbits approximations and estimated the orbital precession, 
period corrections and time variation of the orbital radius which are well within the existing 
experimental bounds for the solar system. 
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1 Introduction 



In 1929 by analyzing the red-shift for radiation received from nearby galaxies Hubble [T] found 
experimental evidence that the universe was expanding at a rate given by the Hubble constant 
H = a/a, where o is the scale factor of the universe. Later, in 70's, it was realized that the 
expansion was not constant, at that time it was believed to be decelerating at a rate given 
by —qoHQ [2\, where qo is the deceleration parameter. More accurate measurements have been 
accomplished and today it is known that the universe is in accelerated expansion, hence qo < 0. 
These observations raised several theoretical (and phenomenological) puzzles both for small and 
large spatial scales. 

In particular, for small spatial scales, the definition of a metric that describes both an expanding 
background (the universe) and local matter distributions is today an open problem. As far as 
the author is aware there are only two distinct metrics that have been proposed. Not long after 
Hubble's work, in the 30's, McVittie derived a metric in isotropic coordinates [3] that preserves 
locally spatial isotropy interpolating between the Schwarzschild (SC) metric [4j which describes 
point-like massive bodies for small spatial scales and the Friedman-Lemaitre-Robertson- Walker 
(FLRW) metric [5] which describes the expanding universe for large spatial scales. It is known 
to be the only metric that maintains local spatial isotropy. Allowing for local spatial anisotropy 
(still maintaining spatial isotropy for large scales) has been recently proposed a Cosmological- 
Schwarzschild metric |6j. As has been discussed in [7i for the McVittie metric and will be analyze 
in detail in this work, both these metrics have singularities at the event horizon (the Schwarzschild 
radius). In this way the respective space-time is not complete and the metrics do not converge 
asymptotically to the SC metric near the point-like mass. For further discussions in this topic 
see also [U [9] [TOl [TTj [12] . Our main objective in the present work is to address this problem and 
try to solve it by deriving (or building) a new metric ansatz. We will take a phenomenological 
approach not having a fundamental theoretical basis for the results derived, the derivation of the 
new metric is only based in the assumptions that a metric describing matter in an expanding 
background (the universe) must interpolate between the SC and the FLRW metric maintaining 
space-time singularity free except for the SC mass pole at the center of mass. We will allow for 
local anisotropy as long as isotropy is recover for large spatial scales (specifically that the FLRW 
metric is recovered at spatial infinity). We note that this characteristic is not unwelcome, there is 
experimental evidence that local matter distribution does generate local anisotropy |13[ll4j. hence 
the imposition by McVittie for a fully isotropic metric seems excessive. 

With respect to large spatial scales, the same observations concerning the expansion of the universe 
led to today's cosmological inflationary models [15] . The deceleration is commonly attributed to 
the existence of dark energy which constitute 72.6% of the total mass-energy of the universe. The 
most common description of dark energy is through the cosmological constant [16] which describes 
a constant background (vacuum) energy with negative equation of state. Alternative approaches 
do exist that conciliate inflationary theories with the existence of a non-null vacuum energy. 
Such examples are quintenssence [IT], phantom energy [18], modified theories of gravity [19J and 
string based/inspired theories [20]. For a review in these topics see [21]. Although not originally 
related to expansion there is also for large and/or very massive astrophysical systems experimental 
evidence of non-radiating matter [22] moving at non-relativistic speed, usually named cold dark 
matter (CDM), which is responsible for significant deviations with respect to the usual General 
Relativity predictions, namely causing the flattening of galaxy rotation curves [23], increase of the 
gravitational weak lensing effects by astrophysical large structures [24] (as is the case of the bullet 
cluster |25] ) and the modulation of the weak lensing effects near heavy astrophysical bodies [26J. 
The most direct interpretation of these effects is the existence of physical non-radiating matter 



4 



that as not been accounted for, hence named dark matter. Alternatives to this interpretation 
are based in modifications to the usual General Relativity Newton law, specifically there are two 
distinct approaches, either known as MOND [27] for which the General Relativity Newton law 
becomes proportional to the inverse of the radial distance for galactic scales and the Scalar-Tensor- 
Vector theory [28] which considers extra fields interacting gravitationally. For reviews in these 
topics see |29| . It is interesting that the metric we derive in this work, although based in very 
simple assumptions, also allows at galactic scales for a significant deviation from the usual General 
Relativity Newton law, hence in our framework these effects may be interpreted as a consequence 
of the universe expansion. 

It is today well accepted that both dark matter and dark energy contribute to the universe evolu- 
tion. In particular this relation led to the most successful model for today's universe, the KCDM 
model [30] which accounts for the contributions of these forms of matter and energy to the total 
universe matter-energy density and pressure. Based in this model and the combination of distant 
supernova of type la red-shifts (SN) measurements [3l], the baryon acoustic oscillations (BAO) 
measurements [32] and the five-year Wilkinson microwave anisotropy probe (WMAP) measure- 
ments |14] the parameters for the expansion of today's universe can be estimated. According to 
these estimative the universe is flat and we are taking the following experimental values for the 
Hubble rate 



H= a - , H =« 
a a 



= h Mpc' 1 Km s" 1 = 2.28 ± 0.04 x 10~ 18 s' 1 , (1.1) 

t=to 



and the deceleration parameter 



a I a\ a 1 



, 10 — TTo 



(1Q 



t=to a H 



-0.589 . (1.2) 



This last value is computed from the equation of state for the universe and we will address the 
details of this calculation later on. Similarly the third derivative of the scale factor can be defined 
in terms of a parameter s 



a I a\ I a\ a I 

so 



a \a J \a I aqn 6 



t=t 



Q-o 1 / 
•— r , (1-3) 



which encodes the variation of the deceleration of the universe. As far as the author is aware there 
are no experimental estimative for the current value of so- We also refer that the characteristic 
lengths of the expansion effects are the Hubble time 



tH= H ' tHo = ~k 

and the Hubble length 

c , c 



1 

— 1 

t=t 



H 



4.39 x 10 i7 s , (1.4) 



H ' J1U H 



c 



Ih = TF > Ih = T7 = -FT ~ L32 x 10 m ■ ( L5 ) 



t=t Hq 



Depending on the model being considered to describe inflation these correspond approximately (or 
at least give an order of magnitude) to the age of the universe and the distance to the cosmological 
horizon. 

The information from beyond a distance of the value of the Hubble length cannot reach the 
observer within the life-time of the universe, hence it is causally disconnect from the observer. In 
cosmology this is known as the horizon problem which, together with the flatness problem (why the 
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universe is flat today?) led to theories where the fundamental constants vary with time [33j GS1 ES] • 
In particular theories of varying Newton constant G may be relevant for planetary system physics. 
For General Relativity in flat backgrounds, the variation of this constant is equivalent to a time 
varying orbital radius, hence experimental measurements of either of this variations are equivalent. 
When considering planetary orbits in an expanding background this is no longer the case and the 
relation between a varying G and a varying orbital radius is distinct. In particular, as we will 
show, for an expanding background the orbital radius varies with time for a fixed G. We note 
that this discussion is relevant when interpreting experimental data or setting bounds on these 
variations. 

We have organized this work as follows: 

In section [2] we review the description of an expanding universe in terms of the FLRW metric. 
In particular we discuss in detail the relation between the choice of coordinates and the physical 
observables defining the several coordinate choices used in the remaining of the work. We also 
discuss other relevant issues as the derivation of the stress-energy tensor for the FLRW metric 
and the derivation of the deceleration parameter from the cosmological parameters. 

In section [3] we review how matter is usually described in terms of a metric both for flat and ex- 
panding backgrounds. We start by analyzing the Schwarzschild metric deriving the most relevant 
quantities for this work. In particular we review how to derive the General Relativistic corrections 
to the precession and period of planetary orbits. We also introduce the McVittie metric and the 
Cosmological-Schwarzschild metric for expanding backgrounds discussing their properties. 

In section |4] we build a locally anisotropic metric ansatz interpolating between the SC and FLRW 
metric. We start by a one parameter metric and analyze the properties of the respective space- 
time. In order to regularize the essential singularities at the origin it will be required to refine the 
ansatz by introducing a second parameter. Once this is accomplished we derive the stress-energy 
tensor explicitly analyzing its properties, in particular we discuss the local anisotropy and define 
the allowed relative bounds of the parameters that ensure mass-energy positiveness outside the 
event horizon. 

In section[5]we derive and analyze the modified Newton law for this new metric. First we study the 
radial acceleration component due to a central mass and how the transition between the General 
Relativity Newton law and the Newton law for the expanding background (obtained from the 
FLRW metric) is affected by the metric parameter values. We further analyze the orbital motion 
for the new metric both for circular orbits and static elliptical orbits approximation estimating 
the time variation of the orbital radius and the corrections to the orbital precession and orbital 
period comparing these with available experimental data. 

In section [6] we resume the results obtained in this work and discuss further directions of research 
in this topic, namely the possible relevance of our results for dark matter effects. 

We have gathered the technical details of the main calculations in the appendixes. In appendix IA1 
we list the conventions and definitions considered in this work. In appendix [B] we present the 
connections, curvature and Einstein tensor for the FLRW metric and derive the lower order series 
expansion for the time evolution of the Hubble rate. In appendix[C]we present the McVittie metric 
for several coordinate choices. In appendix[D]are listed the connections, curvature, Einstein tensor 
and curvature invariant for the new metric, both for the one parameter and two parameter ansatze. 
Are also derived the differential equations for stationary orbital motion. 
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2 Reviewing the Background Cosmological Metric 



In this section we review the theoretical effects of space expansion for free travelling massive 
particles and radiation interacting only with the cosmological background. These effects are well 
known and can easily be found in the literature [3SI 133 EE]. We try to approach this subject in 
a pedagogical fashion from a non-expert perspective. We discuss distinct coordinate systems by 
rescaling the spatial coordinates by the universe scale factor a, in particular we consider expanding 
coordinates (for which spatial coordinates expand over time) and non-expanding coordinates (for 
which the spatial coordinates do not expand over time) as defined in the appendix [AJ This dis- 
cussion is relevant whenever we want to compare theoretical results with physical measurements. 
Also, for technical simplification of the calculations, different coordinate choices are often consid- 
ered independently of the physical measurable and respective interpretation. For these reasons it 
is relevant to discuss what coordinate and unit choice is more adequate when studying a particular 
physical system and relating theoretical and experimental results. 

It is known that the physical laws do not depend on frame nor coordinate choice in the sense that 
the relations between different quantities remain the same, however the way we measure those 
quantities depend on the units and coordinates and consequently our interpretation of physical 
phenomena, as well as the definition and values of the fundamental constants, does depend on 
our measurement procedure [35J. As simple examples, considering to use natural units c = H = 1 
is not suitable when trying to measure in the laboratory the speed of light and radiative decays 
with a clock that measures time in seconds. Also if we have a rod that measures lengths f = 1/r 
(being r the usual distance for which the area of the sphere is At: r 2 while for f it is An/f 2 ) the 
classical Newton law of gravitation would be stated as that the gravitational acceleration has 
direction along growing f being proportional to the cubic distance plus the velocity over distance. 
Specifically we would obtain 

GM r 
r = 5- f = +GM f 3 + - . (2.1) 

Hence the physical interpretation does depend on the coordinate and unit choices which reflect 
the way we interact with reality and how we perform our measurements of the physical quantities. 
Nevertheless the laws of physics are invariant under a change of coordinates, in the sense that 
once we obtain a law, it is still valid under any transformation of coordinates applied directly to 
the specific physical relation, as we have just exemplified in the above equation for the Newton 
law. If necessary, for technical simplification or any other reason, we can use any coordinate 
system we wish for calculation purposes, as long as we recall to transform the final quantities (or 
expressions) back to the coordinate system that is directly related to our reality. Of particular 
relevance to the present work is the definition of measurement for space and time lengths which 
dictates how we perceive our universe and the coordinate choice(s) which translate between our 
theoretical constructions and reality (or the other way around). 

In General Relativity the specific expressions for the metric components depend on the coordinate 
and unit choice. In particular when we are living in an expanding space it is relevant to know if 
our measurement rods are also expanding or not. A material rod is surely expanding with space, 
hence we will not be able to measure any expansion with it. However radiation is sensitive to 
space expansion, it was through the red-shift of radiation emitted by distant objects that the 
Hubble rate (II. ip was originally measured pQ. 
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2.1 Working with Non-Expanding Coordinates 



We will take the global geometry of space-time to be null n = which, based in experimental data, 
is today the most probable case [13] . The usual metric taken to describe a flat homogenous and 
isotropic universe for an observer at rest at the origin of the coordinate frame is the Friedmann- 
Lemaftre-Robertson- Walker type metric 

ds 2 = c 2 dt 2 - a 2 dr 2 - a 2 r 2 (d9 2 + sin 2 0d<p 2 ) = c 2 dt 2 - a 2 5 ij dx i dx j , (2.2) 

We will henceforth refer to it as the FLRW metric (note that, for k = 0, this metric is often 
referred only as Robertson- Walker metric [5]), here expressed both in spherical coordinates and 
Cartesian coordinates. For this coordinate choice we have the geometrical identification with a 
time-dependent area of the sphere A(t) = 4tt a 2 r 2 (with commoving radius ar). We refer to this 
coordinate choice as non-expanding coordinates, although the three-dimensional space lengths do 
expand, as will be shown next, the coordinates (as given by the coordinate equations of motion) for 
massive particles do not. Here a = a(t) stands for the time dependent scale factor of the universe 
as appearing in the definition of the Hubble rate (jl.ip . a is dimensionless and we assume that at 
some reference time to corresponding to the most recent estimative for the size of the universe its 
value is normalize to a(to) = clo = 1- This is achieved either by rescaling the spatial coordinates 
by this constant or, alternatively, to consider a ratio a/ao (instead of a) in the original metric. 
However, unless otherwise stated, we will keep it in the equations to keep track of the several 
technical steps and identify how the scale factor a = a(t) affects the results discussed in this work. 
Also, as much as possible, we will derive model independent results such that a is not a metric 
degree of freedom being some generic time-dependent function describing the background. The 
technical details for this metric including connections and the Einstein equations can be found in 
appendix [Bj 

Let us start by defining the three-dimensional geometrical lengths corresponding to this metric. 
This length is defined through the internal product corresponding to the three-dimensional (spa- 
tial) hyper-surface for some fixed time t. Specifically the three-dimensional intrinsic metric for 
the FLRW metric (|2.2[) is ^'Qij = — g%j = a 2 <% and the coordinates in this spatial hyper-surface 
coincide with the four-dimensional spatial coordinates, ^x l = x % . We use the suffix '( 3 )' to dis- 
tinguish between the four-dimensional and three-dimensional metrics and coordinates. Then, in 
the three-dimensional hyper-surface, the square of the geometric length / ge0 m is 

^jj'^geom^geom — ^ ^ 9ij ^ ^ ^ ^ — Q 5{jX X^ . (^-3) 

As for time lengths, neglecting the relativity effects (7 = dt/dr ~ 1), the proper time and 
the coordinate time coincide. Then, given this approximation, any measurement of coordinate 
time coincides with measurements of proper time. Hence we obtain the following definitions for 
geometrical lengths, geometrical velocities and geometrical accelerations 

^geom OX, 

v geom = CLX , (2-4) 



The geometrical length is the physical observable length for the FLRW metric (|2.2p . We do not 
measure coordinates lengths, once a metric is considered we measure these geometrical lengths. 
However, as we will analyze in detail, neither the geometrical velocity, neither the geometrical 
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acceleration coincide with the physical measurable velocities and accelerations. Instead these will 
be given, as expect in the framework of General Relativity, by the geometrical covariant velocities 
and acceleration. 

To derive the classical limit, for which an Euclidean metric 5ij can be used, the definition of three- 
dimensional coordinates (in the spatial hyper-surface for each fixed t) is not enough. We need 
to consider the above definitions of geometrical spatial lengths that play the role of the classical 
local Euclidean coordinates in the Newtonian limit. Here local means that this construction is 
only valid in the neighborhood of the space-time event we are considering. Specifically we deal 
only with time-independent metric components, hence our results will be valid in a neighborhood 
of some time to- More generally when the spatial curvature is significant we can, at most, define a 
Euclidean metric in the neighborhood of some spatial point x l (the most common example is the 
earth surface, being a sphere, only locally it is approximately flat). In the following we use the 
notation l l for the exact expressions and l l for the respective classical quantity (an approximation 
to ^mass) that represents the Newtonian limit. 

To further proceed without specifying a cosmological model that would fix the expression for the 
scale factor a, let us consider a series expansion of a p in the neighborhood of the reference time to 

a* = <fl+A*-*o) + ?f^ + (p-l)f^yV*-*o) 2 V (* 3 ) 

V «o 2 ya \a J J J ( 2 .5) 

« ag(l+pflo(t-to)-f(go-P+l) Hl{t-t Q f) +O{(H (t-t )) 3 ) . 

Here ao = a(to), do = d(to) and do = d(to) are the scale factor and its time derivatives evaluated 
at the time t = to- In the last approximation we have replaced these derivatives by the Hubble 
rate Hq = do/ao and deceleration parameter q$ = — do/(ao Hq) also evaluated at t = to- Taking a 
phenomenological approach we assume that the values for these parameters are known at time to 
and it is not necessary to actually specify a cosmological evolutionary model to study and estimate 
the effects of the expansion today. This is a valid approach as long as we are dealing with short 
time scales t < 1/H ~ 10 17 s ~ 10 10 year and spatial scales x l <C Ih = c /Ho ~ 10 m. 

We can now compute the time evolution for a coordinate point x with null coordinate velocity 
x l . The evolution of the respective geometrical spatial lengths l in the neighborhood of t = to is 
given by 

l (t) = a(t)x « l (t) = J* (0) + l* m H (t - t ) - \qol l 0(0) Hi(t - t ) 2 + 0(t 3 ) , (2.6) 

where = clqx is the spatial length at t = to, we have considered the series expansions (j2.5[) for 
p = 1 and explicitly wrote the time dependence of the several quantities. The underscore indexes 
'0(0)' represent respectively the fixed coordinate index (x ) and the evaluation time t = to- Then, 
as expected, the length will increase with a velocity given by the product of Hubble rate Ho by 
the initial length l l ^ and an acceleration given by the product of the deceleration rate —qoH^ by 

the same initial length /q(o)(*°)- 

To compute the geodesic path for massive particles let us first compute and solve the equations 
of motion for a particle travelling with some given initial velocity x l as observed by a static 
observer at the origin of the coordinate frame, ^ bs = (0,0,0). We assume that the only physical 
interactions are due the expanding background (given by the metric) and static is meant with 
respect to that background such that the observers coordinate velocity is x* bs = (0,0,0). In 
this section we use Cartesian coordinates. The connections are given in equation (|B.1[) and the 
generic expression for the relativistic factor 7 = dt/dr is given by equation ()A.10p . which for the 
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metric (12.21) is 



dt / a 5j4X l x J 



^=f T = {i-~-^) ■ (") 

The equations of motion for the spatial coordinates are given by the usual geodesic equation as 
expressed in (|A.lip . For the FLRW metric (|2.2p we obtain 



& = -2cT i Qj ± j - j~ 1 jx i » -2- » -2x* H . (2i 

a 



Here we have approximated the time derivatives of the scale factor by the Hubble rate evaluated 
at the reference time t = to, Hq = do/ao (jl.lp . and we have taken the non-relativistic limit 
considering 7 ~ 1 and 7 ~ 0. This limit corresponds to non-relativistic speeds i 1 < c and short- 
scale spatial distances such that the spatial coordinates are well bellow the characteristic length 
scale of the system x l <C Ih (jl-5j) . An estimative for the relativistic corrections is 

1 i fda5jhX^x k a 2 5jhX^x k \ 9 ,■ 

X \ el = - 7 -v* = + (^! — +—^ 2 — wv 

(2.9) 

a 5 jk x j x k 2 2 • / a 2 x 2 \ a 2 x 2 rr ., / x 6 H 



c 2 \ c 2 / c 2 V c 6 



where in the last line we have replaced the second time derivative x % by the equations of mo- 
tion (|2.8p . expanded the relativistic factor in first order in x 2 /c 2 and used the short hand notation 
K For non-relativistic coordinate velocities (x 2 <C c 2 ) and small-scale spatial coordi- 



nates {x l <C Zff ) this correction is negligible with respect to the above expression (|2.8 
Integrating the equations of motion (|2.8p we obtain 



- = -2- **= - Mo- (2-10) 

x l a V o / 

Here we have set the integration constant to a^o = a (*o)^ J (*o) which corresponds to the equation 
evaluated at the reference time t = to- Given the expansion (|2.5p with p = — 2 we can further 
integrate this equation in from the reference time to to some generic time t obtaining the coordinate 
path for massive particles 



dx' 

XQ J to "~ J t 



ft d+' ft 

/ — « i / dt' [1 - 2F (t' - to) + (90 + 3) F 2 (t' - t ) 2 

Jtn a Jtn 



*Lss(t) ~ 4 + 4((*-*o)-#o(*-*o) 2 + ^(<7o + 3) ^o 2 (i-*o) 3 ) + O(ff 3 t 4 ) • 



(2.11) 



We use the notation x' and t! to distinguish between the integration variable and integration 
limits x and t. We note that the coordinate acceleration (|2.8p cannot be directly interpreted as 
the Newtonian acceleration and, as we have already discussed, similarly to the coordinates and 
coordinate velocities, cannot be directly measured, we measure the geometrical acceleration ()2.4p . 
We will return to this discussion later. 

Finally we can define the geometrical lengths Z^ass f° r a massive particle describing its trajectory. 
In order to do so we take the solutions for the coordinate equations of motion (|2.1ip and the series 
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expansion for the scale factor (|2.5p obtaining 



mass 



*L« = ^o(o) + (* " *o) + \a c (t - to) 2 + O (t 3 ) 



/_ 0(0) = a x = Z geom(*o) , (2.12) 
= ao (x l + Xq flb) = ^gcom(^o) + ^geom(^o) #0, 
«c = Ceom(*o) -^0 • 



We interpret = Z l m ass as the classical velocity and a l c = l l mass as the classical acceleration in 
the neighborhood of to- Only the spatial coordinate lengths at the origin coincide with the 
geometrical spatial lengths evaluated at the initial time to- Neither the v l c nor a* coincide with 
the respective geometrical quantities, instead are given by the sum of the background velocity and 
acceleration (|2.6j) with the particle geometrical quantities (|2.4p . For the example just presented 
only the geometrical velocity Ugeom(^o) = °o^o * s P resen t, more generally we can consider as well 
some external acceleration such that a\ = a ge om(*o) — Qo ^o(o) Hq with a geom (to) = oo^ext(^o)- 

This result was actually expected, we recall that in General Relativity the physical observable 
velocity and acceleration are given (at all times) by the projection of the respective covariant 
quantities. Hence let us show how these quantities are related with the classical quantities (|2,12[) 
computed above. The covariant velocity and acceleration are computed as usual by considering 
the parallel transport derivatives 



Dx l Dx 



Dt ' Dt 



=l{x i + cV 



= 7 2 ( x l + 2-x i + -x l + 7 _1 7 ( x i + -x i 
\ a a \ a 

i (a a _i . a -\ a 
= 7-2+7 7-xM Ri -x % . 

\a a J a 

In the final expressions we used the equations of motion for x l (|2.8[) and x % (|2.10p and in the final 
approximations we have neglect the relativistic corrections setting 7 ~ 1 and 7 ~ 0. Considering 
the projection for these vectors into the spatial hyper-surface at fixed time to, we obtain the 
expressions for the respective classical quantities valid in a neighborhood of to 

( 3 ) ffii (3) v i(3) v i = a 2 S . jV i v j ^ „i = now* fa) , 

(2.14) 

(3)0y (3V ( 3 V = a 2 SijSfa? => a\ = a a*(t ). 

In this way we confirm that at the reference time t = to the classical spatial Euclidean coordinates 
(the spatial lengths) l l , velocities v\ and accelerations a\ do correspond to the spatial coordinates 
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x 1 , covariant velocity v* and covariant acceleration a 1 projected to the spatial hyper-surface at 
t = to- Taking in consideration the corrections due to the expanding background, the usual 
laws of classical mechanics are valid in a neighborhood of to as expressed by (|2.12p . However it is 
important to stress that these values are only valid as an approximation. More generally neither of 
these quantities are constant over large time scales (t ~ I/Hq) neither large spatial scales [x ~ ljj ) 
for which it is necessary to properly take in consideration the dynamics of cosmological evolution 
of the universe [15, 21J. From the equations of motion (|2.8p we obtain the expression for the 
coordinate velocity to be x l = x l /a 2 and, as we have just discussed, due to the internal product 
being given in terms of the intrinsic metric, we obtain the usual expression for the coordinate 
peculiar velocities t> geom = ax 1 = ao (aox )/a = ao v georn (tQ) l /a observed experimentally between 
astrophysical objects [35] • These are due to the expansion of space being exact within General 
Relativity and have no classical analogy. In addition, either assuming that we can estimate 
and subtract the coordinate velocity ax 1 , or in the limit of relatively large distances and small 
coordinate velocities x 1 <C l l geom or simply by having a large enough statistical sample for which 
the coordinate velocities average to zero, from the expression for the covariant velocity (12.131) . we 
retrieved the Hubble law for lengths (or distances). The relation between the length Z geom = ax 1 
and its time variation Sl l geom /St is 



geom i a ri /<-> -i tr\ 

— = a— =av ^-/ geom . (2.15) 

The approximation in the last equality is only due to have neglecting the peculiar velocities and 
the variation 5, as stated in the first equality, considers both the parallel transport of spatial 
coordinate and the projection to three dimensions. We stress that this relation is exact within 
the framework of General Relativity being valid at all times, only for short time scales and spatial 
scales the rate of the scale factor can be approximated by the constant Hubble rate Hq = do/ a o- 
Taken the classical limit, as given in (|2,12[) . this relation is (exactly) valid only at the time t = to, 
for which l l {t = to) = v l c ~ H{t = to)^o- We remark that this relation is coordinate dependent in 
the sense that, for a different coordinate choice, the respective spatial lengths will be expanding 
at different rates. We will discuss this issue later on. 

So far we have not discussed any measurement procedure. Let us discuss two common methods for 
measuring distances. Specifically we address two kind of light measurements, range measurements 
where the travelling time is measured (usually the proper time using atomic clocks, see for instance 
section 16.4 of [37] for further details) and the distance computed by considering the speed of light, 
and red-shift measurements which allow to directly measure the variation of the wave length due 
to the expansion of the light path from the emission point to the reception point. The later was 
the original technique used when setting the Hubble law. 

Let us start by defining the line element for light travelling in the background metric (12. 2\ along 
some coordinate direction x. A light-like trajectory corresponds to the shortest possible path 
ds 2 = 0, then from (|2.2I) we obtain the infinitesimal relation and respective differential equation 

cdt = adi^dx 1 ^^- = -e\ (2.16) 
dt a 

where e l stands for a unity vector with respect to the Euclidean metric 5ije l e :) = 1 and is introduced 
to define the propagation direction of radiation. This equation can be reduce to one single direction 
x, for which the internal product x.e reduces to an overal ± sign 5ije l x J = ±|x|. Again, to keep 
the derivation model independent, let us consider the expansion for the scale factor (12. 5p with 
p = — 1 valid in the neighborhood of to and integrate the expression from the reference time to to 
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some genric time t obtaining the light coordinate at time t 

4 gM (*) « 4 + — T fl " #o (*' - *o) + J (90 - 2) F 2 (f - t ) 2 + O ((H$(t - t ) 3 ; 
ao Jt L 2 

= 4 + — ((t - t ) - \Ho (t - t ) 2 + Uq - 2) H 2 (t - to) 3 + O {H*(t - to) 4 ] 

(2.17) 

We recall that we are considering an observer at the origin of the coordinate frame and our aim is 
to measure the geometrical distance between the point x l (at t = to) and the observer for which 
the coordinate is null (at some time t). Hence the coordinate distance for light is 

Alight = x — ^Hght (*) • (2.18) 

From this coordinate length we can infer the respective light length corresponding to the point xq 
from where the light signal was emitted 

lf ight = aAx\ ight = - (ce ! ) (t - t ) - \ (c^o) (t - to) 2 + 0(t 3 ) . (2.19) 

We note that due to light being travelling to the origin of the coordinate frame, e* has the opposite 
sign of the value of the coordinate x l which is correctly offset by the overall minus sign. 

The above light length does not correspond to the geometrical spatial lengths (|2.6p . this is eas- 
ily explained by noting that the spatial expansion affects the geometrical lengths by a factor 
+xHq (|2.6I) proportional to the coordinate x, hence for light travelling from to x l = 

we obtain an average effect such that Ho (x + 0)/2. Also we note that the geometrical lengths 
are time dependent such that they do not match at the time of emission and time of reception 
of the light signal, respectively are l (t = to)^o(o) = a o x b anc ^ ^o(*) = a (t) x h- The above equation 
estimates the lengths at the time the light is received (t in the previous equation) since it is already 
corrected (projected) by the scale factor a. 

Therefore when computing distances by considering the light travel time At there will be a cor- 
rection to the usual classical light length Z (o) = a o x o = c ^t both for the geometrical length 
corresponding to the light signal emission (|2,19p as well as for the estimative of the geometrical 
length distances to the massive particle (|2.12|) . Specifically, to first order in Ho, we obtain 

Vo) = , 



flight - V^light^lighi 



i„ 



l m + -H ^f- + O(tf 2 /J 0) / C ) , (2.20) 



r 

= +l m+ Ho^ + o(H 2 l 2 Q) /c) . 

Due to the very small value of the ratio Hq/c = 1/Ih ~ 2-7 x 10~ 27 m _1 , these corrections are for 
most purposes negligible being usually below experimental accuracy. In figure [1] are represented 
the evolution of the several spatial lengths during a light range measurement. 
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Figure 1: Evolution of spatial geometrical lengths I in a range measurement corresponding to the spatial 
coordinates Xq (fixed at time to ), light coordinates x\ ight and massive particle coordinates x l maBS from the 
emission time t e to the reception time t r as given by equation i2.20\) . 



For the case for which the signal is emitted by the observer, reflected by the massive particle and 
received back at the origin, the problem can be symmetrized with respect to the time at which 
the reflection takes place, both the geometrical lengths given in equations (|2.12p and (|2.19|) are 
symmetric with respect to time inversion. Then, taking the reflection time to be to, we have that 
the emission time t e and reception time t r are spaced apart from to by the same time length, 
to — t e = t r — to- Hence, for fixed Hubble rate Hq = H(to), the measured distance is the same in 
both directions and we get twice the travelling time corresponding to twice the light length Zg. 

We proceed to briefly compute the cosmological red-shift. The most straight forward way to do 
so is to consider the 4- momentum for travelling radiation with frequency cv/2ir and wave number 
k l , k^ = (u>/c,k l ). Its square is null g^k^k" = lo 2 /c 2 — a 2 k 2 = (light is massless), such that by 
energy-momentum conservation we obtain the following optical dispersion relation and respective 
solution for the wave number 

u = ca(t)yj %kW => k l = ^ . (2.21) 

The direction of propagation for the light is given by the wave vector k which corresponds to the 
evaluation of k l at some reference time to- Alternatively, with out specifying a reference time it 
can be set to k = u e l /c where e l is a unit vector {5ije l e 3 = 1) which defines the direction of light 
propagation. By further noting that the wave vector is inversely proportional to the wave length 
A k % = 2-7T e l /X, we obtain the expression 

X(t) = ^ = ^A , (2.22) 
c a 

where Ao and ao stand for the wave length and scale factor values at some time to. Then the usual 
cosmological red-shift z is straight forwardly defined as 

X r — A e ~a~^ e ^ e Or Ap — O e A e Ou T * ,„ „„ s 

z = — — - = ^— = _LJ_J_J_ = JL - i ) (2.23) 

A e A e ct e A e a e 

where the indexes e and r stand for emission and reception of the radiation. Considering the series 
expansion of the right-hand side of this equation we obtain, to first order, the red-shift Hubble 
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law 

z « H At - igo^o A * 2 + 0(H$At 3 ) 

^ (2-24) 

C £ c 

In the last equality we have replaced the time interval At = t r — t e by the classical length interval 
Zq(o) = cAt (|2.20p . The corrections to the classical lengths are mostly relevant only for observation 
of very far objects and, otherwise, are usually below the experimental accuracy. 

Although these deductions are helpful, both as a reference and to introduce several techniques, 
we note that the historical path has been the reversed. The Hubble law was set experimentally [I] 
implying that the universe was expanding which led to the derivation of the FLRW metric (|2.2|) [5] . 



2.2 Working with Expanding Coordinates 

We have just shown that for non-expanding coordinates x corresponding to metric (|2.2p . the three- 
dimensional coordinates x l do not expand, for a fixed coordinate point x l = the coordinate 
velocity (|2.10p and acceleration (|2.9j) are null. The physical measurable quantities correspond to 
the respective geometrical quantities and do expand with the background, the geometrical velocity 
is corrected by the Hubble velocity due to expansion, +H x 1 , and the geometrical acceleration by 
the acceleration of the background, —qH 2 x l (j2.14j) . Classically it is only possible to define local 
corrections that match these background effects at some fixed time (|2,12p . These results raise 
the question whether new four-dimensional coordinates can be defined that coincide with these 
geometrical lengths which would have a more intuitive physical meaning. We also expect this new 
coordinate choice to simplify the technical details and interpretation of the several results. 

The new set of three-dimensional spatial coordinates x\ that we are considering in this section is 

x l = ^ , r = ^ . (2.25) 
a a 

We name these expanding coordinates due to their value being changing over time with the scale 
factor a. Here the underscore index '1' stands for a new coordinate set following the convention 
describe in appendix For spherical coordinates this coordinate transformation is equivalent to 
transformation (|A.13p . We note that the transformation (|2.25p corresponds to the (simple) three- 
dimensional coordinate transformation. The respective generalized coordinate transformation 
h^v for four-vectors and other tensorial contravariant quantities is computed from the linear 
transformation for the infinitesimal forms dx° and dx l 

dx a = dx\ , tf = ^-^^, (2.26) 
a a a c 



such that we obtain 



A°o = S£ = l , A" = S = 



dx\ ' 1 dx\ 



a i ® xl _ d x\ 1 ^ _ dx l , d xl 1 

" /~)n"Q ft ft C 3 £i™7 ^ CL CL C ^ 



(2.27) 



dx\ a a c 3 dx{ 



where, in the last component, no summation over the repeated indexes j is implied. As usual the 
transformation A for covariant tensors is obtained by inverting the transformation matrix A. 
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In particular, the metric for these new coordinates, is obtained either from the infinitesimal forms 
transformation ([2.26|) applied to the infinitesimal length square ds 2 or the transformation A applied 
to the covariant metric components g^ u . Specifically it is 





s 5ijx\x{ j dt 2 + 2c ^— $ijx\^j dx\ dt — 5ijdx\dx{ . 



ac 



Where we are omitting the index '1' on the x° coordinates, the coordinate transformation is trivial 
for the time components (x° = x®). 

When dealing with this metric, the main technical difference with respect to the previous sub- 
section, is due to the non-null metric component g^. In particular when considering a projec- 
tion to the three-dimensional hyper-surface it is necessary to decompose the four-dimensional 
infinitesimal square length into a time length minus a proper spatial length, specifically ds 2 = 
c 2 dt 2 — (dxi — d/adt) 2 . It is this spatial length that corresponds to the three dimensional mea- 
surable length 



geom.l = - \9ijx\x\ - gojx\x{ ) = 5ijx\x{ ( 1 - -t ) . (2.29) 



a 



a 



2 



However this result is not enough to proceed. It is desirable to have a projection definition that can 
be applied to generic vectors, including velocity and acceleration as was considered in the previous 
section. There are several ways to define this projection (see for example section 21.4 of [37]). We 
use the construction originally suggested in [39J such that the definitions for the intrinsic metric 
^9ij an d the three-dimensional spatial projection for any vector is given in equation (|A.12p of 
the appendix El For the specific case of the spatial coordinates we obtain 



(3) <^ = <% , (3) x\=xi[l--t) . (2.30) 



a 
a 

We also note that the off diagonal term g^ usually represents a frame velocity ff rame and, when 
this velocity is constant, a diagonal metric can be considered by shifting the space coordinates 
x\ = x\ — v l ilSLme t. In this case this velocity is both time and space dependent, nevertheless we can 
expand this term recursively and, by successive shifts of the spatial coordinates x\, to obtain the 
same definition for the three-dimensional coordinates to any desired order (in powers of Hq). 

For this coordinate choice, neglecting Special Relativity effects, the proper time coincides with 
the coordinate time, r = t. Then, given the above projection (|2.30p . the geometrical lengths, 
geometrical velocities and geometrical accelerations are defined as 

/,' , - ( .1 - -t 



geom.l 



a 



«Ui = *i - ~ iS = *i " ' (2 ' 31) 



j ..j U X^ „Q 

a geom.l — x l ~ 



In the last equalities of the last line we explicitly expressed the derivatives of the covariant coordi- 
nate Xi = ct by the respective time coordinate expressions X® = 1 and = 0. As before, in order 



16 



to obtain the classical limit, it is necessary to expand the projection factor 1 — a/at. However we 
note that in a neighborhood of a generic time to, it will contain factors depending linearly in this 
value, specifically we obtain 

l- a -t = (l-^t ) + (-^-% 0+ (^) 2 t )(t-t ) 
a V a J \ a ao V a o/ / 

(2.32) 

+ (-« + (*) 2 -i%+™^-(*)\) (t-t f+o(t 3 ). 

\ ao \ a o/ 2 2 a a 2 \a J J 

The dependence on to is not desirable since this value is, in principle, arbitrary. Depending on 
the physical system we are considering it can be either the time we start a specific experiment or 
one century ago. There are several ways to deal with this issue, the simpler is to shift the time 
coordinate by to and perform the expansion in the neighborhood of t = 0. We note that none of 
the metric components change under this coordinate transformation due to the time dependence 
being encoded in the dimensionless scale factor a. Then the expansion in the neighborhood of 
t = is 

1 - ^ t » 1 - H t + (q + l)H^t 2 + O((H t) 3 ) . (2.33) 
Then the geometrical length l\ for a spatial coordinate point x\ (fixed to the background) is 

i\.o(t) = 4.q(i-^) «iU*) 

= 4.0(0) " 4.0(0)^0 1 + (q + 1) x\ m H 2 t 2 + O((H tf) (2-34) 

= 4.0(0) + «i.o(o)ffo t + (go + 1) l^H 2 1 2 + O((H tf) . 

The two underscore indexes 'i.o' in the coordinate stands for the coordinate label (xi in this 
case) and the evaluation time (to), xi.o = xi(t = to), while in the geometrical lengths we use the 
notation 'i.o(o)' to specify the coordinate choice, the type of length (as in the previous section '0' 
for a coordinate fixed to the background, 'mass' for a massive particle and 'light' for radiation) and 
the evaluation time (inside the brackets). In the last expression we have replaced the coordinate 
quantities by the respective geometrical quantities as defined in (|2.31j) with x\ ^ = 0, the 

length at the initial time coincides with the geometrical length evaluated at t = 0, Z| ^ = 
Zg eom l (i = 0) = x\ and the velocity corresponds to the geometrical velocity also evaluated at 
t = 0, v\ 

0(0) — u geom.i(* — 0) — —x\ Hq. As for the geometrical acceleration it corresponds to 
the spatial expansion acceleration for this coordinate choice. 

The remaining of the calculations closely follow the ones of the previous section. However we repeat 
them here. The main reason to do so is to clearly identify the relation between the several results 
derived in both expanding coordinates x and non-expanding coordinates x\ allowing to identify 
its relation with our world physical measurable as well as to identify the technical advantages of 
each choice. 
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The relativistic factor 7 and its derivative with respect to the time coordinate 7 are 
dt 



1= Tr 



1 



d\ 5jjx\x{ + 2 fa\ Sijx\x\ J 5ijx\x\ J 
a I c 2 \a 



7 




1 /. a 

1 + 2? ( Xl -a Xl 




a 

Xx Xx 

a 



(2.35) 



The generic spatial coordinate equations of motion are given by the geodesic equation (jA.llj) of 
appendix[A]and the connections for this metric in Cartesian coordinates are given in equation (|B.4p 
appendix [Bj Hence we obtain the following equations of motion 



-c 2 r 



00 



2cP 



or 



T\ jk ±ix k 



— 1 ■ • i 

7 7X 



+ -x\ 
a 



a 

+ -x 

a 



-go -Ho 4 



Xi 



-Xx 



— 1 ■ ■ i 

7 7X 



(2.36) 



In the last equality of the last line we replaced the second derivative of the scale factor by the 
deceleration parameter value at the reference time to, go (|l-2p and from the second to the third 
line we took the approximation for non-relativistic limit, coordinate velocities much smaller than 
the speed of light, x l <C c, and spatial distances much lower than the characteristic length (Hubble 
length) x l <C c/Hq = Ih (jl.5j) such that 7 w 1, 7 w and (x — Hqx) 2 /c 2 <C go. The relativistic 
corrections corresponding to the second and third terms of the second line can be evaluated by 
replacing the expression for 'x in the time derivative for 7 (|2.35p 



'l.rel 



Id/. a 

( Xx Xx 



cr a 



a 



a 

—x 
a 



H n \ X -l 



Xx 
l H 



x a . a 

+ x i ~ ~ x i 



fM (x\-H x\)+^H (^-^) +0(H$t) . 



(2.37) 



2 V c Ihq 



Again, in the last line, we have replaced the time derivatives of the scale factor by the Hubble 
rate Hq and length lu at the reference time to = 0. 

The equations of motion in the non-relativistic limit (I2.36P can be integrated directly such that 
we obtain 

x\a — x\a = x\ ao — ^1.0^0 j (2.38) 

where the integration constant was set to the value of the functions evaluated at the reference 
time to = 0. As expected, under a coordinate transformation x\ = ax, this equation coincides 
with the differential equation (|2.10p of the previous subsection. Hence taking in consideration the 
solution for x (|2.1ip and the series expansion in the neighborhood of to f° r the scale factor (|2.5p 
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with p = 1, we obtain the geometrical length (|2.12p for the coordinates x which coincides with 
the coordinate x\ solutions for the path of massive particles 

x l.mass = a x mass = ^rnass = X X.Q + x 1.0 * ~~ 9 ^° X l-0^0 ^ + ^(* 3 ) ■ (2.39) 

In the last equality we have used the coordinate transformation to define the following relations 
x\ = cioXq and x\ = ao(± + x l H ). 

Considering the expansion of the projection factor (|2.33|) and the solution for x\ (|2.39|) we obtain 
the geometrical length for massive particles 

4ass.l = rlo + 4c*+^l.ci 2 + 0(i 3 ) , 
^1.0 = x \.0 = Ceom.l(O) ; 

(2.40) 

^l.c = ^1.0 ~~ x 1.0 Hq = ^geom.l(0)> 
«l.c = 90^1.0^0 = 90 ^geom.l (0)^0 • 

We note that both the velocity and acceleration corrections due to the time dependence of the 
background have the opposite sign of the ones obtained for non-expanding coordinates (|2. 12j) . This 
means that for this coordinate choice, maintaining the same experimental values of the Hubble 
rate Hq and deceleration parameter qo, space would be in decelerated deflation. 

Again we note that these classical quantities correspond to the respective covariant quantities 
evaluated at the initial time t = 0. The covariant velocity components are 

UJj \ ( A,i , ..-Y' / . ,.r' 

a 



7 {x\ + cT m t + cT\ j x{ 



= 7 ( h -x\ + * l ° a ° ~ xUho - %H - x\ 4 (-X (l - -t 

\a a a cr \a J \ a 

a i x\ a - x\ do (2.41) 
~ —x 1 H , 
a a 

v° = ^ = c 7 (i + c 2 r V + crVi) 

f 2 / • \ 2 , 



c 7 1-^i - l--t 



c 2 \a l \ a 
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and the covariant acceleration is 



= 7 (^ + c r Oj v^ + c r oo vO) 

(2.42) 



a i 
--x 
a 



a° = ^ = 7(v + cr 0l v^ + cr° 00 v ) 
« 0, 

where we have used the equations of motion for x\ (|2.36j) and x l (|2.38p to replace the time 
derivatives in these equations and, in the last approximations, we have considered the limit of 
non-relativistic velocities and distances keeping only the lower order terms in the Hubble rate. 
By considering the projection of these quantities into the three-dimensional hyper-surface at some 
fixed space-time point (to,xi_o) it is straight forward to obtain the respective classical quanti- 
ties moon . 

The Hubble law corresponding to the geometrical lengths ^ eom .x = x\(l — (<V a )^) <12.31h is di- 
rectly derived from the expression for the covariant velocity (|2.4ip by projecting it into the three- 
dimensional hyper-surface 



^geom.i (Dx\ ax*Dx^\__ 

ft a x l q 



St \ dr a c dr J \ a c 

®0 d() d CLq ; . » 9 tj' , q o 

— zi.o ^1.0 xittt — v ml {0) +qH L coml t + 0{H t 

a a a a 6 B 



(2.43) 



In the last approximation we have use the definition of the geometrical velocity and length to 
replace the x\ and x\ such that we obtain an extra term of order H 3 . The first term correspond 
to the peculiar velocities discussed in the previous section, a particle moving in the absence of 
any other interactions will decrease its geometrical velocity, v geom ,i = v ge0 m.i(o)( a o/a)- However, 
for expanding coordinates x%, the linear term in the Hubble rate H = d/a is absent, instead we 
have a second order term proportional to the deceleration parameter. 

Let us complete our analysis by computing the light path lengths and radiation frequency-shift. 
From the metric (|2.28p we define the respective line element corresponding to ds 2 = 

c 2 dt 2 = Sij (dx\ - x\ dtj [dx{ ~^4dtj =}► ^ = (ce' + ^-A , (2.44) 

where e l is a unit vector defining the propagation direction for radiation. Under a transformation 
of coordinates x\ = ax this equation coincides with the respective equation for non-expanding 
coordinates x derived in the previous section (I2.16p . Then, multiplying the solution (I2.17p . by the 
scale factor a we obtain the solution for the above equation in expanding coordinates x\ 

1 



4ght.i = *i.o + ( ce * + ^i.o^o) t + - {ce l H - q xl H 2 ) t 2 + 0(t 3 ) . (2.45) 

Where we have also transformed the integration constants taking in consideration the coordinate 
transformation x\ = ax. We note that for light travelling from x\ to the origin of the coordinate 
frame, similirarly to the case of the previous section (|2.19p . we obtain the coordinate distance 
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Axj ightl (t) = x\ — x l hghtl (t), such that considering the projection into the three-dimensional 
hyper-surface (|2.33p we obtain the geometrical light length 

^ = Axl ghtA (l - ±t\ « -ce l t + \ce l H t 2 + 0(H*t 3 ) . (2.46) 

Then, given some travelling time At for radiation emitted from x\.q an d received by the observer 
at the origin of the reference frame, we obtain the following estimative for the light length Ziight.i 
and the geometrical length l geom .i expressed in terms of the fixed background classical length 
estimative In 



l 0(0) 



flight. 1 — y^jj^Hght.l^light.l 



l m -±H l -^ + o(HlP {0) /c) , (2.47) 



Y^j^geom^ge 



geom.l — y Oij tg eom tgeom 



T 2 



= +I m -H ^ + o(H 2 l 2 0) /c) . 

We note that the corrections obtained have the inverse sign of the ones obtained for expanding 
coordinates x in the previous section. 

In order to compute the frequency-shift let us consider the square of four-momentum vector 
k v = (uj/c,k l ), g^uk^k" = u> 2 /c 2 — (ki — (a / a)x\uj / c 2 ) 2 = 0, obtaining the following dispersion 
relation and respective solution for the wave number vector k l 

2 x I li ax\uo\ ( . ,■ dXiLd\ . jj I , a 



u^Sijlckl — cfc? — => k{ = -[ e l --x\ . (2.48) 

\ a c J \ a c I c \ a J 

When deriving the solution it is relevant to note that, for light travelling to the origin of the coor- 
dinate frame, k\ and x\ have opposite signs such that x\.e = —\x\\ and uj = c\k\\ + {a/ o)\x\\uj / c. 
Given an observed wavelength Ai. r at the reception (the origin x\ = 0), the same radiation wave- 
length at some generic point x\ is 

Ai = ^ = • (2.49) 

\ki\ i_«M 

a c 

Then, considering radiation from a source at a coordinate length x\. e corresponding to a received 
wavelength Ai. r = Ai. e (l — (a/a)\x\, e \/c) we obtain the following frequency-shift 

Zl = ^^ = ^-l = -^^-HM + ( q o + l)H 2 ^. (2.50) 
Al.e M.e a c c c z 

Hence, maintaining the numerical values of Hq (jl.ip and go (jl.2p we obtain a blue-shift. 

The results obtained in this section are consistent among themselves, when working with expanding 
coordinates x\ (and assuming the experimentally numerical values for Hq (11. ip and qo (II. 2D ) we 
obtain that space is deflating at a decelerated rate (|2.34p such that light trajectories are decreased 
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with respect to the fixed background approximation (|2.47p and radiation is blue-shifted (|2.50|) . 
These results are the opposite to the results obtained for non-expanding coordinates x in the 
previous section and do not correspond to our physical measurable quantities, instead are it is due 
to the coordinate choices which must be properly accounted for each coordinate system. Next we 
discuss this topic comparing the results obtained for both coordinates choices x and x\ and its 
relation to physical measurable as well as the technical advantage of using expanding coordinates 

2.3 Coordinate Choice(s) and Physical Measurable 

From the derivations in the two previous sections we have just concluded that the specific form 
of the Hubble law depends on the coordinate choice. As the Newton law does (|2.ip . or any other 
physical law. Furthermore, in the framework of General Relativity, instead of the usual spatial 
coordinates, velocity and acceleration, the measurable quantities are the respective geometrical 
quantities (projected to the spatial three-dimensional hyper-surface) which we derived for two 
distinct coordinate choices (non-expanding coordinates x in section [27T1 and expanding coordinates 
in section [2T2]) . As expected, for a given value of the Hubble rate H and deceleration parameter q, in 
each coordinate system we obtain distinct time evolutions for these quantities. Recalling that the 
coordinate choice is a representation of our physical measurements, we can safely conclude that the 
coordinate system for which the (General Relativity) geometrical quantities directly corresponds 
to our world are non-expanding coordinates x. The Hubble law, both for distances (|2. 15 j) as for the 
red-shift (I2,24h has been verified experimentally and is directly described only by the coordinates 
x with metric given by (|2.2p . This conclusion was actually expected since, as already mentioned, 
the metric (|2.2p was originally derived to describe the experimental Hubble law [5]. 

Employing any other coordinate system predicts different expansion rates as we have just exem- 
plified, for expanding coordinates x\ radiation travelling towards an observer is blue-shifted (|2.50p 
which implies that the light path expressed in this coordinate system is shrinking. This result 
is actually consistent, as long as we recall that we have to correctly map the results obtained 
back to the coordinate system which corresponds to our measurable quantities. Depending of 
the quantities being discussed this mapping is achieved either by inverting the simple coordi- 
nate transformation (|2.25p (for scalar and non-covariant spatial coordinate quantities), inverting 
the generalized coordinate transformation (|2.27p -for covariant tensors- or account for the three- 
dimensional projections (for three-dimensional projected quantities). In particular we note that 
the Special Relativity corrections obtained through the factor 7 (a scalar) are directly mapped by 
applying the coordinate transformation (|2.25p to the respective expressions. It is also important 
to remark that this discussion is equivalent to a frame choice. In General Relativity there is no 
preferred frame and one of the cornerstones of the the theory is to be frame independent. There 
is no inconsistency between our discussion and this statement, simply a given observer lives at a 
specific frame and, although the physical laws can be inferred for any other frame, the physical 
measurements are necessarily obtained at the observer frame. 

Depending in the specific problem being solved, it is common to employ coordinate systems 
that do not correspond to our physical measurable. In particular working with expanding co- 
ordinates x\ has a clear technical advantage. These coordinates exactly match the measurable 
three-dimensional projections for the non-expanding coordinates, x\ = ax = I, hence they corre- 
spond to the physical lengths. Consistently the equations of motion for x\ (|2.36p (the Newton law 
for an expanding background) match the respective equations for the geometrical quantities for 
the coordinates x. This result can be explicitly checked by directly considering the non-covariant 
coordinate transformation x = x\/a (|2.25p in the respective equation of motion (12. 8p . Consis- 
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tently also the classical solutions for the x coordinates for the massive particle lengths Z mass (|2,12[) 
coincides with the classical solutions for the massive particle coordinates x\ mass (|2.39|) . to show it 
is enough to consider the coordinate map x = x\ja ()2.25j) including the the mapping of the inte- 
gration constants such that x\ = oqXq and x\ = oqXq + dox l . Consistently with out discussion, 
to compute the physical distances squared from the x\ coordinates (as well as for velocities and 
accelerations) we must employ the Euclidean metric 5ij, we recall that for non-expanding coordi- 
nates the distance is defined as gijX t x^ = a 2 SijX l x^ = 5ijx\x\. We note that the time coordinate is 
the same for both coordinate systems, this ensures that quantities with time-derivatives (such as 
velocities and accelerations) as well as the Special Relativity corrections can be directly mapped 
by the spatial coordinate map. 

Hence the conclusion is that we can perform most of our calculations using expanding coordinates 
x\ such that the physical measurable results are given directly by the coordinate expressions. 
However it is important to note that this mapping is only valid for physical relations (laws) 
that are expressed in terms of the coordinates x\ and its derivatives. For non three-dimensional 
coordinate quantities it is necessary to account for the specific map between them and the coordi- 
nate (or geometrical quantities) transformations in each distinct coordinate system. Specifically 
the red-shift z (|2.24p for non-expanding coordinates x and the blue-shift z\ (|2.50p for expand- 
ing coordinates x\ are written in terms of geometrical wavelengths A (I2.22j) and Ai (12.49h . not 
of coordinate wavelengths. Hence to derive the map between these quantities it is necessary 
to consider a transformation that accounts for the geometrical wavelengths definitions, hence 
A = (Ai/a)(l — (a/a)(\x\\/c)). As for covariant quantities such as the covariant velocity and 
acceleration the mapping is given by the generalized coordinate transformations (12. 27ft . 

In the literature the expanding coordinates x\ are often referred to as physical coordinates, this 
is understood as that the coordinates x\ correspond to the physical observable lengths in our 
specific frame. Given our previous discussion we should not forget that the physical metric cor- 
responds to non-expanding coordinates x (|2.2f) . The generic map between the physical quantities 
corresponding to the coordinates x\, coordinate velocities x\ and coordinate accelerations 'x\ 
and the respective measurable geometrical quantities for the x coordinates as given in (|2.4p . the 
geometrical lengths ^ eom , geometrical velocities v l geom and geometrical accelerations ag eom is 

^phys — — Of — ^geom 

< "fphys = x\ = OX 1 + — (ax 1 ) = fgeom + ~^geom (2-51) 

^phys {aX ) -\- 2 (ax ) — dgeom -1- ^geom ~t~ 2 ^geom 

a a a a 

Once we have settled our world coordinate system (meaning the observer frame) and concluded 
what the measurable physical quantities are, we can compute the estimative for the expansion 
effects in some physical systems for which the spatial scale is much lower than the Hubble length 
x l <C Ih (II. 5h . Based in the results for the x % coordinates obtained for light range measure- 
ments (I2.20p and the red-shift measurements (|2.24p we present some examples in table Q] which 
are clearly below experimental accuracy (the accuracy of Cassini spacecraft measurements [JT] 
are still one order of magnitude below these effects). 

In the introduction we mentioned the existence of the Hubble horizon at the Hubble lenght Ih 
(when considering the observer at the origin of the coordinate frame). A specific proof for the 
existence of an horizon can be obtained from the infinitesimal proper distance for light, ds 2 = = 
c 2 — (f\ — H ri) 2 . From this expression we obtain the following expression for the speed of light 
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System 


Av. Distance 
(meter) 


Range Correction 
iit/-|N jZ.Zvl (meter) 


z 

z |2.24t 


Method Used 


Accuracy 




Reference 


Earth-Sun 


1.50 X 10 11 


+ 1.71 x 10 -4 


1.14 X 10 -15 










Earth-Moon 


3.85 X 10 8 


+ 1.13 X 10 -9 


2.93 X 10" 18 


Laser Range 


Al ra 2.00 X 10" 


" 2 m 


46 


Pioneer 10-Earth 


> 2.85 X 10 12 


> 6.17 x 10" 2 


> 2.17 x 10~ 14 


Dopplcr shift 


Av/v tss 10" 


12 


40 




< 6.60 X 10 12 


< 3.31 X 10" 1 


< 5.02 x 10" 14 










Cassini- Earth 


1.11 X 10 11 


9.44 x 10~ 3 


8.47 x 10 -15 


Dopplcr shift 


A.i//i/ m 10" 


14 


41 



Table 1: Range corrections and cosmological red-shifts for Earth-Sun and some current experiments. Not 
enough accuracy is today achieved to detect expansion effects in these range and Doppler shift measurements. 



in an expanding background 

ri = c — Hr\ 44> f\ = for n = — = Ih ■ (2.52) 

H 

Hence, as expected from an horizon, travelling radiation freezes at r\ = Ih such that information 
cannot be exchange between the two causally disconnected regions, r\ < Ih and rx > Ih- So far we 
have not discussed the large spatial scale Newton law for the expanding background. The higher 
order corrections on the Hubble constant to the classical acceleration that we have neglected for 
the short spatial scale approximation (|2.36|) . will become relevant near the cosmological horizon 
r i ~ Ih- F° r spherical expanding coordinates n we obtain 

h « -c 2 T\ « -qH 2 rx-^rf. (2.53) 
Assuming q < 0, the specific radial distance for which the net effect is null is 

ri = V=ql H . (2.54) 

Taking in consideration the estimative for the deceleration parameter go = —0.589 (jl.2p we obtain 
that for small distances the first term is dominant and expansion has a repulsive effect (massive 
bodies will increase its relative distances) while for large distances the second term is dominant 
and expansion has an attractive effect (massive bodies will decrease its relative distances). Also 
we can obtain a rough estimative for the above value of the radial coordinate in today's universe 
to be n ~ 0.77 Ih- However we note that an observer at the present time located at the origin of 
the coordinate frame is observing distant events (fx ~ Ih) that occurred a long time ago t ~ 1/H. 
This is due to the time that information takes to travel from the event location to the observer 
location. Let us note that we have considered for most purposes the reference time to for which the 
Hubble rate Hq and the deceleration parameter go are evaluated to coincide with the initial time at 
which we measure the initial distances and speeds. Usually the expansion parameters H and q are 
known exactly only for the time t C xp corresponding to their experimental measurement. Taking 
a phenomenological approach we can consider a time series expansion for these parameters such 
that their value at the time to is related to their value at the time t e xp by the relations expressed 
in equation (|B.8P of appendix IB.3I For relatively small scales these corrections are well within the 
measurement error. As an example let us take the galaxy NGC U2885 which, for an observer at the 
solar system, is at a distance of r\ ~ 118 Mpc = 3.63 x 10 m (about r\ ~ Ih/W). Then any event 
observed will have happen at a time to in the past such that to — i e xp = — r i/c = —1.21 x 10 16 s. 
Taking this negative time span we obtain the following relations 

# = #exp (1 + 0.011) , q = g exp (l + 0.027 x (s + 0.178)) . (2.55) 
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These correction are not negligible however they are within the experimental error. For the current 
value of the Hubble rate (jl.ip the error is ±0.04. As for the deceleration parameter an estimative 
for the error is ±0.15 (we will briefly address how to derive this value in the next section f2.4j) 
and the author is unaware of any estimative for the value of the parameter s, here we assume 
that it is of the same or lower order of magnitude than q. Hence in the remaining of this work, 
when computing quantitative estimative, we will assume today's values for the Hubble rate Hq 
and deceleration parameter go given in (jl.ip and (|1.2p . We note that although the universe is in 
accelerated expansion (Hq > and qo < 0) the Hubble rate H = a/a decreases with time. This 
is not inconsistent, it means that the derivative of the scale factor a increases slower with time 
than the scale factor a. As for the time evolution of the deceleration parameter q it decreases 
or increases depending on the value of the parameter sq. Today's value for t/o is negative, hence 
for so > 2qo + 1 it increases (its absolute value decreases), for so = ^Qo ± 1 is constant and 
for so < 2go + 1 it decreases (its absolute value increases). A more rigorous definition of the 
time evolution of these parameters for large spatial scales requires to consider a cosmological 
evolutionary model to account for the large time scale evolution of these quantities. 

Next we will discuss the stress-energy tensor and shortly address how the parameter q can be 
estimated from the equation of state for today's universe. 



2.4 The Stress-Energy Tensor and the Deceleration Parameter q 

In cosmology it is usually assumed that the background corresponds to a homogeneous and 
isotropic perfect fluid. Hence the diagonal components of the stress-energy tensor correspond- 
ing to the gravitational pressure are equal to each other and the mass-energy density is isotropic. 
This tensor is related to the Einstein tensor by the usual relation T^ u = c /{Sir G) which, for 
the FLRW metric in expanding coordinates (|2.28p . is given in equation (|B.6P of appendix[Bj Spa- 
tial isotropy necessarily implies spherical symmetry (the opposite statement is not valid, spherical 
symmetry does not imply spatial isotropy), hence for technical simplification we will use spher- 
ical coordinates in the following discussion. To work in a orthogonal locally Lorentz frame (the 
observer locally flat frame) let us consider a tetrad relating the coordinate metric with the 
Minkowski metric, = r]f l 0e fM fJ- e u u , where the hatted indexes represent the indexes of the flat 
space-time. Specifically for the FLRW metric the non-null tetrad components are 

e d = e \ = l , e\ = -^ , e 2 2 = r , e% = rsm9. (2.56) 

The stress-energy tensor components in the Cartan frame and in the coordinate frame are re- 
lated by this tetrad as = e-^e^T^, where the inverse components of the tetrad eJf are 
straight forwardly computed by raising and lowering the indexes with the respective Minkowski 
and coordinate metrics. Considering the stress-energy tensor for a background perfect fluid in the 
commoving frame as given in equation (|A.8p of appendix |A] with fluid velocity = (c, 0, 0, 0) and 
the Einstein tensor for the FLRW metric (|B.6p . we obtain the following density and pressure 

PH = t -=-^g{{ 

where no summation over the repeated indexes i is implied. We also note that the non-null 
components Tqi = T\q do not correspond to a measurable physical stress nor momentum flux 



3H 2 
8^G ' 




c 2 (l - 2q)H 2 
8VG 



(2.57) 



25 



(shear), their value is due to the coordinate choice, in the Cartan frame Tgj = 0. Consistently, for 
non-expanding coordinates r there is no non-null off-diagonal terms in the coordinate stress-energy 
tensor (the Einstein tensor for non-expanding coordinates given in equation (|B.3P of appendix IB1 
is diagonal). 

These equations allow us to estimate today's deceleration factor q = — (do/ao)^ as given in 
equation (jl.2p . The combined WMAP+BAO+SN data [H] gives a cosmological constant (dark 
energy) relative density of Q\ = 0.726 ± 0.015 with an equation of state p\ = c 2 u>aPa, where 
uj\ = — llto'14 ■ Assuming that the background is a perfect fluid, that the remaining matter is 
pressureless (u m = p rn = c 2 uj m p m = 0) and neglecting the radiation contributions (p r « 
and p r w 0) we obtain the universe equation of state 

2 

PH = c louPh , u H « = -0.726 , (2.58) 

\l H 

where the total relative density is unity by definition, = 1. Then, taking the expressions for 
the FLRW pressure and density (|2.57p . and solving the equation of state for qo we obtain 

q « i + ^u;t = -0.589 . (2.59) 

This derivation can be found, for instance, in section 29 of [37] . The error associated to this 
estimative can be inferred directly from the values for f^A and uj\ being the upper error bar 
+0.167 and the lower error bar —0.146. However these values have no statistical significance, 
the value for u\ is quoted to 95% confidence level while Q\ is quoted to 68% confidence level. 
Also we note that this value for qo was computed assuming that the FLRW metric (|2.28p does 
describe the universe today and within the ACDM model which takes into account baryonic matter 
(the usual heavy matter), dark energy (the cosmological constant) and cold dark matter (missing 
non-interacting or weakly interacting matter). 



2.5 Observer not at the Origin of the Coordinate Frame 

So far we have always considered the observer at the origin of the coordinate frame. For this case 
the geodesic equation for the observer is trivial, ri = at all times. The origin of the coordinate 
frame is the only spatial point not affect by expansion. Generally we can consider a commoving 
observer at some generic spatial point ri. b s 7^ 0. The equations of motion for such an observer 
are given by the geodesic equations which have the solutions (|2.34[) . All space is expanding and 
the observer necessarily is attached to it (see for instance the reviews [H E] for a discussion in 
this topic). Hence a commoving observer at a generic spatial observing an event at some other 
spatial point ri will not perceive the geodesic path of the event, instead will be measuring the 
differences between its own geodesic path and the one of the event. Then, from the perspective 
of a commoving observer, instead of the geodesic equations it would be the geodesic deviation 
equations that describe the evolution of the event 

d 2 (x{ - x s l ohs ) _ s p p d(x^ - x^ obs ) d{x\ - %i obs ) . . 

^2 " K W&l *l.obsJ ^ ■ 1 

These equations are valid when the geodesic separation £^ — {%\ — ■^l'obs) ^ ^udi lower than 
the characteristic lengths of the physical system being studied. For the expansion effects these 
are the Hubble time tn ~ 10 17 s (|1.4|) and the Hubble length ljj ~ 10 26 m (|1.5|) . Hence for 
most astrophysical observations the geodesic deviation equations, in the absence of any other 
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interactions, are valid. To lowest order in H in the non-relativistic limit (fi <C c) and assuming 
only radial motion {6 = (p = 0) these equations are 

* ~ *obs , h ~ -qH 2 (n - r Lobs ) , (2.61) 

hence matching the lowest order equations of motion for the expanding coordinates (|2.36[) which 
correspond to the Newtonian limit of General Relativity. The relativistic corrections are substan- 
tially distinct from the expressions (|2.37|) obtained in the previous sections depending both in the 
geodesic separation e M and the observer coordinates x^ ohs such that an observer at ri will have to 
perform simultaneous measurement of its own position with respect to the origin of the coordinate 
frame, ri, bs> as well as of the distance to the event e. 

We will not develop in full detail these calculations here, it is relevant to mention this issue to 
alert for the need to correct the observables due to the geodesic of the observer as well as to recall 
that only for geodesic separation much lower than the characteristic length of the physical system 
it is a valid approach to consider the geodesic deviation equation. A possible alternative to this 
construction, as has been carried out in the previous sections, is to consider a fixed observer at the 
origin of the coordinate frame which has a trivial geodesic path such that the relative motion of 
the event being observed is directly given by its geodesic path. When considering matter effects, 
depending on the physical measurements we want to estimate, this issue may be relevant. However 
for the main results presented in this work this discussion plays no role. 

Next we review how matter in flat an expanding backgrounds is usually described. 



3 Reviewing Matter in an Expanding Background 

In this section we review how the gravitational interactions of matter are described at the level of 
the metric. We start by resuming the relevant characteristics of the Schwarzschild metric [1] for this 
work which, within the framework of General Relativity, is considered to be an exact description for 
space-time in the neighborhood of a point-like mass (central mass) in a flat background (Minkowski 
space-time). We also resume the existing metrics and the main characteristics of the respective 
space-times for matter in an expanding background given by the FLRW metric (|2.28j) , namely the 
McVittie metric [3] and the Cosmological-Schwarzschild metric [6]. 



3.1 Schwarzschild Metric 



The gravitational potential for a central mass (a point-like massive body) has spherical symmetry 
with respect to the center of mass, hence both for technical simplification and easier physical 
interpretation, it is often assume that the mass is centered at the origin of the coordinate frame 
and employed spherical coordinates. Specifically we will use the following definition for the di- 
mensionless Schwarzschild gravitational potential Use 

U sc = -^. (3.1) 

Here M is the gravitational mass of the massive body, G is the Newton gravitational constant 
and c is the speed of light. Semi-classically it is often considered the dimensionfull gravitational 
potential i^sc = c2 ^Sc/2- The Schwarzschild (SC) metric |4J for such a point-like mass is written 
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in terms of the gravitational potential Use as 
2GM\ 

C 2 T\ ) 



ds 2 = ^fi_^^__^__ r f(dff2 + 8m 2 ^) 



2 / 2GM 



c z | 1 == | (it 

c 2 y5ijx\x{ 



2 



( \ (3-2) 

2GM5 ik 6 jl x\x[ 



5u + 



dx\dx\ 



where we wrote the infinitesimal length square element both for spherical coordinates (for which 
the metric is explicitly spherically symmetric) and Cartesian coordinates. No effects for space 
expansion are yet considered, however let us consider expanding coordinates r\ which correspond 
to the physical measurable distances as discussed in section [231 The Ricci scalar (scalar curvature) 
and the curvature invariant are 

^sc = -Svr ^*W(n) , ^sc = Rap***** = , (3.3) 

C C T-y 

where 5^(r\) stands for the three-dimensional Dirac-delta. For Cartesian coordinates it is the 
product of three one-dimensional Dirac-delta 5^((x,y,z)) = 8(x)5(y)5(z), while for spherical 
coordinates it is given only by a one-dimensional Dirac-delta on the radial coordinate divided by 
the area of the sphere such that its volume integral is normalized to unity 



£(3) (ri) = ^ d6 d(p dr 1 ^S^(r 1 ) = l . (3.4) 



Therefore the SC metric has a singularity at r\ = which corresponds to the mass pole at the 
center of mass, the origin of the coordinate frame. The respective mass density can be computed by 
considering the Einstein tensor components written in the Cartan frame. For the SC metric (|3.2p 



the Cartan Tetrad is = yjg^J (without summation over repeated indexes implied) and the 
Ricci tensor components are R^y = R^y/g^u = 4:7r(GM/c 2 )5^ (ri) (also without summation over 
repeated indexes implied). Assuming the stress-energy tensor for a perfect fluid we obtain the 
mass density from the Tgg component (|A.8P 



c 4 



T 06 C = ^g G 66 =c 2 M^)( ri ) ^ p sc = M5^(r 1 ). (3.5) 



Hence, for the SC metric, M is both the gravitational mass (|3.ip . the mass pole value and the 
total mass of the space-time which is obtained by integrating the mass-energy density psc (j3-5j) . 

There is a relevant remark to the remaining of this work, to have a mass pole at the center of mass 
(here at the origin) the curvature invariant 1Z (|3.3|) must be at most divergent by the inverse square 
of the volume element, i.e 1Z ~ V~ 2 ~ r^ 6 . This requirement implies that the '00' component of 
the stress-energy tensor is at most divergent by the volume element, T Q0 ~ rj~ 3 , such that the pole 
value is finite and can be obtained by considering the volume integral of the mass-energy density 
inside a shell of fixed radius r\ and taking the limit r\ — > 0. For more severe divergences, 1Z ~ r _n 
with n > 6, we obtain a essential singularity which no longer can be interpreted as a mass pole. 
The volume integral of the mass-energy density inside a shell of radius r\ is divergent (as well as 
its limit n — ► 0). 

With respect to the metric components for the Schwarzschild metric (|3.2p are singular at the 
Schwarzschild radius, the well known coordinate singularity, 

2GM (i «\ 

n.sc = — n— , (3.6) 
c z 
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which is interpreted as an event horizon for an external observer. 

For astrophysical systems, in order to work perturbatively in the gravitational field, it is often 
considered the (gravitational) weak field approximation 

2GM , . 

C/sc = -g— « 1 • 3.7 



This inequality is obeyed for distances much larger than the Schwarzschild radius (|3.6p . Taking 
the solar system as an example and considering the gravitational field at the surface of the sun, 
corresponding to the mass Mg un ~ 1.98 x 10 30 Kg and radius rs U n ~ 6.96 x 10 s m, we obtain 
U ~ 10~ 6 <C 1 (for further example see box 19.2 of |37j). 

With respect to the proper time rsc an d proper lengths /i.sc> they do not generally coincide 
with the coordinate time t and coordinate distance r\. However for non-relativistic velocities 
ri,d,(p <C c and weak gravitational fields (|3.7p they are approximately the same. Considering an 
observer at ri. b S) far enough of the mass M, observing events near the mass at r\ ~ such that 
?"i.obs 3> ri, we obtain 



2GM 1 , . 

t « t , n = — ^^^=^^^= r\ « n . (3-8) 



V Tl.obs-nl Jl 2GM 

V kLobs-^ll 

Again, this is a valid approximation for most astrophysical measurements. 

As for the classical Newton law for a test particle of mass m and the respective radial acceleration 
are as usual obtained directly from the gravitational potential 



c 2 GM 

F N ewton = m — VU ^ 7=1 = 5- . (3.9) 

2 r z 



In the Newtonian limit the General Relativity corrections to the Newton law can be obtained 
directly from the geodesic equations. Neglecting the corrections due to the velocities of the test 
mass the relevant equation is the radial component of the geodesic equations, specifically we obtain 

1 GM 2(GM) 2 
fi « -c 2 T\ = + \ 3 > . 3.10 



The first term coincides with the classical Newton acceleration (|3.9p . The second term is due 
to General Relativity correction and, outside the Schwarzschild radius ri se = 2GM/c 2 (j3.6|) . is 
never dominant being its value always below the first term. At the Schwarzschild radius (the 
event horizon) both terms have the same value such that the net gravitational acceleration is null. 
Either when the test mass is relatively close to the central mass or when a higher accuracy for 
theoretical estimative is required the geodesic deviation equations are employed instead of the 
geodesic equations. 

Being relevant to the present work, we next discuss what is understood by local anisotropy in the 
presence of a massive body. Light travelling in a gravitational potential generated by a massive 
object is either blue-shifted or red-shifted depending on weather it is travelling away or to the mass, 
respectively, while if travelling (approximately) tangentially to the radial direction its frequency 
is not affected. This is only an apparent spatial anisotropy, although the speed of light is not 
constant along all directions, space-time is still isotropic. Specifically from the SC metric (|3.2[) we 
obtain directly from the infinitesimal line element for travelling radiation the following expressions 
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for the speed of light 



dri _ ( 2GM 

dt \ n c 2 



d6 dip / 2GM 

ce = c ip = n— = nsmd— = c\l T 

dt at V T\ c z 



(3.11) 



Nevertheless there is a specific coordinate choice r 2 for which the speed of light is isotropic. These 
coordinates are usually called isotropic coordinates and the coordinate transformation from the 
coordinates r\ to the coordinates r 2 is given in equation (|A.14|) of appendix |Aj For the new radial 
coordinate r 2 we obtain the metric 



(l _ GM \ 4 

ds 2 = c 2 A ^H-dt 2 - ( 1 + ^L) (dr 2 + rldO 2 - r 2 sin 2 pdp 2 ) , (3.12) 

{ i i GM \ \ 2r 2 c l ) 



1 + 277^ 

for which the speed of light is constant along all spatial directions 

1 _ 2 GM 

dr 2 1 \r 2 -r 2 . \c^ , s 

Cr = Ce = c = — = c _ . (3.13) 

Ml 2 GM y 

Although not directly corresponding to the physical coordinates (the Schwarzschild coordinates) 
these are widely employed in astrophysics, in particular in the PPN formalism [42] . 

A general definition of spatial isotropy is with respect to mass-energy and momentum fluxes which 
are encoded in the energy-momentum tensor T^. Isotropy of space implies that, in the Lorentz 
frame with Minkowski metric rj^ related to the coordinate metric by = e^e^T^p through a 
tetrad the mass-energy flux is null, Taj = 0, and the momentum fluxes are the same (isotropic) 
along the spatial dimensions being identified with the matter pressure p, TV. = S^p. For the case 
of the SC metric (|3.2p the energy- momentum tensor is identically null except for the mass pole at 
the origin (13. 5ft . hence the matter density and pressure are null everywhere else, p = p = 0, which 
is consistent with the assumption of space being empty of any sort of matter. 



3.2 Planetary Orbits 

Here we shortly review how to derive the orbit solutions for planetary motion around a central 
mass and the theoretical orbital properties both for Keplerian and General Relativistic orbits, 
namely we discuss orbital precessions and periods. 

We will take the usual approach for two body orbital systems by considering a variational func- 
tional given by the Lagrangian 

C dx^dx^ 2 

— = Sixv—, -j— = c , (3.14) 

171 (IT CLT 

where m is the orbiting mass which factors out from the Lagrangian (see for instance section 8 
of [38] or section 25.1 of [37] ). For the Schwarzschild metric (|3.2p . considering an orbit lying in 
the plane of constant coordinate 9 = ir/2 such that d9 = and sin# = 1, we obtain the following 
equation 

£SC " (1 " ^sc) ( c #V - ^7- ( P-X ~ * (f) 2 = <? ■ (3.15) 



m \ dr J 1 — Use \dr J \ dr J 
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The Lagrangian does not depend explicitly either in the time coordinate t nor on the angular 
coordinate <p, hence there are two constant of motions. Respectively E and J defined as 

2E 1 5C . TT . / dt 

2(1 - U SC ) [c 



mc m 5(cdt/dr) \ dr 

(3.16) 

23 1 SC ^ 2 / dip 

m d(d<p/dr) 1 \dr 

We note that these constants of motion correspond to angular momentum conservation (J) and 
energy conservation (E). Replacing both of them in the Lagrangian (13.150 . multiplying by a factor 
of (1 — Use) ^ gathering the constant terms on the left-hand side of the equation, considering the 
derivative transformation from proper time r to the angular coordinate ip and redefining the radial 
coordinate to u = 1/n such that 

dri(r) dr^tp) dipjr) = r[(ip) 

' u VP) = 1 TT = 2TT ' ( 3 - 17 ) 



u(ip) 



1 dpn(ip) rf(ip) 



ri(<p) 

we obtain the following differential equation on the function u = u(ip) 



c Ml_( M ) =_J2( U ')2_ j2 u 2 + 2GMu + J 2^GM u 3 (318) 

mc z J I c z 



In the variable and function transformations (j3. 17h we have explicitly written the dependence 
on proper time r and <p while in the differential equation (|3.18p we have omitted the variable 
dependence. The first three terms in the left-hand side of this equation correspond to the Keplerian 
orbit equation while the last term is the General Relativity correction. 

The differential equation (|3.18p is non-linear, the usual linear second order orbital equation is 
obtained by differentiating with respect to the variable (p and factoring out an overall factor of 
2v! such that we obtain 

GM 3GM 2 
u + u =—p~-\ u ■ (3.19) 

Again, the first term on the right-hand side of the equation match the classical Keplerian orbits, 
while the second term is the General Relativity correction which has as effect a small perihelion 
advance with respect to classical solutions (which is equivalent to an orbital preccession effect). 

The classical (Keplerian) orbital ellipse solution is 

GM 

< + u = -j r , (3.20) 

with solution given by 

l + ecosv9 

n = , (3.21) 

which corresponds to an ellipse with the mass M on one of its foci where e € [0, 1[ is the orbit 
eccentricity. As for the parameter d and the constant of motion J 2 are related to the semi-major 
orbit axis r l or b as 

d = r Lorb (l - e 2 ) & J 2 = r Lorb GM (1 - e 2 ) . (3.22) 

We note that both e and ri. or b must be set as initial boundary conditions for each orbit and are 
not derivable theoretically, however the constants of motion J and E are not independent. J 
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is set from the relation between the differential equation (|3.20p and the physical (geometrical) 
interpretation of the parameter d (|3.22p . As for the constant of motion E can be approximated 
for a classical orbit by replacing the solution no (|3.21[) in the non-linear differential equation (|3. 18j) 
such that neglecting the last term on the right-hand side, we obtain 

E 2 « El = (m c 2 ) 2 ( 1 - -^s) ■ (3.23) 

V r l.orb C 2 J 

The period To of this classical solution is known to be independent of the orbit eccentricity e being 
given by Kepler's third law 

3 

2-7T r 2 

T = . (3.24) 

Vgm 

The General Relativity corrections are obtained by considering the differential equation (|3.19j) 
with the last term evaluated with the classical solution no (13.211) 

n = n + n G R , 

2 (3.25) 



U GR + U GR 



3GM f 1 + e cos tp 



c 2 V d 



We note that this is a valid approximation due to the GR correction term being much smaller 
than the classical term (assuming large radial coordinate, hence r > 1 and u = 1/r <C 1). The 
solution of the differential equation (|3.25p is 

UGR = (( X + y) ~ -^ cos2i P + ei P shnp\ , (3.26) 

where we have set the integration constants so that no simple oscillating terms are present con- 
taining either cos tp or sin (p. The reason to do so is simply to maintain the same values for the 
eccentricity e and parameter d than in the unperturbed solution (|3.2ip . Specifically we could have 
set the integration constants such that a term e&GR cosy? with /cqr = 3GM/(c 2 <i 2 )(l + e 2 /2) is 
present which maintains the same value for the orbit eccentricity. Then together with the first 
constant terms in the above solution (|3.26p . it could be included in the unperturbed solution no 
by re-defining the parameter d by a small shift 



1 + &GR ^ 3GM ( e 2 

u = n + nGR = ^ (1 + e cos ip) + ^ ^ I — — cos 2(p + e ip sin tp 

d « 



(3.27) 



Here we have used the notation ncR to distinguish this solution from ncR (|3.26p . For this exam- 
ple this corresponds to a redefinition of the semi-major axis t\ or b by the amount /cgr- A similar 
construction can be carried for which either the eccentricity e or both e and fi.orb are corrected. 
However we note that these correction, besides being negligible when compared with the remain- 
ing terms, have no particular physical meaning, we recall that both the orbital eccentricity e and 
semi-major axis ri. or b are set as boundary conditions (or initial value) for the differential equa- 
tions, hence we are simply re-defining our equation parameters and its relation to the physical 
(geometrical) quantities of the physical orbit. 

Hence, with respect to the several terms in the correction (I3.26p . we have that the first term is a 
constant that, following the previous discussion, can be neglected, the second term is a periodic 
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deformation of the ellipse (the orbit) that slightly increases the orbital path contributing a small 
deviation to the period with respect to the classical one, To (|3.24|) . and the last term increases 
steadily over time (here with the angular coordinate) and is responsible for the well known orbital 
precession (the ellipse rotates steadily around the center of the mass M). To analyze the effect of 
these terms let us define the correction parameter 

3GM 

&GR = c2 ^ , (3.28) 

and consider the following approximation to the sum of the term cos <p on the solution uo (|3,21l) 
and the term ogr V srn <P m the solution ugr (I3,26h by considering the respective lower order 
series expansion 

(f 2 if 2 if 2 

cos (p + otQR ip sin tp sa l- — + 2a GR — rs 1 - (1 - 2a GR ) — 

1 1 1 (3.29) 



pa cos(Vl - 2a GR ip) pa cos((l - a GR )(f). 

From the first to the second line we have approximate the series expansion by the cosine func- 
tion and in the last equality we have expanded the square root to lowest order in the correction 
parameter «gr- Hence we obtain the following approximate solution of the full differential equa- 
tion flgjgH 

u pa - (1 + e cos ((1 - a GR )(p)) + u osc .gk , 

(3.30) 

a GR 2 / \ 
bd 

The precession amount per turn of the orbit is directly taken from the argument of the cosine 

A(£GR 



2vr 



otGR , (3.31) 



where we are using radians to measure angles. As for the period correction it can be evaluated 
from the term m sc.gr by noting that the infinitesimal time displacement dt is given in terms of 
the tangential velocity to the orbital path, v± = dx±/dt, as 

dt = — — = dip , 

v± uv± 

(3.32) 

dip I dip f dtY 1 I 2GM 

where v±, n = 1/u and u are functions of the angular coordinate ip. In the first line we have 
used the relation between the infinitesimal spatial displacement on the orbit path dx± and the 
infinitesimal angular displacement dip, dx± = r\ dip = dip/u. In the last equality we have used the 
definitions of the constant of motion J (|3.16p and considered the limit of non-relativistic velocities 
i^Cc such that {dt/dr) -1 = 7 -1 ~ yT— 2GM u/c 2 . Hence, considering the integration of the 
infinitesimal time displacement dt (|3.32[) over one turn of the orbit, we obtain, to lowest order, 
the period correction due to the term u sc.GR 

T 

dt = T = T + AT GR , 

(3.33) 

- Wosc.GR /, GM 







ATgr - J q *p-q- (1-— wo 
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To derive the final expression for ATqr we have expanded the expression for dt (|3.32p with 
u uq + iiosc.GR to first order in u sc.GR and integrated the term containing this oscillatory 
correction. The term containing only the solution uq corresponds to the period To. Later on, 
when discussing both the precession and period corrections due to the expanding background, for 
comparative purposes of the magnitude of the several effects we will give numerical estimative for 
equations (|3.31j) and (13.330 . 

The remaining correction to the orbital motion are due to solar oblateness (Sun quadrupole mo- 
ment J2) which contributes significantly to the orbit precession (see for instance section 40.5, 
box 40.3, of [37]) and the interactions between the several celestial bodies. It is only possible to 
estimate the ephemerides of the planets by the use of extensive numerical calculations in the PPN 
formalism [32] including most of the known bodies in the solar system, see for instance |43] for 
further details. We will not give further details on these topics and will compare the estimative 
obtained when considering an expanding background with the General Relativity corrections for 
flat backgrounds discussed previously. We note that the precession effects due to General Relativ- 
ity are usually significant and are taken in consideration in theoretical calculations, however the 
orbital period corrections are usually negligible when compared with the other corrections due, 
mostly, to the interactions with other bodies, hence dealt only numerically. We considered these 
corrections here for comparison purposes only. 

Next we briefly resume the existing metrics to describe matter in an expanding background. 
3.3 McVittie Metric 

The McVittie metric [3] was derived originally in non-expanding isotropic coordinates r% demand- 
ing both spherical symmetry and spatial isotropy such that the stress-energy tensor is shear free. 
For this coordinate choice this metric is given in equation (jC.lj) of appendix O For expanding 
coordinates r\ (|C.3P it reads 



where we explicitly factorize the length square infinitesimal element into proper-time and proper- 
distance. This metric, for expanding isotropic coordinates r<i and non-expanding coordinates r, is 
given respectively in equations (1G2|) and (IC.4j) of appendix O 

The McVittie metric is a sobreposition of the Schwarzschild metric f)3.2[) describing small spatial 
scales with the FLRW metric (]2.2p describing large spatial scales. As can directly be inferred from 
the metric expression (|3.34p it has the following properties: 

1. asymptotically, it coincides at spatial infinity (t\ — ► 00) and in massless limit (M — > 0) with 
the cosmological FLRW metric ([2.280 : 

2. in the static limit (a — * 1) it coincides with the Schwarzschild metric (|3.2D ; 

3. for expanding coordinates r\ it has Lorentzian measure which corresponds to the spherical 
coordinate measure to be \J—g = r\ sin#. 




In addition it is isotropic as well as stress and shear free. The non-null spatial components of 
the stress-energy tensor in the orthonormal Lorentz frame (Cartan frame) are 7« = Ta/ga = 





(3.34) 
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Tqi/<7oi = P (for z = 1, 2, 3), where no summation over repeated indexes is implied. Later we will 
write the explicit expressions for this tensor as a particular case of a more generic metric. The 
non-null off-diagonal component Tq% is not a measurable physical stress, it is due to the coordi- 
nate choice as we have discussed earlier in section 12.51 for the FLRW metric (|2.28p in expanding 
coordinates r±, specifically in the Lorentz frame we have that Tgj = 0. 

If we insist in that space isotropy must be maintained it is accepted that there is no other possible 
choice [El[9]. 

At the origin (n = 0) the McVittie metric has the same pole of the Schwarzschild metric, the 
dominant singular term of the curvature invariant in the neighborhood of the origin is 

^Mcv(n ~ 0) ~ 48( ff 2 . (3.35) 

In addition it has also a singularity at the Schwarzschild radius r\ = rgc = IGMjc? [7], in its 
neighborhood the dominant divergent term of the curvature invariant is 

12(l-g) 2 ff 4 

7<-McV(>l ~ ri.SCj r 2GM • ^ 3 - 36 - ) 

This is an extended singularity corresponding to a two sphere. Hence, although in the static limit 
a — > 1 the McVittie metric converges to the Schwarzschild metric, in the limit r\ — > 2GM/c 2 
(near the Schwarzschild radius) the McVittie metric (|3.34p does not converge asymptotically to 
the SC metric (13. 2D having a clearly distinct behavior due to the singularity at the horizon. We 
note that it is reasonable to expect that a metric describing both expansion and matter effects 
asymptotically converges to the SC metric at the SC horizon, in this limit the matter effects should 
be dominant with respect to the background fluid effects describing expansion. Also within the 
framework of General Relativity, the SC metric has been extensively and successfully employed 
for different coordinate choices being in agreement with experimental data for most astrophysical 
systems (see, for instance, section 38 of [37J). From a more theoretical perspective this metric 
describes a non-complete space-time and the total mass inside a shell of finite radius r\ > ri.sc is 
divergent due to the divergence at the SC radius ri.sc- These are clearly unwelcome properties. 

We recall that the original derivation of this metric [3j has been carried for isotropic coordinates 
which are employed for large radial coordinate (usually in astrophysical systems), in particular 
the map between these coordinates and the usual expanding Schwarzschild coordinates (|A.14|) is 
only defined outside the Schwarzschild horizon {r\ > 2GM/c 2 and T2 > GM/{2c 2 )), hence no 
analysis was carried for the asymptotic limit near this horizon. 



3.4 The Cosmological-Schwarzschild Anisotropic Metric 

Global space isotropy is commonly accepted as a fact, both due to theoretical reasoning, Poincare 
invariance is only strictly maintained for isotropic space-times, as well as due to large-scale ob- 
servations of our universe, the background radiation in our universe is globally isotropic |14| . 
However there is experimental evidence for local anisotropy corroborated both by local deviations 
of the Hubble flow in nearby astrophysical systems [13], as well as due to the local anisotropics of 
the background radiation |14j . 

Hence it is not physically unconceivable that a description of matter in an expanding background 
may generate a local spatial anisotropy. Assuming that this may be the case there is one metric 
that has already been considered by several authors employing distinct technical approaches [6] . It 
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can be justified by noting that we have been measuring spatial lengths without taking in account 
spatial expansion. Then, being the classical Newton law our starting point, for non-expanding 
coordinates r, we can rewrite it taking in consideration the spatial expansion as 

GM 

ar = — . (3.37) 
a z r z 

The extra factor of a multiplying r is due to the projection to the three-dimensional hyper-surface 
as discussed in section [2j also the acceleration expands with space. 

The metric corresponding to the above Newton law (|3.37p . both for non-expanding coordinates r 
and for expanding coordinates r\, is 

2GM\ 2 2 2 2 ( r/i2 , • 2 n i 2\ a 2 dr 2 



ds * = 1 - C J d t z _ a J r ; (ffl* + s in 2 6dy? 



ac 2 r 1 v ' \ — 2 ^M 

' ac z r 

2GM\ 



c 2 r\ 



c 2 dt 2 - r\ (d6 2 + sin 2 Odtp 2 ) ( 338 ) 



1 f+.tud™},* 



I _ 2GM_ \ ac \ c 2 ri 



■ 



Where we explicitly factorize the infinitesimal length square into proper-time and proper-distance. 
We remark that by considering the Schwarzschild metric expression (|3.2p to be written in non- 
expanding coordinates r, this metric is obtained by a direct replacing of the non-expanding radial 
coordinate r with an expanding radial coordinate r->ar and correcting the radial component of 
the metric accordingly dr — > a dr. Due to this relation we refer to this metric as Cosmological- 
Schwarzschild metric (CSC). 

It has the same three properties listed for the McVittie metric maintaining spherical symmetry, 
however space-time is locally anisotropic, the components of the stress-energy tensor in the Cartan 
frame obey the following relations Tjj = T u /g n = T 01 /g 01 ^ T 22 = = T22/522 = T 33 /g 33 , 
hence the radial component and angular components have distinct values which explicitly shows 
the existence of spatial anisotropy. 

As the McVittie metric, also this metric describe space-time with two singularities, at the origin 
of the coordinate frame r\ = and at SC radius r\ = 2GM/c 2 . In the neighborhood of the origin 
the curvature invariant dominant term coincides with the Schwarzschild curvature invariant 



48(GM) 2 

fts(ri ~ 0) VlT^ • ( 3 - 39 ) 



1 

As for the neighborhood of the Schwarzschild radius ri.sc = 2GM/c 2 (|3.6p the dominant term in 
the curvature invariant is 



4(GM) 2 H\l+qf 





Ks(ri ~ n.sc) r 7-^2 — ; • ( 3 - 40 ) 

cr l (\ 2GM \ 



Hence it has the same unwelcome properties of the McVittie metric, does not converge asymp- 
totically to the SC metric near the central mass [r\ ~ n.sc); space-time is not complete and the 
mass inside a shell of finite fixed radius is divergent due to the singularity at the SC radius. 

Next we will consider a more generic metric ansatz(e) that has, as particular cases, both the 
McVittie metric ()3.34p and the CSC metric (|3.38[) . We will also obtain a parameter range for 
which only the mass pole at the origin is present such that the SC horizon is singularity free and 
the SC metric ([3.2p is asymptotically obtained at the SC radius. 
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4 A Locally Anisotropic Metric for Matter in Expanding Space- 
Time 

Once we consider a background fluid, as for the case of the FLRW metric (|2.2p . are generally 
present a non-null background gravitational density and pressure. It is a very conservative and 
well accepted assumption to expected that, by considering a point-like massive object in such 
a background, we obtain a spherically symmetric pattern deformation asymptotically vanishing 
at spatial infinity. We also remark that for a central mass, the direction of the gravitational 
interaction is radial, hence there is a preferred spatial direction. Hence it is not unconceivable, 
even being physically intuitive, to consider the fluid deformation to be anisotropic with respect to 
the radial and angular directions. Also, as already discussed, local anisotropy is consistent with 
experimental observations [I3j[l4], as long as global spatial isotropy is preserved. 

In this section we will build a metric ansatz describing matter in an expanding background that 
generalizes the McV metric (|3.34p and the CSC metric (|3,38p and interpolates between the FLRW 
metric (12.28P and the Schwarzschild metric (13. 2p maintaining space-time free of singularities except 
for the Schwarzschild mass-pole at the origin. Although, for the metric ansatz, space-time is 
locally anisotropic, spatial isotropy is recovered at spatial infinity. We start by the simpler case 
of a one parameter metric analyzing the singularities for the several values of the parameter such 
that the space-time is regular at the Schwarzschild radius, hence obtaining a complete space-time 
asymptotically flat at the SC horizon. However the singularity at the origin is more severe than the 
SC mass-pole such that the mass inside a shell of finite radius is divergent. To strictly maintain 
the SC mass-pole at the origin we refine the ansatz by considering one further regularization 
parameter obtaining a space-time with finite total mass inside a shell of finite radius. We further 
analyze the stress-energy tensor defining the range of the parameters for which the mass-energy 
is positive outside the event horizon. 

4.1 The Ansatz I: A First Approach 

In the following we consider that the expansion of space is a global effect due to the background 
matter and energy spread across all universe such that cannot be locally eliminated. This is not 
a widely accepted assumption, except for the cosmological constant effect (dark energy) which 
accounts for 72.6% of all background gravitational effects [13] and, due to be a constant, cannot 
be physically excluded from any system. A common argument concerning the remaining 27.4% of 
matter and energy is that a stationary low density dust (matter background) in the neighborhood 
of a massive stellar object is attracted and rapidly aggregated by that massive object. Nevertheless 
we note that, generally, also stars and other massive objects are responsible for matter and radia- 
tion emissions (hence also contributing for the background matter and energy density), moreover 
expansion is a global effect mostly due to long range gravitational interactions, the matter and 
energy density far from a massive object also contribute to the local expansion and, in the same 
fashion are also affected both by the gravitational field of the several local massive objects and 
their long range interactions. Hence our assumption is that there is no natural mechanism to set 
a local cutoff for which expansion is locally eliminated. 

We next proceed to generalize both the McV metric (13.341) and the SCS metric (13.380 . We 
consider expanding spherical coordinates r\ and will build an ansatz intended to describe local 
matter in an expanding background interpolating between the FLRW metric (I2.28P which describes 
the cosmological expanding background and the Schwarzschild metric f|3.2|) which describes local 
matter in a flat background. With respect to the metric properties, when setting up this ansatz, 
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we consider the following assumptions: 



1. in the massless limit M — > and at spatial infinity converges asymptotically to the FLRW 
metric (^281) : 

2. in the static limit a — ► 1 (d — ► 0) coincide with the Schwarzschild metric (|3.2p ; 



3. has Lorentzian measure ^J—g = 1 for Cartesian coordinates x\ which corresponds to the 




spherical coordinate measure \f—g = r\ sin such that the area of the sphere is A = 4ir r\ . 




These characteristics are shared by both the McV metric (|3.34p and the CSC metric (|3.38p dis- 
cussed in the previous section. The first two assumptions are common to both these metrics, 
we wish to interpolate between two distinct metrics (the FLRW and SC metrics) and these as- 
sumptions ensure that in the limiting cases we retrieve the two original metrics. The second 
assumption can become stronger, we could demand the ansatz to converge asymptotically to the 
SC metric at the origin. This case is actually relevant to strictly maintain the SC mass pole at the 
center of mass and will be addressed in the next section. In addition we will also obtain bounds 
for the metric parameter for which the ansatz converges asymptotically also to the SC metric 
at the event horizon (the SC radius), hence for which space-time is asymptotically flat near the 
point-like mass. In order to both maintain some anchorage with already existing results, namely 
the McVittie metric (]3,34j) and CSC metric (j3.38[) . and to keep track of the several steps in the 
ansatz(e) building and the respective physical interpretations, we will proceed with this milder 
requirement. 

As for the last assumption it may, generally, be lifted and it is clearly coordinate dependent. 
Nevertheless for expanding coordinates r\ (which correspond to the physical measurable lengths 
as discussed in section I2.3|) it is expected that we obtain the Lorentz measure which is shared 
by both the FLRW metric (|2.28[> and the Schwarzschild metric (j3.2|) . This is not by chance and 
it reflects the way we perceive our world, in particular with the way we perform measurements 
and with Local Lorentz invariance. Specifically Euclidean geometry is the basis of all our spatial 
measurements such that our measurable spatial sphere has area of A = An r\ and the geometrical 
generalizations to four-dimensional space-time manifolds are necessarily Locally Lorentz, at least 
in a patch including the physical system being studied. In simple terms this means that a space-like 
vector remains space-like and a time-like vector remains time-like being each orthogonal to each 
other such that an orthonormal Lorentz basis (corresponding to the Cartan frame with Minkowski 
metric) can be considered locally at each space-time point (see, for instance section 13 of [37] for 
further details). 

We also remark that, when working with metrics having Lorentz measure, we can directly infer 
from the splitting of the infinitesimal length square into proper time and proper spatial length 
whether space-time is locally Lorentz or not, considering an ADM parameterization of the metric 
we obtain that dr = Ndt 2 must be time-like and dl 2 = —gij(dx l + N l dt)(dx^ + N^dt) must be 
space- like (N is the lapse function and N l are the shift functions). With respect to the specific 
metrics discussed here, we have that the SC metric (|3.2p describes local Lorentz space-time outside 
the event horizon r% > n.sc an d the FLRW metric (|2.28p describes local Lorentz space-time 
within the cosmological horizon r% < Ih- Also considering this splitting allows us to build an 
ansatz for matter in an expanding background following a very simple procedure. We note that 
for expanding coordinates r\ the deformation of the Minkowski metric due to spatial expansion 
which corresponds to the FLRW metric (I2.28P is given by an additive radial shift function 





a n r\ 



(4.1) 



FLRW Q~g2 
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such that the radial infinitesimal element is shifted to dr± — > dr\ — iVp LRW cdt, where Ijj = c/H 
is the time-dependent Hubble length (jl.5|) . To show it explicitly let us rewrite the FLRW met- 
ric (|2.28p factorizing the infinitesimal length square into proper-time and proper-length obtaining 



ds 2 = c 2 dt 2 _ ^ -°£cdt) -r\ (d9 2 + sin 2 9dy 2 ) . (4.2) 

The Lorentzian measure corresponding to the spherical coordinates measure \f—g = r 2 sin 9 is 
maintained whenever the corrections are encoded in an additive radial shift function. Hence, 
considering a deformation of the Schwarzschild metric (|3.2p by a generic shift function depen- 
dent on the radial coordinate and time N (r,t) such that the infinitesimal radial element is 
shifted to dr\ — ► dt\ — N^dt, ensures that both spherical symmetry and the metric measure are 
maintained. To describe expansion effects N 1 has to necessarily contain the expansion factor 
^flrw = r iA# (|4.ip which ensures that at spatial infinity we retrieve asymptotically the FLRW 
metric (|4,2p . As for the SC metric, it converges asymptotically to the Minkowski metric due 
to its dependence on the factor (1 — Use) = 1 ~~ 2GM/(c 2 n) which converges asymptotically 
to is unity, lim ri _^ 00 (l — Use) = 1- Furthermore, at the Schwarzschild radius this factor is null 
lim n _* ri sc (l — C^sc) = 0. Then, in addition, considering the shift function iV 1 to depend also 
on a multiplicative positive power of the factor (1 — Use) < 1 simultaneously ensures that near 
the massive object the expansion effects decrease being exactly null at the SC radius and, at 
spatial infinity, the FLRW shift function A"p LRW (|4.ip is recovered. In this way the value of the 
exponent of the factor (1 — Use) fine-tunes the intensity of the expansion effects felt near the 
mass. Specifically increasing the positive exponent will decrease the value of the shift function for 
relatively small radial coordinate n, still maintaining the convergence to the asymptotic FLRW 
metric for large values of the radial coordinate. It is also sensitive to note that both the isotropic 
McVittie metric (|3.34p and the anisotropic CSC metric (|3.38[) discussed in the previous section are 
obtained from the SC metric by considering a deformation by such a shift function, respectively 
with powers of the factor (1 — Use) of exponents 1/2 and 1. 

Given the previous discussion we are considering the deformation to the Schwarzschild metric, 
due to expanding background, to be given by a radial shift function of the form 



i 



rl _ an f 2GM~\ 2+2 
a c \ c 2 



such that we obtain the metric ansatz 



ds 2 = ^l_^^ c 2 dt 2_ r 2^2 +sin 2^ 




„.i X 2 M) 
an ( 2GM\ 2+5 \ 

dr\ 1 t, I cdt 

ac \ cr n 

Here we are considering a to be a (real) constant. This metric ansatz is asymptotically isotropic 
at spatial infinity, being only globaly isotropic for a = which corresponds to the McVittie met- 
ric (|3.34|) . As for the CSC metric (|3.38[) it corresponds to a = 1. By direct inspection this metric 
has the desired properties, either in the massless limit and at spatial infinity converges asymp- 
totically to the FLRW metric (|2.28p . in the static limit (d — * 0) coincides with the Schwarzschild 
metric (j3.2|) and it has Lorentzian measure y/—g = r 2 sin 9 coinciding, for all space-time, with the 
measure of both SC metric and FLRW metric. Hence this metric describes local Lorentz space- 
time in between the Schwarzschild horizon n = n.sc = 2GM/c 2 and the corrected cosmological 
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horizon corresponding to the solution of the equation IhIt\ = (1 — 2GM/(c 2 ri))( Q+1 )/ 2 . Assuming 
the weak field approximation (|3.7p and taking the solution for the cosmological horizon to first 
order in the gravitational field Use (|3.ip . we obtain that space time is locally Lorentz in the range 



n g 



2GM (« + l)GM 

o— i «H H o 



(4.5) 



There is one more property of this metric that has relevant physical implications. We note that 
by taking the constant parameter a to infinity we recover the Schwarzschild metric. As already 
pointed out this feature allows to fine-tune the expansion effects in local systems, outside the 
Schwarzschild radius the factor (1 — Use) increases monotonically with r\ up to spatial infinity 
being strictly less than unity (being null at the SC radius and unity at spatial infinity). Hence for 
relatively large positive values of a the expansion effects are highly suppressed near the massive 
object being still relevant for large values of the radial coordinate r±, at spatial infinity the metric 
still converges to the FLRW metric (I2.28f) . It is therefore physically intuitive to expect relatively 
high values for this parameter, expansion effects are for small scale astrophysical systems (for 
instance the solar system) negligible within the experimental accuracy and the SC metric (|3.2[) 
describe to a very high accuracy these systems. 

Next we analyze the singularities, as well as the asymptotic behavior of the scalar curvature 
R a and the curvature invariant lZ a at the Schwarzschild horizon n = ri.se an d at the center 
of mass 7*i = 0, computing the allowed ranges of the parameter a for which the SC horizon is 
singularity free, space-time is asymptotically flat at this horizon (such that the ansatz converges 
asymptotically to the SC metric) and the SC mass pole is maintained at the origin (such that the 
mass inside a shell of finite radius is finite). 



4.2 Singularities and Curvature at the Schwarzschild Radius: Lower a Bounds 

Depending on the value of the parameter a, the locally anisotropic metric (|4.4p has a space-time 
singularity at the Schwarzschild radius r\ = 2GM/c 2 , this is the case for a = and a = 1 
corresponding to the metrics (|3.34[) and (|3.38p already discussed in section [3l Although in the 
previous section we have assumed a positive exponent of the shift function, for completeness of 
our analysis, we are proceeding with the analysis for the full possible range of the parameter 
a g] — oo,+oo[. The Ricci scalar R a (scalar curvature) and the curvature invariant lZ a are 
given, respectively, in equations (1D.4P and (ID.7P of appendix [Dj By direct inspection of these 
expressions, we obtain the regularity and asymptotic leading expressions of these quantities at the 
Schwarzschild radius, r\ — ► ri se = 2GM/c 2 , as listed in tableland table EJ respectively, for the 
several distinct ranges of the parameter a. 

For a < 3 the Ricci scalar R (the scalar curvature) has a singularity at the Schwarzschild radius, 
while for a = 3 it is finite and for a > 3 it is null. It is explicitly checked from the curvature 
invariant 1Z that for a < 3 the Schwarzschild radius is singular. Based in these results we conclude 
that only for a > 3 space-time is free of singularities at the event horizon. Neither the McVittie 
metric (|3.34p nor the CSC metric (|3.38[) discussed in the previous section [3] obey this bound. 
A space-time singularity at the Schwarzschild radius is clearly a significant deviation from the 
Schwarzschild metric meaning that, near the massive object, the expansion effects are dominant 
with respect to the matter effects described by the original Schwarzschild metric. Such result is 
clearly in disagreement with the existing experimental evidence, the Schwarzschild metric (|3.2H 
has been widely applied in small scale astrophysical systems being in very close agreement with 
experimental data as well as with most tests of the General Relativity Newton law, in particular 
in the solar system and earth based experiments. Hence it is expected that near the massive 
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a R a regularity R a leading asymptotic expressions 
for ri -* ri.sc for n ~ r^gc 

a G ]+3, +00 [ finite ^qOlsc) = ^ScOlSc) = 



a = +3 finite R a (n.sc) = 3 ((A) 2 - JL 



a E [— 1,3[ divergent 
a S] — 00, — 1[ divergent 



a 3 



2GM\ 2-2 
ric 2 / 

2GM\ 



a-1 



ric 2 



Table 2: Regularity and asymptotic values of the Ricci scalar R a &D-4\ ) (the scalar curvature) for expanding 
coordinates near the Schwarzschild radius r% ~ 2GM/c 2 . 



objects we should obtain, at most, small corrections to the Schwarzschild metric, a divergence is 
clearly not a small correction. Also we expect that, to exist, any space-time singularity to lie at 
the origin (as it does for the SC metric). Given this discussion we conclude that the bound a > 3 
must be obeyed. 

This bound can also be justified at the level of the metric by noting that in the limit r\ — ► n.sc 
the asymptotic leading expressions for the metric components are 



500 



90r 



9r 



r\ c 2 J a 2 c 2 \ r\ c 2 J 



a c 



1 



an A 2GM\ 2 , A „s 

' 1 1 02" 2 (4.6) 



n c 2 / 



1 



1 _ 2GM Ql 



A necessary condition for the metric to converge to the SC metric at the Schwarzschild radius is to 
demand that, in this limit, all terms containing corrections due to expanding background vanish 
at least as fast as the terms corresponding to the original metric components of the Schwarzschild 
metric (13. 2ft . From the component <?o we obtain the lower bound a/2 — 1/2 > 1 which coincides 
with the above bound, a > 3. 

For all values of a, at spatial infinity (n ~^ +°°)> the curvature converges asymptotically to the 
FLRW curvature 

lim R a = -6 [ (—) + -A.) = i?FLRW = 6(? - 1)H 2 , (4.7) 
ri-+oo \\acj a& 

The expression for -Rflrw is given in equation (|B.2|) of appendix [Bl and we have replaced the time 
derivatives of the scale factor by the time-dependent Hubble rate H and deceleration parameter 
q. This result is expected, by construction the locally anisotropic metric (|4.4p converges asymp- 
totically to the FLRW metric (12.280 . The asymptotic value of the curvature at spatial infinity 
will be positive for q > 1, null for q = 1 and negative for q < 1. 
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a lZ a behaviour TZ a leading assymptotic expressions 
for r\ -> ri.sc for r x ~ r^gc 

a G ]+3, +oo[ finite ^ a (r LSC ) = ftsc (?"i.sc) = 



a = +3 finite 7^a(ji.sc) = 12 



c 2 
2GM 



1 



^(2^?) ((ac) ac 7 ) + ^ ( ( ac ) ac 



a G [+1,3[ divergent 
a G] — 00, l[/{0} divergent 
a = divergent 



2GM\ 2a " 4 



ric 2 



2GM\ 
r\c 2 J 



-1 



Table 3: Regularity and asymptotic leading expressions for the curvature invariant lZ a \D. 7| ] for expanding 
coordinates near the Schwarzschild radius n ~ 2GM/c 2 . 



Related to the previous discussions, there is one more relevant issue that we want to address, the 
sign of the curvature near the Schwarzschild radius and the value of its spatial derivative. This 
analysis allows to determine the asymptotic behavior of the curvature near massive bodies which 
is related to the gravitational interactions for small spatial scales (small astrophysical scales as 
planetary systems). We will only consider the case for a > 3 for which the SC radius is non- 
singular as we have shown. Close to the massive body, there are several distinct asymptotic 
behaviors which depend on the value of the deceleration parameter q. To identify these distinct 
asymptotic regimes let us consider the leading term for the spatial derivative of the curvature near 
the SC radius 



dr\ 



3{l + g )H 2 (GM) 2 fa _ 3\ / (6 - a) GM \ f _ 2GM \ £~i 



The regularity and asymptotic values for this expression at the event horizon, for a > 3, are listed 
in table HI From these asymptotic values we conclude that space-time is asymptotically flat near 
the Schwarzschild radius (-R a ( r i.sc) = ^a( r i.sc) = 0) for a > 5 and all values of q or for a > 3 
and q = — 1 (for this case the curvature is negative up to spatial infinity). 

With respect to the curvature behavior in the radial coordinate range r\ G]ri.sc> +°°[ we obtain 
that: for a > 5, depending on the value of q, we will have distinct behaviors for the curvature , 
for q < — 1 the curvature is strictly negative, for q > +1 it is strictly positive and for q g] — 1, +1[ 
it is positive near the SC radius r\ = ri.sci with growing r\ it grows to a maximum and then 
decreases becoming negative and converging to the FLRW value at spatial infinity; for a G]3,5] 
the curvature is null at the SC radius, being its derivative negative or null for q < —1 or q = — 1, 
respectively (for which cases the curvature is negative up to spatial infinity), or positive for q > 1 
(for which case the curvature is positive until spatial infinity); for the particular case of a = 3, at 
the SC radius ri.sc> the curvature is negative, null or positive for q < —1, q = — 1 and q > — 1, 
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a R' a regularity R' a asymptotic value 
for ri -> n.sc for r x -> r LS c 

a G ]+5, +oo[ finite ^4( r i.sc) = 

a = +5 finite R' a (n . sc ) = 5 {[f,f - ^ 

a G]3, 5[ divergent ^(ri.sc) = — oo for q > — 1 

finite -R'a(ri.sc) = for g = — 1 

divergent i?' Q (ri.sc) = +oo for g < — 1 

« = 3 finite J R' Q ( ri . sc ) = - 3 ((^) 2 + ^ 



Table 4: Regularity and asymptotic values for the spatial derivative of the curvature R' a for expanding 
coordinates near the Schwarzschild radius r\ ~ 2GM/c 2 . 



respectively, and its derivative is negative, null or positive for q < 1, q = 1 and q > 1, hence for 
the deceleration parameter range g g] — 1, 1[ the curvature will be positive near the SC radius and 
negative for large values of the radial coordinate (up to spatial infinity) . 

Specifically in the range of the deceleration parameter q g] — 1, 1[, which contains as a particular 
case the estimative for today's value of this parameter q$ = —0.589, and for a > 3 the curvature 
will always be positive near the massive object becoming negative with growing radial coordinate. 
The value of the radial coordinate for which the curvature changes sign is, approximately to first 
order in the gravitational field Use an d the parameter a, 

2GM 2a(2 - q) + 3(1 + q) „ . . , A , 

R a = => r x w — 2 6(1 -g) g €] - 1,+1[ , a > 3 . (4.9) 

This expression is derived considering an expansion of the curvature in the weak field approxima- 
tion being valid only for aU^c = 2aGM / '(c 2 n) < 1. We note that for relatively high values of 
the parameter a for which a > 1/C/sc it is required to compute numericaly the solution of the 
equation R a = by considering the exact expressions for the curvature. 

Resuming, from the analysis of the space-time singularities at the Schwarzschild radius carried in 
this section, we have obtained the following lower bounds for the parameter a 

a > 3 =>■ no space — time singularities at r\ = ri.se > 

(4.10) 

a > 5 space — time asymptotically flat at r\ = ri se • 

Some examples for several values of the parameter a and the deceleration parameter evaluated at 
the reference time to = 0, go = g(0) are presented in figure [2j 



43 




Figure 2: Curvature R ljD-4\ ) as a function ofr± for several values of the deceleration parameter qo evaluated 
at the reference time to and parameter a (see tabled and table^for comparison) : 

(a) example for a — 6 and several values of qo = — 1.5, — 1, —0.5, 0, +0.5, +1, +1.5, at the SC radius 
R = R' = for all plotted lines and, at spatial infinity, all converge to qo — 1; 

(b) an example for qo > 1 for a = 3, 4, 5, 10; at the SC radius only for a = 3, R > 0, R' > and only 
for a > 5, R = R' = 0, the plotted lines converge at spatial infinity to qo — 1 = 0.5 represented by a dotted 
horizontal line and the curvature is positive for r% G]ri.sC) +°°[/ 

(c) an example for the particular case q a = +1, at the SC radius only for a — 3, R > 0, R 1 = and only 
for a > 5, R = R' = 0, at spatial infinity the plotted lines converge to qo — 1 = and the curvature is 
positive for r 1 £]rx.sc> +°°[/ 

(d) estimative for today's value of the deceleration parameter qo = —0.589 for a = 3, 4, 5, 10 which is 
an example of q £] — 1, +1[, at the SC radius only for a — 3, R > 0, R' = and only for a > 5, R = R' = 0, 
at spatial infinity the plotted lines converge to qo — 1 = —1.589, the curvature is positive near the massive 
body and negative at spatial infinity, changing sign at a finite value of r\ 

(e) the particular case of qo = —1 for a = 3, 4, 5, 10, at the SC radius only for a — 3, R = 0, R' < and 
for all other values of a > 3, R = R' = 0, at spatial infinity the plotted lines converge to qo — 1 = —2, the 
curvature is negative in the range r\ €]ri.sc, +oo[; 

(f ) an example for qo < —1 case for a = 3, 4, 5, 10, at the SC radius only for a — 3, R < 0, R! < and 
only for a > 5, R = R' = 0, at spatial infinity the plotted lines converge to q a — 1 = —2.5, the curvature is 
always negative in the range r% €]rx.sc> +oo[. 
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These results are physically appealing, in particular for today's value of the deceleration parameter 
qo = —0.589 and for a > 5, near the massive bodies space-time is asymptotically flat. Hence 
the usual gravitation laws for flat backgrounds obtained either from General Relativity or the 
Newtonian limit approximation should be asymptotically (or at least approximately retrieved 
to a very good accuracy). We recall that these laws are experimentally well established being 
tested to a very high precision for small scale spatial scales systems such as the solar system. 
Far away from the massive objects, the expansion effects dominate and the curvature converges 
asymptotically to the FLRW curvature which is also experimentally well established result from 
large scale observations (for instance the Hubble law). In between these two asymptotic regions 
there is a transition region for which the gravitational interactions will be modified such that an 
interpolation between the two asymptotic limiting cases is obtained. The parameter a fine-tunes 
the transition between these distinct regimes. 

We will study these three regimes in detail later on, for now let us complete our analysis of the 
space-time singularities by analyzing the curvature R a and curvature invariant lZ a at the center 
of mass, the origin of the coordinate frame r\ = 0. 

4.3 Singularities at the Center of Mass: Mass Divergences and Upper a 
Bounds 

Having already discussed the space-time singularities at the Schwarzschild radius, it is still nec- 
essary to analyze the singularities at the center of mass r\ = 0. Independently of the value of 
the parameter a, the mass pole divergence for the Schwarzschild metric (|3.3p . is also present for 
the locally anisotropic metric (14, 4h . However, depending on the value of a, there will exist other 
contributions for the singularity at the origin. We list the regularity and the asymptotic leading 
terms at the origin for the Ricci scalar R a (curvature) and the curvature invariant TZ a in table [5] 
and table [H 

a R a regularity R a leading assymptotic expression 

for n — > for r\ ~ 



a G ]0,+oo[ divergent ~ — 

rf 

GM ,i\ 

a £]- oo,0 mass pole R a (0) = R S c(0) = -8vr^-(5 (3) (0) 



Table 5: Regularity and asymptotic leading expressions for the Ricci scalar R a &D-4\ ) (curvature) for 
expanding coordinates near the center of mass n *~ 0. 



At the center of mass n = 0, for a < 0, the Ricci scalar R a coincides with the Schwarzschild 
curvature (]3.3|) . while for a > it is divergent by the power of the radial coordinate r^ a . As 
for the curvature invariant lZ a , for a < —9, it asymptotically coincides with the Schwarzschild 
curvature invariant T^sc ~ r \ 6 (|3.3p . while for a G [—9, +3] the leading divergent term in the 
neighbourhood of r% = is still the SC curvature invariant ~ rj~ 6 , however the divergence at the 
center of mass has several other contributions of lower l/r\ powers. For a > 3 the curvature 
invariant leading divergent term is rj~ 2a , hence a more severe divergence than the SC mass pole 
divergence such that for this range of the parameter the total mass at the origin diverges. 

This discussion is not complete without explicitly computing the mass-energy density and the 
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a TZ Q behavior TZ a leading asymptotic expression 

for n —* for n ~ 



a S ]+3,+oo[ divergent ~ 



1 



a G [—9, +3] divergent ~ — 



a S ]— oo, — 9[ divergent Tl a (r\ ~ 0) = 7£sc( r i ~ 0) ~ — 



Table 6: Regularity and asymptotic leading expressions of the curvature invariant lZ a ljD-4\ ) for non- 
expanding coordinates near the center of mass r\ ~ 0. 



mass content of the space-time described by the locally anisotropic metric (|4.4p . In the commoving 
frame of the background cosmological fluid the mass-energy density is, generally, p = Tqq/c 2 , given 
in terms of the component Tqq of the stress-energy tensor in the Lorentz frame (Cartan frame). 
We give the details of this computation in section 14.51 f° r now let us quote the expression for the 
total mass-energy density 



Ptot(a)( r l) = PSC+Pa , 



Psc 



(4.11) 



where we are explicitly considering the Schwarzschild mass pole contribution (|3.5p to the mass- 
energy density such that this expression has both a contribution from the standard point like mass, 
the SC mass pole, a localized distribution at r\ = 0, and a contribution from the background, hence 
an extended mass-energy distribution p a . We recall that the background mass-energy density in 
the absence of the local mass M is pu = 3(d / a) j (8tvG) as given in equation (|2.57p and note that 
for all values of a, p a asymptotically converges to this mass-energy density at spatial infinity, 
lmv^+oo p a = ph- Hence the deformation of the background mass-energy density is encoded 
in p a , such that to evaluate the deformation (or correction) pdef to the background mass-energy 
density due the presence of the mass M is necessary to subtract the background mass-energy 
density pu 

Pdef(o) = PSC + Pa- PH ■ (4.12) 

We note that in this discussion we are interpreting the locally anisotropic metric (|4,4p as a defor- 
mation (or correction) to the cosmological expanding background given by the FLRW metric (12.281) 
due to the presence of a local central mass M . When deriving the ansatz for this metric we have 
considered the opposite interpretation, that it is a deformation (or correction) to the Schwarzschild 
metric (]3.2p . Both these interpretations are correct being compatible with each other, we recall 
that the locally anisotropic metric (|4,4p is interpolating between both the FLRW and the SC 
metrics, hence it is necessarily a deformation (or correction) to both of them. 

Generally the mass (or equivalently, the energy) within a shell of fixed radius is obtained by 
integrating the mass-energy density over the respective volume. Specifically the total space-time 
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mass and the mass of the deformation due to the presence of a point-like mass M within a shell 
of fixed radius r\ , are given by 



M tot{a) (n) = M + AM tot{a) ( n ) , 



Mtef(a)(ri) = M + Atf t ot(a)(n)-^(n) 



(4.13) 



AM 



tot (a) 



(n) 



M a (n)-M a (0) , 



where M is the usual Schwarzschild gravitational mass contribution from the density p$c, AM tot ( a ) 
is the contribution due to the expanding background when the mass M is present from the density 
p a and Mff is the background mass in the absence of the mass M corresponding to the density 
Ph- Specifically we have that M a (r\) is given by the indefinite volume integral of p a 



2- 



M Q (n) -- 

and Mh by the volume integral of pn 
M„{rx) 



1 10 / dip I dri^J—g p 



x 2 3 



2G 



c 2 r\ ) 



dO I dip I dr[* 
-71 Jo Jo 



~9PH = 7^7=7 1 

2G \ ac 



(4.14) 



(4.15) 



The measure is given, as usual, by y/—g = r\ sin#. 

At the origin M Q (0) has three distinct behaviors which depend on the value of a being greater, 
equal or smaller than 3, 



( (-l) Q oo 



M«(0) 



a > 3 , 



a c 



a 



cf 



3 , 



(4.16) 



I 



a < 3 , 



therefore the total mass for a shell at fixed radial coordinate ri, including both the SC mass 
pole contribution M and the contribution of the extended deformation due to the expanding 
background AM tot(a) (ri) is 



^tot(a)( r l) 



( ±oo 



M + 



M + 



M ( GM\ 

72 



2G 



2GM\ a 
1 - —3 — ) 




a > 3 



a 



a < 3 



(4.17) 



Hence, for a > 3, the total mass inside a shell of finite radius r\ is divergent while for a < 3 it 
is finite. This result is consistent with the analysis of the curvature invariant singularities given 
in table EJ for a < 3 the leading divergence of lZ a is a power of the radial coordinate coinciding 
with the inverse volume squared 



~ r 



1 > 



while for a > 3 it has a more severe divergence by a 
power of ~ a . We also note that for a < 3 the mass pole value at the origin is maintained, 
lim 



ri ^ M tot(a < 3) (ri) = M. 
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We recall that at spatial infinity the density p a converges asymptotically to the background density 
Ph however, depending on the value of the parameter a, the respective mass expressions do not 
match this behavior. This is due to the value of the indefinite integral M a at the origin. Hence 
taking the limit of large radius for the volume integration shells we obtain 



^tot(a)Ol ~ +00) = 



For a > 3 the total mass is divergent, for a = 3 it corresponds to the background mass Mg plus 
the SC gravitational mass M of the massive body plus a small correction proportional to H 2 and 
for a < 3 it is exactly given by the background mass Mg plus the SC gravitational mass M. 

Resuming the results obtained in this section, from the analysis of the singularities at the center 
of mass, we have set the following bounds on the parameter a 

a > 3 =>■ space — time singularities at n = diverge with the power ~ r^ 2a , 
total mass inside shell of radius r\ is divergent , 

a = 3 =>■ space — time singularities at T\ = diverge with the power ~ rj~ 6 , 
total mass inside shell of radius r\ is finite 

(4.19) 

being above M + Mn for r\ ~ +00 , 

a < 3 => space — time singularities at n = diverge with the power ~ rj~ 6 , 
total mass inside shell of radius r\ is finite 

being M + Mh for n ~ +00 . 

Then we conclude that for the locally anisotropic metric (|4.4p only for a < 3 the origin is free of 
essential singularities such that the total mass within a shell of finite radius is finite. Hence, in 
order to avoid mass divergences, this bound should be imposed. In addition we note that only for 
a < 3 at spatial infinity we recover that the total mass is the sum of the background mass Mjj with 
the SC gravitational mass M while for a = 3 there is a small positive correction, H 2 G 2 Af 3 /c 3 . 
If we require to exactly maintain the relation between total mass and gravitational mass in the 
universe the stronger bound a < 3 should be considered, however we note that this correction is 
for most purposes negligible. These bounds are the opposite of the ones obtained previously when 
analysing the space-time singularities at the Schwarzschild radius, r\ = n.sc fl4.10j) . In the next 
section we discuss and deal with this (in) compatibility. 

As a final remark let us note that within the Schwarzschild radius the metric (|4.4p becomes, 
generally, complex, except for odd integer values of the parameter a for which a signature flip 
is obtained. The locally anisotropic metric (|4.4p does not describe a local Lorentz space-time 
inside the event horizon (|4.5p . this characteristic is inherited from the Schwarzschild metric and 
it is attributed to the coordinate choice, the Schwarzschild coordinates only describe the physical 



M + M H + M[^j '{^f) 2 , a = 3, 



(4.18) 



M + Mu 



a. < 3 
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space-time outside of the event-horizon. In order to fully describe the horizon inner region and 
properly compute geodesic paths it is necessary to consider other coordinate choice such as Novikov 
coordinates or Kruskal-Szekeres coordinates [33j. We are not further discussing this topic here, for 
further details see for instance section 31.4 of [37J. We also note that the above expressions for the 
total mass M tot r a \ (|4.13p . M a (|4.14p and AM tot ( a \ (|4.1T[) are well defined as long as we consider 
the integration shell of fixed radius r\ outside the SC event horizon, i.e. t\ > ri.sc = 2GM/(c 2 r\), 
and that for a > 3 the limit lim ri _»o -^tot(a) = -^tot(a)(0) dUBD is also well defined being a real 
value which is enough for mass computation purposes and evaluation of singularities. 

4.4 The Ansatz II: Removing Essential Singularities at the Center of Mass 

The results of the two previous sections raise a problem that must be solved. In section 14.21 
we have concluded that only for a > 3, at the Schwarzschild radius r% = rsc, space-time is 
singularity free (|4,10p while in section 14.31 we concluded, from the analysis of the singularities 
at the origin, that only for a < 3 the total mass inside a shell of finite radius is finite (14. 19|) . 
Compatibility between these two bounds leave us with the only possible value for this parameter 
to be a = 3. This value is theoretically consistent however, based in the argument that, when 
compared to experimental data, the Schwarzschild metric describes astrophysical gravitational 
systems to a very high precision, we expect that the deformation of this metric, as encoded in 
the ansatz (|4.4B . has negligible corrections for short spatial scales. Namely, in agreement with the 
General Relativity and the Newtonian gravitational laws, we may expect that near massive bodies, 
space-time is asymptotically flat. This requirement corresponds to the lower bound a > 5 (14.10P 
which are clearly incompatible with the upper bound a < 3 (|4. 19j) . 

We will next solve this incompatibility by modifying our metric ansatz ()4.4p . So far we have 
considered the parameter a to be a constant, generally we may assume it to be a space-time 
dependent function. We recall that in the physical system being addressed, matter in an expand- 
ing background, there are only two dimensionless quantities, the FLRW expansion shift function 
^flrw = r i/lH (|4.ip and the SC gravitational potential Use = 2GM/(c 2 ri) ([33]) . As discussed 
in the previous section 14.21 and section 14.31 we further note that the parameter a fine-tunes the 
expansion effects relatively close to the massive objects, hence its value is mostly relevant for 
relatively small radial scales (specifically for r\ <C Ih) for which the matter effects are dominant. 
As for large radial scales (specifically for n ~ Ih) the expansion effects are dominant, the gravita- 
tional potential Use decreases with the radial coordinate (lim ri ^ +00 Use = 0) and the expansion 
shift function increases with the radial coordinate (-ZVp LRW ~ r±) such that the factor (1 — Usc) a 
converges asymptotically to unity at spatial infinity independently of the value of a (as long as a 
is finite). It is therefore consistent to expect a space-time dependent parameter a to depend only 
on the dimensionless gravitational field. Hence let us assume that its dependent on Use is linear, 
specificaly 

2GM 

a(ri) = a + ai Usc(n) = «o + «i —5 • (4.20) 

cr n 

This exponent maintains spherical symmetry of the metric (|4.4p converging asymptotically to the 
constant coefficient ao at spatial infinity. For relatively small values of the radial coordinate r±, 
near the massive body, the effects of its dependence on the gravitational potential Use become 
relevant. Hence the coefficient a\ fine-tunes the exponent dependence on the gravitational field 
near the massive object. 

Next we discuss the effects of this new term in the exponent. The connections, curvature and cur- 
vature invariant for a space-dependent exponent a are given in equations (jD.lOp . (ID. 121) and (ID. 161) 
of Appendix [Dj With respect to the singularities at the Schwarzschild radius we note that, for 
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this value of the radial coordinate, the exponent is a[r\ = rj.sc) = a o + a i such that the same 
bounds on a as expressed in (|4.1U|) are valid for the exponent evaluated at the event horizon, 

o:(n.sc) = ao + «i > 3 =^ no space — time singularities at n = ri se j 

(4.21) 

a (ri.se) = Oio + ai > 5 space — time is asymptotically flat at n = ri se • 

The proof of these bounds is straight forwardly obtained by noting that lim,^-^ sc (l— Usc) p l°g(l — 
Use) = f° r P > such that all terms in the curvature and curvature invariance containing the 
derivative of the exponent a' vanish at the event horizon. 

As for the singularities at the center of mass, we retrieve a constant exponent for a\ = corre- 
sponding to the case analyzed in the previous section with a = ao, for ot\ > the divergences at 
the origin, independently of the value of ao, are more severe than for the case of constant a such 
that we always obtain a divergence by a positive infinite power of 1/ri 

1 \ + ~ 



ft(ao,«i>o)(ri~0)~ I -j . (4.22) 

Finally, for a\ < 0, the exponent a(r) diverges at the center of mass, the origin of the coordinate 
frame, to lim ri _».o a(rj) = — oo. Hence by direct inspection of the curvature invariant lZ a given 
in equation (|D.16j) of appendix [D] we conclude that all the terms which dependent on the scale 
factor a vanish at r\ = 0. This result is established by noting that, for all p, q and for ct\ < 0, 
lim ri ^ r p (l - t/sc)~ KI/ri = and lim ri _ r p (l - £/sc)~ M/ri log(l - U SC ) g = such that, 
for all values of ao, the only contribution to the space-time singularity is exactly given by the 
Schwarzschild curvature invariant (13.31) 



^(ao,ai<0)(n ~ 0) = ft SC = ^M) 2 (423) 

c r x 

The total mass-energy density for a radial coordinate dependent exponent a(ri) is given, similarly 
to (|4.1ip . by Ptot(a(ri)) = Psc + Pa(n)j where the extended distribution p a ( n ) is 



Pa(ri) 



c 2 6o 8ttG 65 




2GM\ 

c 2 n / 



a-l 



(4.24) 



2GM\ a , (_ 2GM\ 
c 2 ri 



In this expression and ^g? are the Einstein tensor and stress-energy tensor for the 

locally anisotropic metric (|4.4p in the Cartan- frame (excluding the SC mass pole contribution). 
The details of the computation of these tensors are discussed in the next section. 

The total mass inside a shell of constant radius ri is obtained, as usual, by considering the 
respective volume integral. We follow the same definitions for the several mass quantities of the 
previous section l4~3l Similarly to the definition of the quantity M (|4.14p . for the specific exponent 
a{r\) = a>o + 2GM/(c 2 r\) (|4.20p and considering the indefinite integral of the above density (I4.24|) 
we obtain that 



d\ 2 r? r 2GM\ ao+ai 



2GM 
" r l 
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At the origin, this quantity is divergent for ol\ > and null for a\ < 

e tS oo , a\ > , 

(4.26) 

, «i < , 

while for a\ = we retrieve the case of constant exponent a = ao analyzed in the previous 
section 14.31 In the above expression e tS is a generic complex phase. Hence consistently with the 
singularities analysis of the curvature invariant 1Z( a0tai ) as expressed in equations f|4.22j) and (|4.23p . 
independently of the value of ao, the total mass is divergent for a\ > and is finite for a\ < 0. 
Furthermore, for a± < and all values of ao, at the origin the singularity exactly matches the SC 
mass pole singularity such that the mass coincides with the SC mass pole value M. As for the 
total mass within a shell of large radius (r% ~ +oo), for a\ < and all values of ao, is given by 
the sum of the SC gravitational mass M with the background mass (in the absence of the local 
mass M) M H J£IB]> . M tot(a0)ai<0 ) (n ~ +oo) = M + M H . 

Resuming the results obtained in this section from the analysis of the singularities at the center of 
mass, we have concluded that when considering a radial coordinate dependent exponent a(ri) = 
ao + ai Use (14.201) . independently of the value of the parameter ao, we obtain the following bounds 
on the parameter a\ 

a\ > =>■ space — time singularity at r\ = diverges by ~ rf°° , 
total mass inside shell of radius r\ is divergent , 

(4.27) 

ai < =^ space — time singularity at r\ = coincides with SC singularity 

(^a ,oi>o = ~ ^r 6 ) > 

total mass inside shell of radius r\ is finite . 

Then we have accomplished our main objective for this section, by considering a radial coordinate 
dependent exponent a{r\) = ao + a\ Use (14.20P with a\ < we have removed the essential singu- 
larities at the origin (except for the SC mass pole) maintaining the main properties of the locally 
anisotropic metric (|4.4I) . Furthermore for this exponent choice it strictly converges asymptotically 
at the origin to the Schwarzschild metric. 

Also we note that assuming the above bound a\ < f|4.27j) and either of the bounds (|4.2ip we 
obtain the bounds ao > |ai| + 3 and ao > |ai| +5 such that the following inequality applies 

a > |ai| . (4.28) 

As a consequence of this relation the corrections to the several quantities outside the SC event 
horizon (r\ > ri.se = 2GM/cP) due to the coefficient a± are smaller than the ones due to the 
coefficient ao, ao > ai Usc( r i > r sc) such that the results derived in section [4721 with respect to 
the curvature behavior close to the Schwarzschild radius are qualitatively maintained (see figure[2]). 
In particular the bounds for a at the event horizon (|4.10p are valid for the exponent value evaluate 
at the SC radius a(ri.sc) = «o — |cti | such that for ao > 5 + |ai| space-time is asymptotically flat 
at the event horizon. 

Next we compute the stress-energy tensor in the Lorentz frame (Cartan frame) and both the 
mass-energy density and pressures of the background fluid for the locally anisotropic metric (|4.4p 
analyzing their properties. 



^(«o,«i)(0) = < 



51 



4.5 The Stress-Energy Tensor with Anisotropic Pressures and 
Positive Definite Mass-Energy Density 



In this section we analyze the stress-energy tensor for the locally anisotropic metric (|4.4|) with a 
radial coordinate dependent exponent a{r\) = «o + «i t^sc (|4.20[) with strictly negative a\ < 0. 
In particular we give the details of the computation of the mass-energy density P( a0iCri ) and 
the anisotropic pressures p r (a ,ai)} P6(a ,ai) an d Ptp(ao,ai)i showing that these quantities have the 
correct asymptotic leading expressions at the origin and spatial infinity consistently with the 
metric ansatz construction assumptions and interpretation as an interpolation (or deformation) 
between the cosmological FLRW metric (|2.28p and the Schwarzschild metric (13, 2j) . 

The stress-energy tensor is related, as usual, to the Einstein tensor by the Einstein equations 

T )±V 



TjS°' ai) = c 4 /(8vrG) G%°' ai) (see equation (TO|) in appendix|A|). The Einstein tensor G { " {n)) 



Gffl+AG^ (excluding the SC mass pole contribution) for a generic exponent a(r\) dependent on 
the radial coordinate is given in equation (|D.14p of appendix [D] with G^} given in equation (|D.6P 
and AG^ given in equation (ID.15D of the same appendix. Here we use the index notation '(a(ri))' 
when referring to a generic exponent dependent on the radial coordinate n and '(ao; a i)' when 
referring to the specific exponent a{r\) = ao + a\ Use f|4.20|> . Either by direct inspection of the 
Einstein tensor G^0 \ or by noting that the metric (|4.4p is spherical symmetric we conclude that 
this symmetry is maintained both by the mass-energy density P( ao>ai ) and anisotropic pressures 
Pr(oo,ai)> Pe(a ,a!) an d P^(a ,a 1 )- This statement is translated as that the pressures along the 
angular directions are identical for all space-time Pe( ao ,a 1 ) = Pip(a ,ai) while the pressure along 
the radial direction may generally be distinct p r ( a0tai )- We also note that /9( a0)Q;i ), p r (aa,ai)> 
PB(a ,ai) an( i P<p(a ,ai) are scalar quantities, here the pressures do not correspond to a vector along 
the spatial directions, instead generally we have three distinct scalars (only two when spherical 
symmetry is maintained), hence frame independent. 

In the following we will work in the Cartan frame (Lorentz frame) corresponding to the local flat 
frame with the Minkowski metric i]^ for which the stress-energy tensor is diagonal being explicitly 
stress and shear free (T^f '" 1 ^ = 0). Here hatted indexes represent the flat coordinates in the 
Cartan frame. Assuming a background perfect fluid with commoving velocity Un = (c, 0, 0, 0) the 
stress-energy tensor in the Cartan frame is 

T («o,«i) _ 2 T (a ,ai) r («o,ai) T (ao,ai) 2g x 

00 — c P{a , ai ) > 1 n — Pr(a , ai ) , 1 ^ — ^33 P0(a o ,<xi) ' 

where as already discussed, due to spherical symmetry, T^°' ai ^ = Tg?°' ai . To explicitly compute 
T^°' ai ^ let us consider a Cartan tetrad = e^dx^ such that the coordinate metric is related 
to the Minkowski metric by g^ u = e^e^rfru). Then, for the locally anisotropic metric (|4.4I) . we 
obtain the following non-null components of the tetrad 



o 2GM » an ( 2GM\ 2 
e u n = 1/I it , e n = 1 1 



= V 1 2 ' e = 1 2" 



(4.30) 



e 1 ! = — , , e 2 9 = n , e 3 o = r\ sin / 

2GM 



The non-null inverse tetrad components e/ 4 are straight forwardly obtained by raising and low- 
ering the flat indexes '/t' with the Minkowski metric r)^ and the coordinate indexes '//' with the 



52 



coordinate metric g^ v 



1 j _ an/, 2GM\ 2 



e ,° = 

" /-^ 2GAf a C \ 7"! C 2 



(4.31) 



1 _ I 2GM 2 1 - 



i V ri c 2 ' 2 ri '3 ri s [ n ( 



Hence the stress-energy tensor in the Cartan frame is obtained by considering the contraction of 
the coordinate stress-tensor with the inverse tetr; 
equation (|A.8j) we obtain the following relations 



the coordinate stress-tensor with the inverse tetrad elements T^ u = e^e u v T^ v . From the Einstein 



(a(ri)) _ C 4 



oo 8vrG 



(p °p (7 (a(ri)) 4- 2p °p 1 a {a{Tl)) +e 1 e 1 n {a{ri)) \ 

e o oo + ze 6 6 oi + e 6 6 n J 



4 

T (a(n)) _ _2 („ 0. ^Mn)) , Li ^(n)) 

6i ~ 8vrG V 1 01 +e 6 e i u n 

T (a(n)) _ A p ip 1 f7 (a(ri)) (4.32) 
ii 8ttG 1 1 11 

4 

J 22 " 8vrG e 2 e 2 G 22 , 

4 

As already mentioned, the cross component of the stress-energy tensor in the Cartan frame are 
null Tgj = such that in the local flat Minkowski frame the stress-energy tensor is stress and shear 
free. The non-null Tq\ is due to the (global) coordinate system choice and does not represents a 
measurable physical stress nor shear. 

In the following, for compactness of the expressions and direct numerical evaluation, we will 
consider the dimensionless radial coordinate rescaled by the Schwarzschild radius 

n 2GM 

h = — ^- = — ^ n . 4.33 

ri.sc c z 

Hence from the expressions (|4.32p relating the stress-energy tensor Tnp in the Cartan frame with 
the Einstein tensor G^ (ID.14D in the coordinate frame and from the definition of T^y for a 
(anisotropic) perfect fluid in the Cartan frame (|4.29p we obtain, for the particular radial coordinate 
dependent exponent a(r\) = Qo + ai Use ()4.20p the following solutions for the mass-energy density 
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P(ao,ai) and the pressures p r (a , ai ) and Pe{a , ai ) 



H 2 / l\ a0 " 1+? T/ o «o-3 ai 



log 1 



n n , -> £0_I-l-iil. / , s a 1 

2 H 2 / l \ 2 2^2f=! / /q-q «l\ / l\ 2 2 




PS(o!0iO!l) £V(ao>cn) 



c 2 H 2 A ]_\ a »- 2 +^ 



Mg\ 1 -TJ X (4-34) 



fo , 2 («o ~ 3) «o(«o -5) +4«i + 6 (2a - 5)ai af 
X V fx + 2ff + 2ff + 2ff 

c 2 R 2 a x ( ±V°~ 1+ ^V / "A, / 1 
+tt^T -5- 1 - - 1 + «o + log 1 - - 



8tt G r 2 \ T\) \ \ r\ J V ri 

c 2 # 2 ai / 1 \ Q0+ ^ A n ! ( 1 \A , / I 



c 2 # 2 ai(l + g) / lX^-l+^T / 1 



8vrG 2fi V 1 n/ ri 

Here we have replaced the derivatives of the scale factor by the Hubble rate H = a/a (jl.ip and 
deceleration factor g = —d/(aH 2 ) (II, 2p definitions. We note that the expression for the mass- 
energy density corresponds to the equation for generic exponent a(r\) (|4.24j) already considered 
in the previous section evaluated for the particular case of the first order exponent (|4.20p . 

We also recall that, similarly to the analysis carried in the previous section [431 (equation (|4.11|) ). 
the total mass-energy density p to t = psc + P(a ,ai) lS gi ven by the sum of both the Schwarzschild 
mass pole contribution psc A3.5H and the mass-energy density contribution P( aQ)Ql ) (j4.34jl due to 
the expanding background deformation. As for the total pressures correspond to the anisotropic 
pressures just computed, p r tot = P r (a ,ai) anci Petot = P9(a ,ai)> the Schwarzschild pressure is null 
everywhere (including the origin). With respect to the asymptotic limits of these quantities we 
obtain at the origin that, for a\ < 0, both the total mass-energy density and pressures consistently 



54 



coincide with the respective Schwarzschild quantities 



Ptot(O) 
Pr.tot(O) 

Pe.tot(O) 



lim n {PSC + P(ao,ai)) = ^SC = M5 (3) (ri) , 



hm Pe(a , ai ) = 



(4.35) 



Asymptotically, at spatial infinity or at least for large enough T\ such that the gravitational 
potential is negligible f7sc( r i ~ +°o) ~ 0, independently of the coefficients ao and ati, we retrieve 
the cosmological mass-energy density pn and pressure pu A2.57B corresponding to the FLRW 
metric ([238]) 

3H 2 

Aot(n ~ +oo) sa p# 



Pr(ao,ai) 

(ri 



P6»(a ,ai)( r l 



-OO 



-OO 



Ph 



PH 



8vrG ' 

(-l + 2g)c 2 ff 2 
8itG 

(-l + 2q)c 2 H 2 



(4.36) 



8^G 

At the Schwarzschild radius r\ = rj.sC) f° r cti < and ao + «i > 3, the values of these quantities 
are 

/Wn.sc) = , 



Pr(a , ai )(n.Sc) 



)(n.sc) 



(4.37) 



+oo 




a = 3 
a > 3 



As for the asymptotic values of the derivatives of the mass-energy density and pressures we obtain 
that, both at the origin and at spatial infinity, are null coinciding, respectively, with the behavior 
of the SC metric (J32D and the FRLW metric ^M, p'toM = P'r {ao , ai )(°) = P*(oo,ai)(°) = and 
p[ ot (+oo) = p' r ( aQ Q . 1 )(+°°) = P'e(a Oi)^ 00 ) = At the event horizon these quantities are 



Ptot(n.sc) 



Pr( Q0 , Ql )( r l-Sc) 



: 

( c 2 (l + q)H 2 



4vrG 



ckq — |cifi I > 3 



«0 — \ a l\ 



, ao — | ai| > 3 



(4.38) 



P'e(ao,ai)( ri 



SCj 



+oo 

hc 2 {l + q)H 2 



16vrG 







ao — |ai| G [3, 5[ 

ao — | ai| = 5 , 
ao — | ai| > 5 , 



These limits and respective bounds on the coefficient sums ao — |ai| are consistent with the space- 
time being asymptotically flat at the event horizon for ao — |ai| > 5 as expressed in equation (|4.21l) 
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based in the analysis of the asymptotic behavior of the curvature R. Here the same conclusion 
is obtained from the asymptotic behavior of the mass-energy density and pressures, all these 
quantities and its derivatives are null at the event horizon for a® — \ a\\ > 5 matching the respective 
quantities for flat space-time. In figure [3] are graphically represented the mass-energy density and 
pressures for several values of a\ and oq — \ai\. 




L , 1 1 , 1 1 , 1 1 1 1 1 , 1 1 , 1 1 , U -0.005 Lj ' ' ' ' 1 ' ' ' ' 1 ' ' ' ' 1 ' ' ' ' □ 

1.0 1.5 2.0 2.5 3.0 1.0 1.5 2.0 2.5 3.0 



Figure 3: Reseated mass-energy density P( Qo , Ql ) = (8irG/H 2 ) p( ag iBtl \ and pressures p r ( ao ,a 1 ) = 
{8itG/(c 2 H 2 ))p r(a0:ai) , Pe( ao , ai ) = (8ttG/(c 2 H 2 )) p 6 { aa , ai ) M-3$ as a function of n = n/n.sc U.33\) . 
The assymptotic values for large f\ are the respective quantities for the FRLW metric pn = (SttG / H 2 ) pu = 
3 and pu = (SttG/ (c 2 H 2 ))ph = (2go — 1) \2.5 r fy . Are given the examples for: 

(a) ao — \ai\ = 3 with ao — 4 and ct\ = —1, the pressure pe( aQ ai ) and its derivative P'g^ ao ai ) diverge at 
fx = 1; 

(b) ao — |ai| = 4.5 with ao = 5.5 and ct\ — —1, the pressure derivative p' s ^ ag ai \ diverges at r% = I; 

(c) ao — |ai| = 6 with ao = 7 and ct\ — —1, both mass-energy density, the pressures and their derivatives 
are null at f ± = 1 ; 

(d) ao — |ai| = 3 with ao = 30 and ai = —27, characteristics similar to the case of figure (a), however 
outside the SC event horizon fi > 1 exists a region of negative mass-energy density shown in the blow up. 
For ao = 30 the mass-energy density is positive definite for ct\ > 3 — 30A Q (30) ~ —23.805 

In addition, to preserve causality, it is required to demand the energy density to be positive 
outside the SC event horizon r\ > n sc- This condition is not verified for all negative values of 
the parameters a\ < and ao > 3 + \a%\. By direct inspection of the expression for /?( a0jai ) given 
in equation (|4.34|) we note that the sign of the full expression for f\ > 1 is set by the sign of the last 
multiplicative factor within brackets. The zeros of this expression cannot be found analytically, 
however a direct numerical inspection can be carried to define the allowed range for the coefficient 
a\ as a function of ao, simultaneously considering, at the SC event horizon, either the criteria of 
the absence of space-time singularities or space-time being asymptotically flat (|4.2ip . Hence it 
is further necessary to constraint the coefficient ot\ to be either in the range ol\ £ [3 — ao,0[ or 
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a\ S [5 — ao,0[, respectively. We note that here we are considering the bounds (|4.2ip keeping ao 
as generic as possible and constraining a\ such that for ao > 3 and ao > 5 we obtain, respectively, 
the following allowed ranges for the coefficient a± 



ao G ]3, +oo[ 



space — time at SC event horizon is singularity free 



a i 



G -4 {3) («o) 
a G ]5,+oo] 



(4.39) 



space — time at SC event horizon is assymptotically flat 



a 



i e 



a ) 



The sets *4 (3) («o) and ^ 5 )(ao) are negative ranges depending on the specific value of ao- Specif- 
ically are 



a J 



a J 



^S 3) («o) 

I 4 3) («o) 



4 5) («o) 



4 5) (a ) 



]3-a ,0[ , a G ]3, 28.347] , 

]3 - A Q (a ) x a ,0[ , a G ]28.347, +oo[ , 

]5 -a ,0[ , a G ]5, 45.226] , 

]5 - A Q (a ) x a ,0[ , a G J45.226, +oo[ , 



(4.40) 



where the values 28.347 and 45.226 figuring in the ranges for ao were computed numerically 
and the function A Q (ao) has no analytic expression such that, for each particular value of ao, a 
numerical evaluation must be carried out. The values of this functions for some values of ao are 
listed in table El 



Q: 



A a (a ) 



30 



50 



10 2 



10 3 



10 5 



+oo 



0.8935 0.8887 0.8848 0.8812 0.8808 0.8808 



Table 7: Values of the functions A a (ao) f or several values of ao. 



Hence, resuming the final results with respect to the bounds on the coefficients ao and a\ for the 
radial coordinate dependent exponent a{r\) = ao + ol\ 2GM/(c 2 n) (|4.20p . we have obtained the 
properties of space-time described by the locally anisotropic metric (]4,4j) for the following ranges 
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of these coefficients 

q > 3 , qiG -4 (3) (ao) : 

only contribution to singularity at n = is given by SC mass pole , 

mass — energy density is positive definite for r% > rx.sc , 

space — time is singularity free at r\ = ri.se \ 

(4.41) 

a > 5 , a\ € ^ (5 ^(ao) : 

only contribution to singularity at n = is given by SC mass pole , 
mass — energy density positive definite for r% > ri.se j 

space — time is asymptotically flat at r\ = ri.se • 

Hence we have finished to build our ansatz. As a final remark we note that when considering 
relatively large distances when compared with the SC radius (n 3> ri.se > which is the case even 
at planetary scales), for most purposes the parameter ax works simply as a regulator of the 
singularities at the origin and it can be considered as close to null as wished such that its effects 
can safely be neglected away from the origin of the coordinate frame. We will take this approach 
in the following analysis. 

Next we will analyze in detail the modification to the General Relativity Newton law due to the 
ansatz for the locally anisotropic metric (|4.4p . 

5 The Modified Newton Law 

The main objective in this section is to compute the corrections to the General Relativity Newton 
law corresponding to the locally anisotropic metric (|4.4p and analyze the physical effects due to 
these corrections. We compute the Newtonian limit for the physical acceleration acting in a test 
mass from the geodesic equations considering a massive body with center of mass at the origin 
of the coordinate frame independently of considerations concerning corrections to the observables 
due to the location of the observer which we will address somewhere else |45j . In the following 
we employ spherical non-expanding coordinates rx , as discussed in section 12. 3| these coordinates 
correspond to the physical spatial lengths. 

5.1 Radial Gravitational Acceleration due to a Central Mass 

In this section we consider the gravitational acceleration of a test mass in the gravitational field 
of a central massive object. We assume the limit of non-relativistic velocities, hence neglect 
the contributions due to the test mass velocities such that the only non-null component of the 
acceleration corresponds to its radial component, fx ~ — c 2 ]? 1 ^ = — c 2 (( a )r 1 00 + Ar 00 ), where 
the connection is given in equations (|D.2|) and (jD.lip of appendix ffD|) . 
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The corrections to the General Relativity Newton law contain powers of the Hubble rate H (on 
H 2 and H 4 ), hence we consider a decomposition of this modified Newton law into three distinct 
factors 

h = F GR + F H 2 + F h a , (5.1) 

where Fqr is the usual expression obtained within General Relativity for flat Minkowski back- 
ground (|3.10p . F H 2 the second order correction in the Hubble rate (~ H 2 ) and F h a the fourth 
order correction in the Hubble rate (~ H 4 ). Specifically these factors are 



For 

F H 2 



GM 2{GM) 2 



cr r 



-n 1 



2GM\ a 

c 2 n / 



(1 - a)GM 



c 2 ri 



(1 + 



2GM\ 

C 2 T\ J 



log 1 



c 2 n / 

2GM\ , 
~ — a 

& n j 



H 



(5.2) 



2GM\ 

c 2 n / 



2a- 1 



2aGM 



+ i 



2GAA 

c 2 n / 



2 + n log l 



2GM\ 



c 2 ri 



where we have replace the derivatives of the scale factor by the time dependent Hubble rate 
H = a/a and deceleration parameter q = —a/(aH 2 ). Next we consider and discuss both the 
perturbative regime and non-perturbative regime with respect to the gravitational field Use- 



5.1.1 Perturbative Regime 

Commonly, when deriving the classical Newtonian limit for General Relativity quantities, it is 
assumed the weak field approximation f|3.7[) and considered an expansion on the gravitational 
field Use = 2GM/(c 2 ri). Here, when applicable, we consider a third order expansion on the grav- 
itational field Use such that the several powers of (1 — Usc) p an d the logarithm are approximated 
by 

/ 2GM\ P 2GM 1 . . [2GM\ 2 

\ c z r\ J & r\ 2 \ c z r\ J 

1 , w ,/2GM\ 3 / 4 /2GM\ 4 \ , 

- 6 P(v-^-^{ 7 r K ) +°^ 4 (^rj j ■ (5-3) 

\ c z r\ J 2 \ c z n J 6 \ c z r\ J \\ cr i/J 

Generally the first of these series is convergent independently of the value of the exponent p (as 
long as it is finite). The exponent p = a = ao + aiUsc (|4.20p is generally greater than unity so 
that to attain any desired accuracy for the series expansion it may be necessary to consider more 
than the lower order terms, often it is simpler to directly evaluate the exact expressions. With 
respect to the first order series expansion we note that is only a valid approximation when the 
first order term is less than unity, i.e. for radial distances greater than the Schwarzschild radius 
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times the exponent p 

2a GM 

n > — 2 — • ( 5 - 4 

cr 

Hence, assuming a small parameter ot\ and interpreting this result only with respect to the value 
of the parameter ao, the series expansion is no- longer a valid approximation for large values of 
this parameter, oq > c 2 ri/(2GM). In the following discussions we will distinguish between the 
perturbative regime and non-perturbative regime in the gravitational field which depends both on 
the value of the radial coordinate r% and the value of the parameter ao in the region of interest. In 
the perturbative regime we will be at most working to order rj~ 2 in the radial coordinate, however 
it is necessary to consider an expansion to one higher order to correctly compute the coefficients of 
the several powers of r\. This is due to the connection T^q, within the same order in the Hubble 
rate powers, having the factors with a powers multiplied by distinct powers of n, specifically r\ 
and r^ 1 . 

Hence, considering that the above bound (|5.4p is obeyed, we can consider a series expansion of 
the modified Newton law (|5.1[) such that the factors Fjj2 and Fjji (|5.2p are approximate by 



F H2 « - q r l H 2 + (a + l) q ^H 2 + {l + ±a 1 -al)(q + l)^^H' 1 

c z 2 cr n 



l H l H \ C L H 



In these series expansions we have explicitly considered the radial coordinate dependent exponent 
ot{r\) = «o + 2ai GM/(c 2 rj.) (14.201) . The last expression is written in terms of the Hubble length 
Ijj = c/H due to these terms being only relevant for large radial coordinate r\ ~ Ijj. As for the 
factor -Fgr is exact within the framework of General Relativity (in the non-relativistic velocities 
approximation) being already a series expansion on the gravitational field 

The lower order terms —qr\H 2 and —c 2 rf/ljj in the expansion of the factor F H 2 and F h a cor- 
respond to the Newton law for the expanding background (|2.53p described by the FLRW met- 
ric (I2.28p . The first of these terms correspond to a repulsive gravitational interaction (we recall 
that qo < 0) which becomes dominant with respect to the classical attractive gravitational inter- 
action, Fgr ~ —GM/r 2 , for the value of the radial coordinate r*/ pert ^ 

=\ qo \ r \H 2 o r\ = , (5.6) 

V r l(pcrt)J VI 90 1 -Ho/ 

where we derive these expressions assuming today's value of the Hubble rate Hq (jl.ip and de- 
celeration parameter go 02]). We note that the value r*/ pert N obtained perturbatively is a good 
approximation as long as we assume a small negative parameter a\ ~ (as discussed in the pre- 
vious section 14.41 this is enough to regularize the singularities at the center of mass maintaining 
the Schwarzschild mass pole) and a relatively small parameter «o such that the higher order in 
the series expansion (I5.2p are negligible. 
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Figure 4: The modified Newton law r\ = Fqr + F H 2 + F H 4 (TO]) /or the Sun (M = 1.98 x 10 30 Kg) and 
the large galaxy UGC2885 (assuming as approximation a point-like core mass of M = 2 x 10 12 Ms U n = 
3.96 x 10 42 Kg) for ao = 10 such that the perturbative regime in the gravitational field is valid. The lower 
order terms in the series expansion 115.5]) of \5.2\) are dominant for each of the factors, For ~ —GMjr\, 
F H 2 w —q Q riH^ and F H i w ~^ r \I^H- F° r relatively small scales the factor For is dominant, in the 
intermediate region the factor F H 2 is dominant and for relatively large scales the factor F H i is dominant 
such that the effects due to the central mass are only felt at relatively small scales: 

(a) relatively small values of the radial coordinate above the Schwarzschild radius r%,sc = 2940 m for the 
Sun, the modified Newton law is the continuous line such that the gravitational acceleration changes sign at 
r\ — r\ — 3.51 x 10 18 m \5.6\) . the factor Fqr is the dashed line, the factor F H 2 is the dashed-dotted line 
and the factor F H i is negligible; 

(b) relatively large values of the radial coordinate up to the cosmological horizon r% = lg = 1.31 x 10 26 m 
for the Sun, the modified Newton law is the continuous line such that the gravitational acceleration changes 
sign at r% = \J — q In = 10 26 to {2.5$ , the factor Fqr is negligible, the factor Fg2 is the dashed-dotted line 
and the factor Fgi is the dashed line; 

(c) relatively small values of the radial coordinate above the Schwarzschild radius ri se — 5.88 x 10 15 m 
for the large galaxy UGC2885, the modified Newton law is the continuous line such that the gravitational 
acceleration changes sign at r\ = r\ = 4.42 x 10 22 m \5. 6\) . the factor Fqr is the dashed line, the factor 
Fjj2 is the dashed-dotted line and the factor Fjja is negligible; 

(d) relatively large values of the radial coordinate up to the cosmological horizon n = Ir = 1-31 x 10 26 m 
for the large galaxy UGC2885, the modified Newton law is the continuous line such that the gravitational 
acceleration changes sign at r\ = y/—q Ir = 10 26 m \2.5J$ , the factor Fqr is negligible, the factor Fgs is 
the dashed-dotted line and the factor Fh* is the dashed line. 
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Hence for t\ < r i( port ) the General Relativity Newton acceleration factor -Fqr is dominant such 
that the gravitational interaction is attractive, at n = ?"iV pert ) both the factors Fqr and F H 2 terms 
cancel such that the gravitational interaction is null and for r\ > r\ r peli \ the factor F H 2 is dominant 
such that the gravitational interaction is repulsive. As for large distances, at n = \/—qo Ih (|2.54p . 
the gravitational acceleration is null and, for r\ > yj—q® Ihi the gravitational interaction becomes 
again attractive due to the term —c 2 rf/ljj in the factor F h a becoming dominant as already 
discussed in section 12.31 It can also be concluded straight forwardly that, although the value of 
r i(pert) Q5.6P increases with the value of the mass, for large enough radial coordinate the mass effects 
are negligible such that, independently of the value of the mass M, the modified Newton law (|5.ip 
converges asymptotically to the Newton law for an expanding background (|2.53|) described by 
the FLRW metric (|2.28p . We recall that this result is consistent with the assumptions considered 
when constructing the ansatz for the locally anisotropic metric (|4.4p in section 14.11 

Therefore, for small values of the parameter ao, we obtain the lower order perturbative modified 
Newton law (the weak field classical limit of the relativistic acceleration) interpolation between 
the General Relativity Newton law (|3.10p obtained from the Schwarzschild metric (|3.2p and the 
Newton law for an expanding background (12.53P obtained from the FLRW metric (12. 2p . As 
examples of typical profiles of the perturbative Newton law for a point-like central mass are 
plotted in figure [4] the gravitational accelerations felt by a test mass in the gravitational field of 
the Sun and the large galaxy UGC2885 for ao = 10. 



5.1.2 Non-Perturbative Regime 

For higher values of the parameter ao the remaining terms on the series expansion (on the gravita- 
tional field) become relevant and contribute significantly to the modified Newton law. Instead of 
directly consider the validity limit imposed to the first order term by equation (|5.4p obtained from 
the series expansions of the individual multiplicative terms (15.3P of the factors Fjj2 and Fjja, (|5.2p 
we can carry out a similar analysis directly in the series expansions for these factors (|5.5p . Hence, 
for spatial scales for which the factor F H 4 is negligible, by directly comparing the terms on the 
series of the factor F H 2 (15. 2p . we can conclude that, for some given n, the series is no- longer 
accurate when the first order term in ao, (a + l)go GMHl/c 2 , is of the same order of the zeroth 
order term in ao, —qo T\ Hq. When the parameter ao is above the value that satisfies this equality 
the series expansion is no longer valid in the neighborhood of the given value for r\ and the exact 
expressions for F H 2 and F H 4 (|5.2p must be considered. In particular, for large values of the param- 
eter ao, the exact value of the radial coordinate r\ for which the gravitational acceleration is null is 
above the value given in equation (15, 6h . r* > r*( pert A, and must be computed by solving the exact 
expression for the modified Newton law. For the specific case of n in the neighborhood of rj/pert) 
the zeroth order term — <Zo r i( pcrt ) Ho coincides with the classical Newton law term GM/r^ pert ^, 
hence we can directly compare this term with the first order term (ao + l)qo GM H 2 /c 2 such 
that the specific value of the parameter ao = aj separating in between these perturbative and 
non-perturbative regimes in the gravitational field is 

GM _ GM 2 ^ c 2 -\g \l{GMH )l _ c 2 



[r 



l (pert) J ^ \qo\s(GMH )s \q \3(GMH ) 



As examples are plotted in figure[5]the values of r\ computed numerically from the exact expression 
for the modified Newton law (|5.ip as a function of the parameter ao for the Sun and the large 
galaxy UGC2885 as well as the respective values aj above which the perturbative regime in the 
gravitational field is no longer valid. 
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Figure 5: The values of r* for which the modified Newton acceleration \5. 1\) is null computed numerical 
from the exact expressions of Fgr, Fg2 and Fjj4 115.2]) for the Sun (Ms-an = 1-98 x 10 30 Kg) and the large 
galaxy UGC2885 (assuming a point-like core mass of M — 2x 10 12 Ms U n = 3.96 x 10 42 Kg): 

(a) the case of the Sun for which r*i vert \ = 3.51 x 10 18 m and — 2.39 x 10 15 ; 

(b) the case of the large galaxy UGC2885 for which r* { t) — 4.42 x 10 22 m and a* Q — 1.50 x 10 7 . 



We further note that, the higher the value of the parameter ao, the closer the modified Newton 
law (|5.ip is from the General Relativity Newton law for small values of the radial coordinate 
(roughly for n < r \ / ert ) ) . This result is consistent with the results derived in section 14.21 when 
analyzing the curvature close to the Schwarzschild horizon, for high values of the parameter ao 
space-time near this horizon is flat. Furthermore, for large enough values of this parameter, 
there is an attractive acceleration contribution that becomes significant for values of the radial 
coordinate in between r*/ peTt \ and r\ increasing the Newtonian acceleration towards the central 
mass. In figure [6] are presented the examples of the modified Newton law (|5.ip for the Sun and 
the large Galaxy UGC2885 for values of the parameter ao above ctg (j5.7j) . 

In particular these results imply that, for large values of the parameter ao, while for planetary 
scales the effects of expansion may be for most purposes negligible, for galaxy scales (and above) 
we may obtain significant deviations from the General Relativity Newton law. The most straight 
conclusion is that exists a well defined cutoff corresponding to r\ = r\ such that for radial coordi- 
nate above this value the gravitational interaction is repulsive due to the expanding background 
and no stable orbits exist. In order to further analyze large scale (meaning at least galactic scale 
and r\ < r*) orbital motion, let us assume the approximation to circular orbits such that the 
orbital speed is as usual computed by noting that the gravitational acceleration exactly matches 
the centrifugal acceleration 



-h = — <^ v orh = y/-fi n . (5.8) 

n 

When the acceleration is given by the modified Newton law (|5.ip the orbital velocity will signifi- 
cantly deviate from the orbital velocity for the General Relativity Newton law (|3.10p for values of 
the radial coordinate above J*5V pert v As an example are presented in figure [7] the rotation curves 
for several values of the parameter ao f° r the large galaxy UGC2885 assuming (as a simplification) 
the core to be a point-like central mass. 

Hence in this section we have shown that the modified Newton law corresponding to the locally 
anisotropic metric (|4.4p with exponent a given in equation (|4.20p interpolates between the General 
Relativity Newton law (|3.10p for relatively short radial scales and the Newton law for the expanding 
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Figure 6: The modified Newton law r\ — Fqu + Fjj^ + F#a 15. 1\) for values of the parameter ao above 
a o j5.7\ ) for which the perturbative regime in the gravitational field is not valid. Are shown the cases 
of the Sun (Ms U n = 1-98 x 10 30 Kg) and the large galaxy UGC2885 (assuming a point-like core mass of 
M = 2 x 10 12 Afsun = 3.96 x 10 42 Kg ) with the General Relativity Newton law For represented by the dashed 
lines, the contribution of the factor F^2 is represented by the dotted lines, the factor Fjj4 is negligible and 
the modified Newton law r\ h5.1\) is represented by the continuous lines: 

(a) for the Sun with a = 10 4 , 10 15 , 5 x 10 15 , 10 16 ; 

(b) for the Sun with a a = 7 x 10 16 , 10 17 ; 

(c) for the large galaxy UGC2885 with a Q = 10 4 , 10 7 , 5 x 10 7 , 10 8 ; 

(d) for the large galaxy UGC2885 with a = 2 x 10 s , 3 x 10 s . 



background (|2.53|) for relatively large scales. In between there are significant deviations from 
these limiting laws, in particular for large values of the parameter ao there is an increase of 
the attractive gravitational acceleration towards the central mass near the value of the radial 
coordinate for which the net gravitational acceleration is null. This result is not completely 
unwelcome and could contribute, for instance, for the flattening effect of the rotation velocity 
curves for galaxies [23]. We remark that a large value of this parameter is physically justified by 
noting that it maintains the usual General Relativity Newton law at planetary scales such that 
space-time is approximately flat near massive bodies. 

In order to properly evaluate the effect of the locally anisotropic metric (|4.4p for planetary motion 
and analyze if it is possible to obtain any bound on the parameter ao for the solar system we 
next analyze the orbit solutions for this metric. We start by taking the most simple approach 
by considering circular orbits which allows to estimate the time dependence of the orbital radius 
due to the expanding background. Then we proceed to compute the elliptic orbit solutions to 
zeroth order on time (static solutions) which allows to obtain estimative for precession and period 
corrections due to the expanding background. 
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Figure 7: Circular orbital velocities v ot \, — \J—(Fgr + Fjj2 + F H 4.)r\ h5.1\) for values of the parameter 
ao = 10, 10 7 , 10 8 , 2.5 x 10 s above «q = 1.5 x 10 7 (J.Tp for the large galaxy UGC2885 (assuming a point-like 
core mass of M — 2 x 10 12 Msun = 3.96 x 10 Kg). The orbital velocity for the General Relativity Newton 
law 13. 1 0\) is represented by the dashed line. 



5.2 Circular Orbits Approximation: Time Varying Orbital Radius 

As a first approximation to orbital motion we are considering circular orbits and compare the 
effects of the modified Newton law (15. If) obtained from the locally anisotropic metric (14. 4p . We take 
as starting point conservation of angular momentum. For orbital motion the angular momentum 
is given by the constant of motion J = —rfdip/dr (|3.16j) discussed in section [3~2l when reviewing 
orbit solutions in the framework of General Relativity. We note that the expression for this 
constant is the same both for classical Keplerian orbits, relativistic orbits obtained from the 
Schwarzschild metric (|3.2p and orbits in an expanding background obtained from the locally 
anisotropic metric (14. 4j) as derived in appendix ID.3I Hence, recalling that the orbital velocity 
for circular orbits is generally given by equation (15. 8p and taking the usual definition of angular 
velocity (fi = uj = r\ with -u or b = \J — r\ ¥[ (|5.8p , we obtain the following definition of J c i rc for 
a circular orbit of radius ri. or b 



•4rc ~ -7 2 »"l h 



ri=ri.orb 



GMn. olh -^-F H >r* 



I 2GM ffr Lolb / ^ 2GM 



cf - n- l.orb 

*7 ^ITv 7 (5.9) 



Here 7 = dt/dr is the relativistic factor for the generic metric (|4,4p and we are taking the limit of 
non-relativistic velocities x M <C c for which 



d 2 r\ 
~d^ 



7 2 h 



n 

: 7~2 T^2 (5.10) 

5oo-2 5 oi^-5ii ( r i) -922 (f) 



l _2GM_fHr 1 f ( l _2GM\ a 
c^ri \ c I \ r\ j 



65 



In (|5.9p we are considering the modified Newton law (|5.ip and neglecting the term F h a which is 
only relevant for cosmological scales {r\ ~ Ijj ~ 10 26 m) as discussed in the previous section (here 
we are considering planetary scales, at most r\ ~ 10 13 m) such that the gravitational acceleration 
is given by r% « Fqr + F H a with the factors Fqr and F H 2 given in equation (|5.2p . In the 
last approximation of equation (|5.9j) we have expanded the denominator to lowest order in the 
gravitational field and H 2 . 

For circular orbits, the main effect obtained due to the corrections of the expanding background 
is a time varying radius. Recalling again that J is a constant of motion (which corresponds to 
angular momentum) and that the Hubble rate H(t) ~ Hq — qo Hq t ()2.5[) is time dependent we 
note that necessarily some other parameter of the equation (|5,9p must be varying with time to 
maintain J constant. Assuming that the Newton constant is fixed on time we are left only with 
the possibility of a time -varying radius r^ofb — ^*i.orb 

(t). Hence differentiating equation (|5.9I) and 
solving the resulting equation J c i rc = for ri. or b we obtain, to lowest order, the time dependence 
of the orbital radius 



n.orb 



n.orb 
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where we have neglected the contribution from the parameter a\ and, in the last approximation, 
have expanded the denominator to lowest order and kept only the terms on Hq. We note that 
this expression increases with the orbital radius ri orb and decreases with the mass M. Hence, 
for fixed value of the parameter ao, the values for this expression will change significantly for 
distinct orbits. We note that this result is consistent with the expected results for the FLRW 
metric (|2.28p . the expansion effects are larger for larger distances and taking the interpretation of 
expansion being due to the net effect of the long range gravitational interactions we also expect 
that in the neighborhood of a larger mass M its effects will be dominant with respect to the 
effects of more distant objects wile for smaller masses the net effect of distant objects will become 
relevant. With respect to the value of the parameter ao, the time variation of the orbital radius 
will be positive with a maximum for the lowest allowed value of the parameter ao = 3, will decrease 
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with increasing ao until reaching a negative minimum and them increase towards zero (being still 
negative, it is only null in the limit ao — * +00). We note that both the values of the parameter 
ao for which ri. or b = and for the negative minimum of ri.orb are below ckq (|5.7j) although being 
close to this value, hence still for values of ao f° r which a series expansion in the gravitational field 
can be considered. However due to being relatively close to this value to attain accurate results 
perturbatively it would be necessary to consider higher order terms of the series. Hence we will 
evaluate numerically the exact expression (|5.1ip in the following discussions. As an example of 
the typical values of ri.orb as a function of the parameter ao are presented in figure [8] the values of 
"Ti.orb/fi.orb for the Moon orbit with the central mass being Earth and for Venus, Earth and Mars 
orbit with the central mass being the Sun. 



(a) Earth-Moon 
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Figure 8: Examples of the profiles of the time variation rate of the orbital radius ri.orb /ri.orb 115.11)) 
function of the parameter ao assuming the circular orbits approximation: 

(a) for the Earth-Moon orbit, the maximum positive variation is ri.orb Ai 
sponding to ao = 3, the radius variation is null ri.orb /ri.orb = for ao = 1.72 x 10 11 and the 
negative variation is ri.orb/ri.orb = —2.17 x 10 _44 s _1 corresponding to ao — 2.40 x 10 

(b) for the Sun-Venus orbit, the maximum positive variation is ri. rb/ri. rb 
ing to ao — 3, the radius variation is null ri. rbA*i.orb = Os^ 1 for ao = 
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(c) for the Sun-Earth orbit, the maximum positive variation is ri. or b/ri. orD 
ing to ao — 3, the radius variation is null r l rb/ri. rb = 0s _1 for a = 
negative variation is r l or b/r l or b — —3.86 x 10~ 42 s _1 corresponding to ao 

(d) for the Sun-Mars orbit, the maximum positive variation is ri. or b/ri. or b = 7.34 x 10~ 40 s _1 correspond- 
ing to ao = 3, the radius variation is null ri.orb/ri.orb = 0s _1 for ao = 
negative variation is ri.orb/ri.orb = —1-36 x 10~ 41 s _1 corresponding to ao 
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Depending on the value of the parameter ao, the radius variation of the planetary orbits in the 
solar system can be, either positive, either negative or some positive and others negative. For low 
values of the parameter ao <C 10 8 the radius variation are all positive, while for high values of the 
parameter ao 3> 10 8 the radius variation are all negative and for values of the parameter ao ~ 10 8 
the radius variation for planets closer to the Sun are positive and for planets farther away from 
the Sun are negative. As an example, taking ao = 10 8 , we have that Mercury and Venus are 
approaching the Sun while the remaining planets are drifting away from the Sun. We note that 
this result is not inconsistent with global spatial expansion, we recall that for large spatial scales 
we recover the FLRW metric (|2.28p and the usual properties of spatial expansion, instead this 
behavior should be interpreted as due to the local background deformation in the neighborhood 
of the central mass M. 

In order to compare the values obtained for the expression (|5.1ip with experimental data we note 
that the only available experimental results are the values for the time-varying Newton constant 
G/G (see for instance |35| for a extended review in this topic). Similarly to the derivation of 
equation (|5.9p . for the General Relativity Newton law we have that the constant of motion J for 
a circular orbit is given by 

j 2GM 



Jlc.GR « -7Mn| ri=ri . orb . GiJ « GMn. olh . GR - c y°^ GR = GMrx.ato.GR , (5.12) 

where, as before, the non-relativistic velocity limit has been considered, however no expansion on 
the gravitational field has been considered such that, up to Special Relativistic corrections, this 
result is exact. Hence for a static background no time variation of the orbital radius is possible 
unless either the Newton constant G or mass M are also varying on time. Assuming that the 
mass is not varying on time, M = 0, and that both the radius and the Newton constant are time 
dependent r\. a to.GR = r i.orb.Gi?(^) an d G = G(t), differentiating the equation (|5.12p with respect 
to time, we obtain the following equality 

G fl, mb .GR . s 

r = -~r • ( 5 - 13 ) 

^l.orb.GR 

Therefore when analyzing experimental data (as well as in numerical calculations) employing 
General Relativity in a static background, a varying radius r\. oro .GR is equivalent to a varying G 
which is quoted for several systems in the literature. In particular we list the results obtained 
for the moon Venus [Uj and Mars [48] and compare them with the respective allowed values 
obtained from equation (|5.1ip and presented in figure [8] As far as the author is aware there are 
no estimative based on Earth's orbital motion [35]. Although within the error bars the theoretical 

Orbits h.orb.GR/ri.orb.GR = ~G/G n.orb/Vi.orb 

Moon-Earth —G/G = -(1.9 ± 2.22) x 10" 20 G [-2.17, 1.17] x 1(T 42 
Sun- Venus \G/G\ < 4.76 x 1(T 18 G [-0.0224, 1.2] x 10" 41 

Sun-Mars -G/G = (6.34 ± 31.71) x 10~ 20 G [-0.136, 7.34] x 10~ 40 

Table 8: Experimental orbital radius variations and allowed theoretical range obtained from equation H5.ll]) . 

ranges obtained from equation (|5,lip are clearly many orders of magnitude below the experimental 
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values. This comparison is not conclusive due to the error bars being significantly larger than the 
quoted experimental values and, depending on the sources and methods employed in the data 
analysis, the values are quite distinct. Hence no bounds on the parameter ao can be drawn from 
this analysis. 

It is also relevant to remark that usually it is assumed that a varying G (equivalent to a varying 
orbital radius in GR) is independent of the orbit parameters, this is no longer the case when 
considering the locally anisotropic metric (|4.4p . equation (|5.1ip does depend non-linear ly both in 
the central mass value M and the orbital radius ri, or b, hence it is expected that distinct values 
are obtained for each orbit (as opposed to GR) which may justify the dispersion of experimental 
values in the literature. In table [9] we list the allowed ranges for the remaining planets of the solar 
system. We note that for larger orbital radius the radius variation significantly increases. 

As a final comment we note that from the locally anisotropic metric (14, 4h it is also possible to 
consider either a varying G with fixed orbital radius ?"i or b ° r both varying G and orbital radius 
rx.oib.GR- These approaches may be justified in the context of varying fundamental constant theo- 
ries [35j or extended theories of gravity [49, 5U1G1]. Assuming conservation of angular momentum 
we obtain the generic expression for the time variation of the radius when G / 
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(5.14) 

where ri.orb/ r l.orb|(j = o ^ s given in equation (|5.1ip . It is possible to maintain the orbital radius 
fixed on time by fine-tuning the time variation of G such that this equation is null, ri.orb = 0. 
Otherwise if the time variation of G is given by some extended theory of gravity its effects simply 
add to the orbital radius variation. We remark however that in the present framework, for which 
we consider only the expanding background, our previous analysis is consistent and does not 
require a varying G. 
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5.3 Perturbative Static Elliptical Orbit Solutions 



With the objective of estimating the orbital precession and orbital period corrections to an ellipti- 
cal orbit due to the expanding background we will proceed to derive orbital solutions for the locally 
anisotropic metric (|4.4p . It is hard, if not impossible, to obtain a analytical solution considering 
the differential equations for a time varying Hubble rate H. The main difficulty is that energy 
conservation is no-longer given by a constant of motion, instead we have a non- linear second order 
differential equation on the function t((p) coupled to the differential equation for u(<p) = 1/ri (<£>). 
Hence, for technical simplification purposes, we are taking the static orbit approach by consid- 
ering a fixed Hubble rate H = Hq. Also we note that a estimative for the orbital radius time 
dependence has already been computed in the previous section considering the approximation to 
circular orbits. 

The differential equation describing an static orbit of a test particle in the gravitational field of a 
central mass M for the locally anisotropic metric (|4.4p . to order Hq, is given in equation (|D.26|) 
of Appendix ID.3I for the inverse radial coordinate function u(<p) = l/n(<p) and fixed angular 
coordinate 6 = tt/2. We note that for planetary scales the corrections to order Hq are negligible. So 
far the author failed to obtain a treatable analytic solution to this differential equation considering 
the exact expressions for the corrections due to the expanding background. Here, when refering 
to corrections we mean with respect to the General Relativity orbit equation (13.191) discussed in 
section 13.21 For orbits in the solar system the function u has small values (0.5 x 10~ 12 < u < 
0.5 x 10 m , where ri. or b is the orbit semi-major axis) such that we can consider a series 
expansion on u of the corrections due to the expanding background. However we note that this 
expansion is equivalent to an expansion on the weak gravitational field and, as already discussed 
in section [5.1 .11 this perturbative regime is valid only for values of the exponent parameter ao (we 
are assuming a negligible negative valor for the parameter a\) up to (|5.4p 



Above this value it is either necessary to consider higher order terms on the series expansion or 
to consider the exact expressions. Nevertheless we remark that, for a fixed value of the radial 
coordinate r\, and larger values of ao > ao.max.pert the corrections given by the exact expres- 
sion due to the expanding background will decrease significantly in absolute value becoming, for 
very large values of the parameter ao ^> ao.max.pert) negligible. Considering a rough numerical 
estimative we conclude that the maximum deviation for the corrections for higher values of the 
parameter ao > ao.max.pert is less than 20% of the value obtained in the perturbative regime for 
ao = ao.max.pert- As an example are plotted in figure [9] the values of the exact and perturbative 
correction terms on Hq as presented in the differential equation (|D.26p for the case of Earth's 
orbit. 

Hence, for technical simplification, in order to approximately estimate the orbit precession and 
period corrections due to the expanding background and compare our results with the same effects 
due to standard General Relativity we are considering a perturbative series expansion on u to the 
orbit differential equation (ID.26j) of Appendix lD.3[ We note that the General Relativity corrections 
to the orbit differential equation are second order in u (|3,19p . To match at least the same order 
of these corrections we are considering a series expansion up to second order such that we obtain 
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Figure 9: Plot of the exact (dashed line) and perturbative (continuous line) expressions for the corrections 
to the General Relativity orbit differential equation due to the expanding background given in equation \D. 26V 
as a function of the parameter ao for Earth's orbit. The perturbative regime is valid up to ao. max. pert = 
c 2 ri. orb /(2GM) sa 0.5 x 10 7 fQ3U : 

(a) plot of the exact expressions up to ag = 10 9 , the corrections asymptoticaly vanish for large Oq 3> 

C^O. max. pert j 

(b) plot of the exact and perturbative expressions up to a = 10 s , the perturbative and exact expressions 
approximately match up to ao. max. pert- 



where the correction coefficients are 
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(5.17) 
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As an example of the magnitude of the several correction terms for Earth's orbit we have that 
the constant term in the right-hand side of the differential equation is of order ~ 7 x 10 



-12 



m 



the General Relativistic correction is of order ~ 2 x 10 



-19 



m 



and the corrections on H% are of 



order ~ 9 x 10 34 m 1 for ao = 3 and ~ — 3 x 10 41 m 1 for ao = 10 9 . 

To solve the differential equation (|5. 16|) we take the same approach discussed in section 13.21 We 
start by solving the differential equation considering only the dominant term in the right-hand 
side of (|5.16p obtaining 



GM o 

u o.m + A u o.hi = -jr B 



U 0.H 2 



1 + e cos(\/ r A(p) 



d 



(5.18) 



where d is defined in terms of the ellipse semi-major axis ri, or b and the eccentricity e, being related 
to the constant of motion J and the correction coefficients A and B by the following equalities 



d = n.orb(l - e 2 ) 



J 2 A 



GM B 



— 44> J 



r Lorb GM(l-e 2 )^ 



(5.19) 
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We note that, besides this correction to the relation between the constant of motion J and the 
parameter d given by the ratio of the correction coefficients B /A (|5.17p , there is a small precession 
effect proportional to Hq due to the factor y/A in the cosine argument of the solution (|5.18p . We 
will deal with this effect in detail later. 

Next let us compute the corrections to the solution u jj2 by considering the remaining terms in 
right-hand side of equation (|5.f6p evaluated for the function u H 2 such that the full solution is 

u = u . H 2 + u GR . H 2 + u H 2 . (5.20) 



Here the functions ji GR fi 2 and u^j2 correspond respectively to the (modified) General Relativ 

T 2 



ity and expanding background corrections approximated to order Hq being the solutions of the 



following differential equations 

U GR.H* + A U GR.H* = —3— C u 2 0H2 



3GM „ 2 



3GM (l + ecos(VI^)) 2 
I 2 

2 1 / TT \ 2 



/, „ (HoY I ( Ho\ GM 1 

*" +AU " = -{-J ^-2 +a °{-J —^2 (5 ' 21) 



H \ 2 J 



-ao 



J J (I + ecos(VAip)) 3 
H \ 2 GM 1 



J J c 2 (l + ecos(VA<p)) 2 



such that we obtain 



U GK.H 2 = ^ ( ( X + y ) ~ ~Q cos ( 2 v 7 !^) + v 7 ! e if sin(\/I <p)\ 



d 3 H2 f ao GM (-4 + e 



^J 2 (f-e 2 )l c 2 d 4(f-e 2 )(f + ecos(^^)) ( 5 - 22 




3 _ a GM\ 2e arctan I \l ttt, tan 
2(f-e 2 ) c 2 d ) 

Where acn = 3GM/ (c 2 d) (|3.28p and we have set the integration constants such that no cos(vG4 cp) 
neither sin(V^4y) are present (see the discussion of equation (|3.27|) for the justification of this 
choice). The modified General Relativity solution u GRH 2 has the same structure of the usual 
solution ugr ()3.26p . the first term is a constant that can be neglected, the second term has a 
period that is a multiple of the one of solution u H 2 (|5.f8p contributing a small shift to the 
orbital period and the last term contributes to the orbital precession. As for the solution Ujj2 has 
a similar structure, the first term is a constant that can be neglected, the second term contributes 
a small shift to the orbital period and the last term contributes to the orbital precession. This last 
result is justified by noting that the analytic continuation of the inverse of a function corresponds 
to the argument of the function (in this way arctan(tan ip) = ip increases monotonically with 
tp), here due to the correction being small compared to the dominant term we will expand the 
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functions to lower order as we did in equation (|3.29[) . Therefore we are considering the following 
series expansions for the oscillatory terms of u H s (|5.18p . the last term of «gr.h 2 (|5.22p and the 
last term of Ujj2 (|5.22p 



2 2 

cos(VI^) « l-A 2 ^ + 0((/) , ffan(VAip) « 2\/I^ + 0(y> 4 ) 
arclan ( Ian ( ^ ) | ~ .1 \/ y + C% 4 ) . 




(5.23) 



Hence, neglecting the constant terms in the solutions u GRH 2 and u H 2 (|5.22p . gathering the several 
terms and respective coefficients of the series expansions (|5.23p and approximating these by a 
cosine function as we did in equation (|3.29[) we can rewrite the full solution u (|5.20p as 

If ( ( A(/? G R A<p H i\ \\ 
U " ^ ( 1 + e COS ( ( 1 ^ TT - ) V J J + ""osc.GR + W osc . H 2 , 

2tt 2 (i_ e )(i + e )| \c 2 d 2(1 -e 2 )GM J 

(5 c -6 A )a GR 2 ,„ x . /, « 2 -4 + e 2 + 3e 2 cos(2(^) , 4, 

« = 6 e cos(2^) + (M 4(1 _ e2)2gM (1 + ecosy) +°^) 

where A</?gr/(27t) = «gr is the precession per turn of the orbit due to General Relativity cor- 
rections given in equation (|3.31|) and Aip H 2 / (2ir) is the precession per turn due to the expanding 
background. The factor tt OS c.GR is the oscillatory factor of the correction to the orbit solution 
due to General Relativity given in equation (|3.3U|) and u osc ^2 is the oscillatory factor correction 
to the orbit solution due to the expanding background. For both these factors we considered 
the approximation \^Aip ~ ip in the argument of the cosines. This approximation is justified by 
noting that 5a (|5.17p is a small perturbation (A = 1 + 5 a) such that over one turn of the orbit its 
contribution to the period correction is negligible. In deriving these expressions we have used the 
following equality for the constant of motion J 2 = GM d A/ B (|5.19p and considered an expansion 
on 5a, 5b and 5c (I5.17p . keeping only the lower order terms in Hq. 

In order to compute the period correction to the orbits due to the expanding background we use 
the same method of equation (|3.33p discussed in section 13.21 Then the period is 



T = T + AT G r + AT H 2 , AT H2 = ~ C" dtp ^1 ( X _ ^£ Uq ) , (5.25) 

\J\ Jo u o \ c 



where no corresponds to the classical Keplerian orbits solution (|3.2ip . ATqr is the period correc- 
tion due to General Relativity given in equation (|3.33p and AT H 2 is the period correction due to 
the expanding background with u osc H 2 given in equation (|5.24[) . 

In table [9] are listed the estimative for orbital precession and period corrections both for the 
General Relativity corrections and for the expanding background corrections in the solar system. 
As can readily be concluded the effects due to the expanding background are several orders of 
magnitude lower than the respective GR effects (from 12 orders of magnitude for mercury to 
7 orders of magnitude for Pluto) being, for most purposes, negligible. We also note that the 
quoted values for these corrections are well below the experimentally detected deviations from the 
theoretical predictions (usually obtained through numerical simulations) for the orbital motion in 
the solar system |51j. Hence, no bounds for the parameter ckq can be drawn from this analysis. 
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Planet 


GR precession J3.3H 


Ai GH 13.331 


if 2 precession J5.24I 




r\ 


orb/ri.orb < 5 . 1 1 1 




(arcsec/ century) 


(s / century) 


(arcsec/ century) 


(s / century) 





— ^ / century) 


Mercury 


10.35 


5.66 X 10~ 2 


e [-6.04, -0.122] x 10 -14 


e [1.27, 1.30] X 10 -16 


e 


-0.0706, 3.80] X 10 -32 


Venus 


5.30 


8.01 X 10~ 8 


S [-3.63, -0.0906] X 10~ 13 


= 1.93 X 10~ 15 


e 


-0.0460, 2.48] X 10~ 31 


Earth 


3.84 


3.64 x 10~ 6 


G [-9.59, -0.239] X 10 -13 


= 8.29 X 10~ 15 


e 


-0.122, 6.54] x 10~ 31 


Mars 


2.54 


4.48 X 10 -3 


S [-3.44, -0.0809] X 10 -12 


= 5.62 X 10~ 14 


e 


-0.043, 2.31] X 10~ 30 


Jupiter 


0.74 


6.12 x 10~ 4 


6 [-1.36, -0.0331] X 10~ 10 


= 1.39 X 10 -11 


e 


-0.171, 9.22] x 10~ 29 


Saturn 


0.40 


1.48 X 10 -3 


S [-8.48, -0.206] X 10 -10 


= 2.18 X 10~ 10 


e 


-0.107, 5.76] X 10~ 28 


Uranos 


0.20 


8.97 x 10~ 4 


£ [-6.81, -0.167] x 10~ 9 


= 4.96 X 10~ 9 


e 


-0.0861, 4.63] X 10~ 27 


Neptune 


0.13 


4.18 X 10~ 6 


S [-2.60, -0.0650] X 10~ 8 


= 3.70 X 10~ s 


e 


-0.0330, 1.77] x 10~ 26 


Pluto 


0.10 


1.28 


6 [-6.67, -0.125] X 10~ 8 


£ [1.44, 1.52] X 10~ 7 


e 


-0.0749, 4.03] X 10~ 26 



Table 9: Orbital precession and period corrections due to General Relativity (not including the effects of the 
Sun quadrupole J2) as given by equations iS.31]) and 13.33]) and due to the expanding background as given 
by equations ft5.24\ ) and &5.25\) . It was considered the range for the parameter ao £ [3, ao. max. pert] H5.15\) 
for which the perturbative expansion H5.16\) is valid, when no interval for the values of the correction is 
given is due to the respective upper and lower values differing by less than two precision digits. The allowed 
ranges for orbital radius time variation ri.orb/^i.orb os given by 15.11]) are also quoted for the circular orbits 
approximation and the full range of the parameter ao £ [3, +00 [. 



6 Conclusions 

6.1 Resume of Results 

In this work we have built an ansatz for a locally anisotropic metric ()4.4p describing matter in an 
globally expanding background, our universe. This metric interpolates between the Schwarzschild 
metric (13, 2p near massive bodies and the FLRW metric (I2,28h at spatial infinity where spatial 
isotropy is retrieved. By considering a two parameter radial coordinate dependent exponent on 
the shift function, a = a Q + a\ 2GM/(c 2 n) (fOUl) . we have maintain space-time complete, free 
of singularities except for the Schwarzschild mass pole at the origin and the total mass, within a 
shell of finite radius, finite. The negative parameter a\ < plays the role of a regulator removing 
the essential singularities at the center of mass. Its upper bounds have been analyzed in order 
to maintain the mass-energy density positive, however we note that it can, for most purposes, be 
considered as close to zero as desired being negligible outside the event horizon r\ > 2GM/(c 2 r\). 
As for the parameter «o > 3 fine-tunes the transition between small and large scale physics. For 
«o > 5 space-time is asymptotically flat near the event horizon and for relatively larger values we 
have shown that, although for relatively small spatial scales space-time can be contracting, for large 
spatial scales it will always be expanding consistently with the global expansion assumption. In 
particular this characteristic may contribute to dark matter effects at galactic scales maintaining 
both the usual physics of planetary systems (small scales) and global expansion (large scale). 
Hence we may have manage to solve the long standing puzzle of consistently describing matter in 
a expanding background first approached by McVittie. 

We have analyzed in some detail the orbital motion in the solar system considering the locally 
anisotropic metric H4.4j) and have concluded that the corrections due to the expanding background 
are negligible by many orders of magnitude with respect to the General Relativity corrections as 
well as to the detected experimental deviations such that no bounds for the parameter ao can be 
obtained based in experimental observations on the solar system. This result is welcome in the 
sense that the effects of expansion are negligible for small scale systems such that the standard 
General Relativity in Minkowski background describes these systems to a very good accuracy. Also 
we note that the anomalous Pioneer acceleration [40] towards the Sun of a p ~ 8 x lCP 10 ms -2 
cannot possible be explained by the expanding background effects. In the range of the Pioneer 
distance from the Sun n ~ [2.85, 6.60] x 10 12 m the maximum background correction to the General 
Relativity Newton law towards the Sun corresponding to a value of the exponent parameter ao ~ 
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10 10 is of order f\ ~ 10 25 ms 2 , hence 15 orders of magnitude below the measured acceleration. 

As a final remark it is relevant to note that we have at most considered a test mass and a 
central mass when deriving the analytical results in this work. When considering many body 
interactions the net effect of the corrections due to the expanding background may become slightly 
more significant, mostly for distant astrophysical body interactions (the expansion effects are 
proportional to the distance). In order to implement such corrections in numerical simulations 
it is necessary to derive the modified Newton law in isotropic coordinates for inclusion of the 
corrections in the PPN formalism [32] which is the most widely employed description of many 
body gravitational systems. 

6.2 Outlook: Possible Contributions to Dark Matter Effects 

The corrections to the General Relativity Newton law (|5.ip due to the locally anisotropic met- 
ric (|4.4p discussed in section [5] increase the orbital velocities for galactic scales. This effect is 
similar to the ones attributed either to the dark matter hypothesis [29J or modified theories of 
gravity [271 [28] . In either of these cases the modified Newton law at galaxy scales may explain 
the observed deviation from General Relativity predictions such as the flattening of the galaxies 
rotation curves [23j and the deviation from the predicted gravitational lensing by astrophysical 
objects [241126]. 

As discussed in section [5] the locally anisotropic metric (|4.4p , for large values of the parameter ao , 
also predicts a significant deviation of the orbital velocities for galaxy size scales as exemplified 
in figure [7] for a central point-mass with the value of the core mass of the large galaxy UGC2885. 
As can readily be verified, when considering the full mass of the galaxy to be point-like the radial 
distance at which this effect is verified is too large to significantly change the orbital velocity 
within the radius of the galaxy. However by comparison with the same effect for the Sun (see 
the modified Newton acceleration in figure [6] for comparison) we conclude that when considering 
several massive objects with a large parameter ao within the galaxy we may expect that the 
modifications to the orbital velocity will be felt for lower values of the radial distance, hence 
contributing to the dark matter effects. With respect to large scale observations we note that 
the local anisotropy detected in the background radiation [21J may help to set bounds on the 
parameter ao allowing to properly evaluate to which extend can the local anisotropy contribute to 
the astrophysical dark matter effects. We will give an account of these issues somewhere else [45] . 

As a relevant final remark we note that the exponent a has been introduced phenomenological 
when building the metric ansatz (|4.4j) without any more fundamental reasoning concerning its 
meaning. In particular if it represents some sort of unknown gravitational interaction associated 
to some sort of matter or if it is some kind of fundamental constant identical for all masses. If it 
is to be interpreted as an effective description of dark matter we may expect that for each body 
the parameter ao has a distinct value. For each galaxy (or set of bodies within each galaxy) 
its value can be inferred from the deviation of the usual General Relativity Newton law and 
fitted accordingly. Although our construction can be interpreted as that the locally anisotropic 
metric (|4.4|) represents the net effect of long range gravitational interaction in local masses due 
to the expansion of the universe, the specific nature of the interaction that it may describe or its 
meaning in terms of the known physical interactions (or fundamental principles) is not obvious 
to the author and is here left as an unsolved question. We note however that a dilaton-like scalar 
field [331 EH EH] with local solutions of the form $ = a log(l - Use) such that (1 - U SC ) a = 
e a lo s( 1_f/ sc) — e * could justify this metric. The author is unaware of a specific theory that 
predicts these particular solutions. 



75 



Appendixes 



A Conventions and Definitions 



In this appendix we list the General Relativity and Special Relativity conventions employed in 
this work. For units we use the International System maintaining all the physical dimensionfull 
constants explicitly in the equations. The Einstein convention is considered such that, unless 
explicitly stated, sums are implied when repeated indexes are present with the following convention 
for greek and roman indexes 

H = 0,i , i = 1,2,3 . (A.l) 

For the definition of Cartesian coordinate parameterization we consider x° = ct, x 1 = x, x 2 = y 
and x 3 = z in an infinite interval. As for the definitions of the spherical coordinate parameteriza- 
tion and its mapping to Cartesian coordinates we consider 

*/ oc^ i y^ y 

x° = ct , x 1 = r = J x 1 + y 2 + z 2 , x 2 = 9 = arctan — , x 3 = <p = arctan — , (A. 2) 

z x 

with 9 G [— 7r/2, 7r/2[, <p £ [0, 27r[ and r G [0, +oo[. The inverse mapping between spherical and 
Cartesian coordinates is 

x = r sin 9 cos ip , y = r sin 9 sin (p , z = r cos 9 . (A. 3) 

Bold face quantities r (x) represent spatial vectors for spherical coordinates (Cartesian coordi- 
nates) and, when not specified, a dotted function / represents simple derivation with respect to t 
and a primed function /' with respect to r. 

For a generic metric g^ v we have the following definition of the infinitesimal length square 

ds 2 = c 2 r/oodT 2 = g iiv dx il dx v , (A. 4) 

being r the proper time corresponding to the free falling frame and we adopt the metric signature 
(+,—,—,—). The connections are defined as 

Tfl uS = \d W ( 9pv,5 + 9p8,u ~ 9v8,p) , (A.5) 

where the index ', 5' denotes simple derivation with respect to x s and the covariant derivative for 
a vector is defined as 

Dtp 

dx v 

We note that, considering non-natural units (c^ 1), for derivatives with respect to the coordinate 
time t = x°/c we obtain the expression q^. t = q^ t + cT^ 0S q s . 

The Riemann tensor is defined as 

t>P -pp -pp I -pp -pS -pP -pS { \ 7\ 

U u\p — L up,X 1 u\,p~t~ L S\ L up 1 S p L u\ > 

and, unless otherwise stated, we express the Einstein equations in the following form 

G - 

^ pv — ^4 ± jxv •> 

Gfiu = Rpu-^gpuR, ( A - 8 ) 



i% = -iSr = 1% + r V 5 • (A.6) 
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where G^ v is the Einstein tensor, T^ v is the stress-energy tensor, R^ v = R a mu is the Ricci tensor 
and R = g Sp R a g ap is the Ricci scalar. In the last line we consider the stress-energy tensor for 
a perfect isotropic fluid where p is the pressure, p the mass-density and u p = (c,Vi) the fluid 
velocities. In the commoving frame of the fluid the four- velocity is given by = (c, 0) and we 
obtain the tensor components 

r c om = pc 2 +p(i - 500) , nr = -pgoi , ^ om = -pw , (a.9) 

The sign conventions in these expressions depend on the metric signature choice. We note that 
under a coordinate transformation = x^{x v ) the contravariant velocity components generally 
change to v u = dx^ jdx^v^, hence as long as we consider transformations for which x° = x° we 
remain in the commoving frame of the fluid and equations (|A.9|) are still valid. 

For light-like trajectories we compute the geodesies directly from the distance element (|A.4j) . 
ds 2 = 0. As for time-like trajectories of massive objects we consider the usual geodesic obtained 
from the Euler-Lagrange equations that minimize the length L = V ds 2 . The solution for the 
relativistic factor dt/dr = 7 is obtained either from integrating the geodesic equation for x° or 
directly from (|A.4|) 



2 f dt\ 2 dx % f dt\ 2 dx % dxi ( dt\ 2 



7 



dt / rjoo 



(A.10) 



dr \ n 1 Zgoi* 1 1 9a&& 



900 + ^ + 



c 



Then the geodesic equations for the remaining space coordinates are 

d 2 x % , dx^ dx v 

~d^ = dr dr ' 

a? = -c 2 r i 00 -2cr 0j i; J '-r ifc i J x A: -7- 1 7^ 



(A.ll) 



For projections of generic vectors from the four-dimensional space-time manifold to the spatial 
hyper-surface at a given fixed time we take the following definitions [39] for the intrinsic metric 
^gij and three-dimensional vector components 

(3) <^ = - (4) ^-, 

(A.12) 

(3) w i = (3) g ij (4) w . = (3) glj ^) g .. (A) w i + (4) 5 . Q (4)^0) _ 

In this equation ^'g^v stands for the usual four-dimensional metric components and stand 
for the four-dimensional vector to be projected. In particular this projection is applicable to the 
space-time coordinates ^x^. 

In this work, both in order to simplify the technical details and for easier physical interpreta- 
tion of the derived quantities we consider several coordinate systems corresponding, for spherical 
coordinates, to the following radial coordinate definitions labeled by an underscore index as: 



r - a radial coordinate with integration measure yf—g = a 3 r 2 sin#, and a time dependent 
area for the 2-sphere, A(t) = Ana 2 r 2 , being a = a(t) a generic cosmological solution for the 
universe scale factor. In this work we refer to it as the non-expanding coordinates due 
to the geometrical spatial lengths expanding over time while coordinates lengths are fixed. 
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T\ — a radial coordinate with integration measure \J— g = r\ sin 9, hence with a time-independent 
area for the 2-sphere, A = An r\. These are the usual Schwarzschild coordinates and we will 
refer to it as expanding coordinates due to the coordinates lengths being coincident with 
the physical geometrical lengths, hence expanding over time. r\ is related to r as 

n = ar . (A. 13) 



T2 - the usual expanding isotropic coordinates employed for the Schwarzschild metric for 
which the speed of light is the same along all space directions, r 2 is related to r\ as 

/ GM \ 2 f 2GM fk , 

n = r 2 1 + k for n > — 5— . (A. 14) 

\ 2r2 c z J & 

r3 - non-expanding isotropic coordinates for which the metric is explicitly dependent on 
the scale factor a. r 3 is related to r 2 as 

r 2 = ar 3 . (A. 15) 



In table [TH] are summarized the IS units for the quantities and the values of the fundamental 
constants employed in this work. 

symbol units value name 



c ms- 1 2.998 x 10 8 

G Kg~ 1 m 3 s~ 2 6.670 x 10" 1 
M Kg 
m 

R mT 2 



speed of light 
Newton Constant 
mass 

space-time coordinates 

metric tensor 

metric connections 

Ricci scalar (coordinate curvature) 



Table 10: IS units for the fundamental constants, metric and other derived quantities 



B Cosmological Metric 

In this appendix we summarize technical details for the FLRW metric representing a homogeneous 
and isotropic, globally flat, expanding universe. 



B.l Non-Expanding Coordinates 

In the commoving frame for non-expanding coordinates this metric is given by equation (|2.2p and 
the non-null connections are, for Cartesian coordinates, 

r % = ^^,r 0j = rv = ^^. (B.i) 
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The Ricci scalar (scalar curvature) for this metric is 



%].inv -0 I — ) -6-^ , (B.2) 



and the non-null Einstein tensor components are 

Goo = 3 (-^j , Gij = -S i:i i (a 2 + 2aa) . (B.3) 

B.2 Expanding Coordinates 

For expanding coordinates this metric is given by equation (|2,28p . for Cartesian coordinates the 
non-null connections are 

r° - ^■■T i T j ( — \ r° r° - -A-V f 

1 oo — o tJ x 1 x 1 I I ) 1 oi — L io — ^jn I 

r%- = dyA , r* = x < s jk 4x1 (±Y- X ^, (B .4) 



and, for spherical coordinates are 



3 



bo " r H^J ' r ° 01 " r ° 10 " ~ n w 

1 11 — — ! 1 022 — r l L 11 

r» 22 . r ? si „^r» 33 , i*„ = -E£ + r } (-£.•* 

r oi = t\ = — r° 00 , r 1 11 = — r° 01 (b.5) 

r i 22 = _ ri+r 3^ ( r i 33 = sin 2 ^ 

r 2 12 = T 2 21 = - , r 2 33 = — cos 9 sin 6 , 

r 3 _ r3 — r2 r3 — r 3 — cos ^ 

1 13 — 1 31 _ 1 12 j 1 23 — 32 — sin 6* ' 
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The Ricci scalar (scalar curvature) is given by the same expression than for expanding coordi- 
nates (|B.2p and the non-null Einstein tensor components for spherical coordinates are 




(B.6) 



B.3 Short-Scale Time Evolution of the Expansion Factor 

When the time to considered for the series expansion of the universe scale factor a given in 
equation (12. 5|) does not coincide with the time i eX p corresponding to the experimental measurement 
of the Hubble rate H, 



exp 



a/a\ 



t=te 



and the deceleration parameter —q exp /H l 



2 

cxp 



-a/a\ t __ 



there is a small correction for these quantities when evaluated at io, Hq = a/a\ t=tQ and —qo/H^ = 
— a/a\ t=t . Similarly to the expansion (|2.5p let us take a series expansions to third order on the 
product of the constant Hubble rate by time (Hot) 3 such that in the neighborhood of t eX p 



a cxp 1 + 



A exp 



1 a 



1 a, 



(i iexp) „ (i iexp) ~i~ „ (t iexp) 



•"-exp 



2 a, 



cxp 



6 a. 



cxp 



A cxp 



^1 ~\~ H cx p(t tcxp) ^QexpHgyrp^t iexp) 



+0 ((t - iexp) 4 ) 



Q S expQexpH ex p(t iexp) 

+0 (^ 4 (i - i exp ) 4 ) 



(B.7) 



where in the last line we replace the derivatives of the scale factor by the Hubble rate H (II. ip . 
deceleration parameter q (|1.2p and variation of the deceleration parameter s (|1.3p evaluated at 
i exp . Hence expanding the Hubble rate H and deceleration parameter q in the neighborhood of 
iexp we obtain, for their value at the reference time i = io the following first order (on time) 
relations 



Hn 



H, 



exp 



1 — (<Zexp + 1) -f^cxp (io — iexp)) ; 



t=t 
d 1 

a IP 



(B. 



g exp (1 + (1 + 2q, 



exp S e xp) H ex p (ig 



iexp ) ) 



t=t 



Without loss of generality the above series expansions can also be employed with a negative value 
for i ex p- For the specific case of expanding coordinates x\ the reference time is set to zero, io = 0. 

C McVittie Metric 



The McVittie metric for a point-like central mass M at the origin of the coordinate frame in an 
expanding background was originally derived for non-expanding isotropic spherical coordinates 
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corresponding to r% as denned in equation (|A.15P 

GM \ 2 



ds 2 = r 2 ( 1 2 ~^ 



1 + 



GM 
2 a r3 c 2 



di 2 - a 2 I 1 + 



GM 
2 a r 3 c 2 



(dr 2 + r 2 d0 2 + r 2 sin 2 # V) • (CI) 



In the following our aim is to express this metric in the several coordinate choices corresponding 
to non-expanding coordinates r, expanding coordinates t\ and expanding isotropic coordinates r 2 
as defined in appendix [Aj Hence considering the coordinate transformation (|A.15P we obtain 

a 



T3 



ds 2 



1 



GM 
2r 2 c 2 



1 _l GM 



1 + 



2r 2 c 2 

GM 



c 2 a \ 2 T2 c 2 



, - 2r 2 a ( GM 
d t 2 + _^ _ l + 

c a V 2 r2 c z 



dtdr 2 (C.2) 



(dr| + r^0 2 + r\ sin 2 ^ 2 ) . 



2r2C 2 r 

The transformation (|A.14p is, generally, not bijective. Then, instead of proceeding directly from 
the above metric, we can directly infer the expression for the metric in the coordinates r\ 

2 



T2 



ds 2 



n i + 



GM 
2n c 2 



r\c z & \a ) } a . I \ _ 2GM 

ri c 2 



dt df\ 



dr\ 



Y _ 2 GM 



"\ [dO 2 + sin 2 Odif 2 ) . 



Finally considering the coordinate transformation (|A.13P we obtain 
r\ = ar , 



(C.3) 



ds 2 = c 2 I 1 



2 GM a 2 r 2 
arc 2 c 2 



+2 aar 



a 2 dr 2 



2 gm 1 



2 




2GM 
ar c 2 


i > 




2GM 
arc 1 ) 



+ 



1 



2GM 



,, . u - a 2 r 2 (dO 2 + sin 2 ^ 2 ) . 

1 — 



(C.4) 



We note that the transformation to isotropic coordinates is valid for large radial coordinate as 
expressed in (IA.14p . this may be the reason why McVittie originally did not notice that the 
Schwarzschild radius is singular. 



D Locally Anisotropic Metric 

Except for the curvature invariant lZ a , in the following expressions we are omitting the Schwarzschild 
mass pole contribution which, when required, is considered in the main text. 
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D.l Constant Exponent a 

The non-null connections (a\£ for the metric (|4,4p with a constant exponent a in the shift 
function (|4.3j) expressed in expanding spherical coordinates T\ are, for (3 = 0, 

1 j uu c 2 r! l^ac/ \ c 2 ri J 



c 3 \ a / \ c 2 ri / \ c 2 ri 



(a) 1 01 — (a) 1 10 



GM / 2GM\ _1 



c 2 r\ \ c 2 ri / 



acj \ c T\ I \ cr n 



/aW 2GM\f-| / (a-2)GM \ 



(a) 1 22 



a 1 



r° _ rU^) (i 2GM ^ 



ac J \ c 2 n y 



(a) r °33 = sin26, (a) r °22 > 



(D.l) 
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for 13 = 1, 



! GM f 2GM\ f a \ f 2GM\^ 



a i 1 / a 1 



+ri | A) 2 A _ ™<) i+i f , _ (■ , _ (, + M> GM ) 

acj \ cr ri / \ \ cr ri / \ c z ri / 



+r? ^— Ml- — I ( U 
pi _ ri 

(a) 1 01 — (a) 1 10 



a \ 4 ( 2GM\ 2a ' 1 ( (a-2)GM 



\ ac) \ c 2 r\ ) \ ~* (? r\ 



a 1 



GM ( a \ / 2GM\ 2-2 



c 2 ri V a c 7 ^ c 2 r\ ) 

3a 3 



•:(^)('™) T "'('^)- 



r i GM ( 2GMM 1 

(a) 1 ii - ^ 1 J 



„. / 2GM\ , / a Y f 2GMY" 

(a) rl 33 = s i n2 ^(a) rl 22 , 



(D.2) 

and the remaining connections coincide with the ones for spherical coordinates in flat space-time 

(a) r2 12 = (a)T 2 2i = — , (a)T 2 33 = COS 9 shl 6 , 

(D.3) 

F 3 r 3 _ 1 r 3 _ r 3 _ COS^ 

H 1 13 - (a) 1 31 - — > (a) 1 23 ~ (a) 1 32 ~ • 

In these expressions we have ordered the terms in increasing powers of the Hubble rate (a/a = H 
and a/a = —qoH 2 ). 

The Ricci scalar (curvature) is 

R a = -6R T (— ) - 6 R n (D.4) 
\ac J \ac z J 
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where the coefficients Ri and Ru are functions of r\ only 
Ri ■■ 



2o ( _ 2C/.U V' 2 / _ (7-a) GM \ 



c 2 ri \ c 2 r\ J \ 3c 2 ri / 



a 3 



i _ 2GM ^ "2 ^ _ (6-a)GM j + 2 ^1 - 2GM \ 



c 2 ri 



/ \ 3c 2 r\ ) V c 2 ri / 



(D.5) 



a 3 



\ cr r\ J \ 6c z r\ J 

where we organize the several terms in growing powers of (1 — 2GM / '(c 2 n)). The non-null 
components of the Einstein tensor are 
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3a 1 



2(;.i/\ _r " 



a 

ac J \ r\ 



a 1 



2GM\ / 2(a - 3) GM\ 



x 2-3 1-^— 1 + 



-2r 2 I — | I 1 I 1 - 



a \ 2 f d \ ft 2GM \^~ ■ 



c 2 r\ ) \ 3c 2 r\ / 



3a 1 



1 \acj \ac 2 J \ c 2 r\ J 
r {a) _ 901 r (a) _ "I/, 2GM\t+f (a) 

Lr 01 — <J n — I 1 5 ] Wl ' 



a 3 



a 



ac 2 J \ c 2 r\ 



1 



c 2 n/ (D.6) 



7 



4(a - 3) GM 2(a 2 - 5a + 6) (GM) 2 



x 3 + ^ + 



c 2 ri c 4 r 2 
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The curvature invariant is 



n,+n u [—) +n in (-%\+n IV (—\ + n v ^ I + n vl \ — 



ac I \ac 2 ) V a c I a 3 c 4 Vac 2 



where the several coefficients TZj to TZyi are 
48(GM) 2 



c 4 r 6 > 

16a (GM) 2 f 2GM\%-% ( 2(a-l)GM f 2GM\^ - ^ 
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c 4 r\ \ c 2 r\ ) c 2 r\ \ c 2 r\ J 

2GM\°~ 1 16a 2 (GM) 2 / 2GM\^~^ 



(a - 9) GM 2(a - 5) (GM) 5 



x 2 + 



c 2 n c 4 r 2 
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and 



8a 2 (GM) 2 ( 2GM\ a ~ 3 16a GM ( 2GM\ a ~ 2 



Kv = W^" 1 



c 4 r? V c 2 r\ 



) c 2 r\ V ° 2 r i / 



2GM\ a ~ 1 t 16a 2 (GM) 2 f A 2GM\^~* ( n , (a - 5)GM 



-24 ( 1 - -g— + j-f-L- 1 - 2 + 



+ 



8n GM ( x _ 2_GM\ ¥"! / _ 12GM\ L 2_GM\ " (D . 9) 



\ ' > - 5 2 +^41- -2 

c z n \ c z ri J \ C A T\ J \ C z T\ J 



4a 2 (GAf) 2 / 2GM\ 



c 4 r 2 c 2 ri ) 



8a GM ( 2GM\ a ~ 2 ( 2GM\ a ~ x 
c z n \ c z ri / \ cr ri / 

where we organized the several terms in growing powers of the factor (1 — 2GM/(c 2 ri)). 



D.2 Space Dependent Exponent a(ri) 

When the exponent a in the shift function f|4.3|) of the metric ansatz (|4.4p is dependent on the 
radial coordinate ri, a = a(r\), as introduced in equation (|4.20j) . the connections are 

(a(n)) r % =(«) + Ar^ , (D.10) 

where (^r^, are the connections for constant exponent a given in the above equations (|D.1|) . (|D-2p 
and (|D.3P and the non-null corrections Ar^ are expressed in terms of the derivative a' = 
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da[ri)/dr\ as 



r?/a\ 3 / 2(?M\t4 / 2GM\ , 
iii iii a 



Ar° 00 = +± - i-^ log i 



/l 2 / & I 2 / 
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aC J \ C T\ J \ C T\ J 

Arl rl(a\ 2 (^ 2GM\ a+ \ / 2GM\ , 



+T ( — ) ( 1 " ^F 1 ) lo § ( 1 - ^ ) «' 

2 \acj \ c z ri / \ c A r\ J 



4 / n/1I ,\2Q 



2GM\ , / 2GM\ 



ArV = Ar\ 



3a 1 



~ I — ) M 

ac , 



V \ c z n / \ crr\ J 



The curvature is 



R„ ( ^\ = -6(Ri + ARj) -6(R u + AR n )(-^) , (D.12) 



l a(ri) 



where the coefficients i?/ and -R77 are given in equation (|D,5P and the corrections due to the radial 
coordinate dependence of the exponent a are expressed in terms of the derivatives a' = da/dri 
and a" = d 2 ajdr\ as 

2GM ( 2GM\ a " 1 ( , / 2GM\\ , 

a 1 



i (\ 2GM\-" I (\ 2GM\ [ \ } f t 2GM\ , 
2 / 2GM\ a , / 2GM\ 



1 \6\ c 2 r\) 3 \ c 2 ri j J \ c 2 r\ J 



) iog(i-^ )(M 2 + («") 2 ) 



6 V. c 2 ri 



ri A 2GM\t"5 / 2GM\ . 



(D.13) 



The Einstein tensor is 

G $r 1 )) =G («) +AG (D . 14) 
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where G^t) is the Einstein tensor for constant a given in (ID. 61) and the corrections AG^ V due to 
the dependence on the radial coordinate of the exponent a are 

f a\ 2 f, 2GM\ a+ \ / 2GM\ , 
AGoo = n — 1 - — — log 1 - a' 
\acj \ cr r\ J \ c z ri J 



Vac/ V c z ri/ V c z ri / 



an ac V c z n J 



( a \ 2 / 2GM\ a_1 1 / 2GM\ , 
AGn = -n — 1 - — log 1 - — — a' , 
\ac J \ c z r\ J \ c z n / 

r\(a\ 2 (±GM( 2GM\ a ~ 1 { ( 2GM\\ . (D.15) 



y c 2 ri \ c 2 n / V ^ V c 2 r\ ) 



2GM\2^, / 2GM\ . / 2GMY 
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c A T\ ) \ c z r\ J \ c z r\ J 



, . , 2GM\ / 6 , / 2GM\ , „ 2 , 

x log ( 1 - -= —a' + log 1 - -= — (a') + a' 

c ri / \ri \ c z ri / 

r? / a \ / 2GM\t-i / 2GM\ . 
1 ' w 1 ' log 1 - -= a' , 



2 \ a c 2 / \ c 2 r\ J \ c 2 r\ J 

AG 33 = sin 2 9 AG22 • 
The curvature invariant is 

= ^ + (^w + A72. w ) (— ) + (7J m + AK IH ) (-^ 

+{K IV + AK IV ) (—) + (K V + AK v )^j + (K VI + AK V i) (A, 
\acj a? cr \ a c z 

(D.16) 

where the coefficients IZi to TZyi correspond to the coefficients of the curvature invariant (|D.7p 
for constant a given in equation (|D,8P and (|D.9|) . The corrections to these coefficients ATZi due 
to the radial coordinate dependence of the exponent a are 



8GM / AGM ( 2GM\ a ~ l ( , / 2GM\\ , 
AR U = -5-5- — 2 1 " — 2 — 1 + alog 1 - — — ) a' 



9 2 
c z rf 
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x log 1 1 f 
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(D.17) 
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and 
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where the coefficients AlZiv-i to ATZjy-4 in the definition of ATZjv are 
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D.3 Planetary Orbits 



In this section we derive the perturbative differential equation describing orbits of a test mass in 
the gravitational field of a central mass M for the locally anisotropic metric (|4.4p with a radial 
coordinate dependent exponent a(r±) = ao + ai Use (|4.20p . 

We take the same approach of section [3721 starting from equation (|3. 14[) . For the specific case of the 
locally anisotropic metric (|4.4p and considering an orbit lying in the plane of constant coordinate 
8 = tt/2 such that d6 = and sm.6 = 1, we obtain the following equation 



c 2 = (l-t/sc-f^Va-t/scr) M) 2 



Q 1 



v,2 



(D.20) 



+2 ( tl^L ) (1 - C7 sc )f-f cih - i _ 1 jJsc - rfip 2 



where we have replaced the derivatives of the scale factor by the Hubble rate H = a/a (II. ip . 
have written the equations in terms of the Schwarzschild gravitational potential Use ()3.ip and the 
dotted quantities represent derivatives with respect to the proper time r. We will further take the 
approximation of static orbits considering that the Hubble rate is a constant corresponding to the 
measurement of this rate at the reference time to, Hq = H(t$) (jl.ip . In this way the Lagrangian is 
independent of the time coordinate and a conserved constant of motion corresponding to energy 
exists given by the functional variation of the Lagrangian with respect to ct 

1 5C _ 2E H 
m S(ci) mc 



2 1 1 / I^Yo-f/scn (ct) • 2f^i) : . n 



a 1 



This equation can be solved for ct 



(D.21) 



. Em. - (Mori) n - u s C )f-2 fi 
ct = ^ V c M , (D.22) 

l-^sc-(^) (I- U sc ) a 



such that replacing this solution in the Lagrangian (|D.20p we obtain 

i-u sc -(^) 2 (i-u S cT 1 ~ u ^ 



o \ mc I V c / V w 1 Ti n r, 

c = ~ 77 2 q 77 r lV ■ ( D - 23 ) 



Multiplying by the factor 1 — Use — ( H °J 1 ) 2 (1 ~~ ^sc) a and gathering the constant terms in the 
left-hand side of the equation we obtain 




2GM o , , / 2GM\ 
r 2 - rfc 2 1 - 



(D.24) 

2 / o^«)i-r\ a 



{—) +{nHo) V'^J 



The terms in the first line of this equation match the usual General Relativity terms, while the 
terms in the second line are the lower order corrections (in Hq) due to the expanding background. 
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In addition the Lagrangian is independent of the coordinate ip such that a constant of motion 
corresponding to angular momentum exists given by the variational derivation of the Lagrangian 
with respect to up 

1 SIL 

2J = -2r\ub . (D.25) 



m Sip 

We note that this expression matches the same quantity for Keplerian and General Relativity 
orbits (|3.16p . Further re-expressing the time derivatives r\ as derivatives with respect to ip, 
dri/dr = dr\/dip x dip/dr and considering the change of variables u = 1/r as given in equa- 
tion (|3.17|) , replacing the respective derivatives in equation (|D.24p , differentiating it and factoring 
out an overall factor of 2u' J 2 we obtain the lower order equation in Hq for the function u(cp) 
which describes an orbiting test mass in the gravitational field of a point-like central mass M 
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u log 1 
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(D.26) 



The terms in the first line match the usual terms in the respective General Relativity differential 
equation (|3.19j) while the terms in the second and third lines are the corrections due to the 
expanding background. 



r l orb (Xl0 a to) 



Mercury Venus 



57.91 
0.2056 



108.21 
0.0067 



149.60 
0.0167 



227.92 
0.0935 



Jupiter 



778.57 
0.0489 



Uranos Neptune 



1433.53 
0.0565 



2872.46 
0.0457 



4495.06 
0.0113 



5906.38 
0.2488 



Table 11: Planetary orbits parameters considered: the semi-major axis ri. or b and the eccentricity e \52\ 



In table[TT] are listed the orbital parameters for the planets considered in section [5] when estimating 
the corrections due to General Relativity and due to the expanding background listed in table 
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